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Aus gesundheitlichen Gründen hat der bisherige Herausgeber dieser Zeitschrift, 
Prof. Dr. Dr.h.c. F. A. Willers, seine Tätigkeit niederlegen müssen. Seinem Wunsche folgend 
habe ich nach Vereinbarung mit dem Akademie-Verlag Berlin und mit Zustimmung des 
übrigen Herausgeberkreises und der führenden Fachkollegen der Gesellschaftfür Angewandte 
Mathematik und Mechanik mit Beginn des neuen Jahrganges die Geschäfte des Heraus- 


gebers und des Schriftleiters übernommen. 


Ich glaube, mich zum Sprecher aller Leser und aller Mitarbeiter der Zeitschrift für An- 
gewandte Mathematik und Mechanik zu machen, wenn ich zu Beginn meiner Tätigkeit 
Herrn Willers für seine aufopferungsvolle Arbeit und alle Sorgfalt, die er zweiundzwanzig | 
Jahre lang, zunächst als Mitarbeiter von Trefftz und nach dessen Tode als sein Nachfolger, 
der Zeitschrift hat angedeihen lassen, den Dank ausspreche. Wenige Jahre, nachdem er 
1937 die Schriftleitung übernommen hatte, stand er vor der schweren Aufgabe, die Zeit- 
schrift auch während des Krieges erscheinen zu lassen, sie auf ihrem hohen Niveau zu hal- 
ten und sie über den Krieg hinweg zu retten. Dies alles ist ihm in vorbildlicher Weise ge- 
lungen, und 1947 war die Zeitschrift für Angewandte Mathematik und Mechanik die erste 
mathematische Zeitschrift in Deutschland, die — von ihm zu neuem Leben erweckt — 
wieder erschien. Seiner umfassenden Sachkenntnis, der Leidenschaft, mit der er sich der 
Arbeit an ‚seiner‘ Zeitschrift hingab, und dem hohen Ansehen, das er bei den Fach- 
genossen im In- und Ausland genießt, ist es zu danken, wenn sich die Auflage seither fast 
verdoppelt hat und die „ZAMM“ zu einem weltweiten Begriff geworden ist. 

Als neuer Herausgeber sehe ich es als meine vornehmste Pflicht an, ihm nachzueifern, 
die Zeitschrift auf ihrer wissenschaftlichen Höhe zu halten, sie ebenso vielseitig zu gestalten 
und sie nach den Richtlinien zu leiten, die schon v. Mises bei ihrer Gründung aufgestellt 
hat und die auch heute noch ihre Gültigkeit haben. Ich bin mir bewußt, daß ich dazu der 
Mithilfe der Fachkollegen bedarf, und ich bitte alle Mitarbeiter und Leser um ihre tat- 


kräftige Unterstützung. 


Dresden, den 1. 1. 1959 H. Heinrich 


2 ZAMM 39 (1959) Heft 1/2, Seite 2—8 


Ein durch allgemeine Massenkräfte beanspruchtes 
unendliches elastisches Medium 


Von Vaäclav Vodicka 


Der folgende Aufsatz bringt die vollständige Lösung des Spannungsproblems für ein unendliches Medium, 
das durch beliebig verteilte Volumkräfte beansprucht wird. Der gelegentlich in der Literatur behandelte Son- 
derfall einer Binzellast erscheint im Rahmen unserer allgemeinen Theorie als eine ganz einfache Aufgabe. 


The paper presents a complete solution of the stress distribution problem in case of an infinite medium 
strained by arbitrary volume forces. The special case of a single load, occasionally considered in the literature 
on the subject, appears as a very simple problem when treated from the point of view of our general theory. 


L’article suivant donne une solution complete du problöme des tensions pour un milieu infini qui est 
influence par des forces de volume distribudes a volonte . Le cas special, traite de temps ü temps dans la 
litterature, d’une charge unique apparait dans le cadre de notre theorie generale comme un probleme tout ü 
fait simple. 


B HacToameü cTaTbe NaeTca UeJIHOE pellieHue 3anayu O HalpsIHKeHUH ÖeCKOHEYHOÄ CPenkl 
monBepraeMoi IIPOU3BeJIbHO pacnpenejIeHHbIM IIPOCTPAHCTBEHHEIM curam. Hacrusä cayyaü 
cocpeNoToYeHHOH HATPY3KM, HHOLNA BCTpeyamINHäcah B IepeNoBofA AurTeparype, B PamkKaxX 
Hamei OoÖMei TeopHu OKAasbIBAeTcA COBePpIIeHHO UPOCTOH 3apayei. 


Das in der Literatur [1] zu findende zweidimensionale Analogon unseres jetzigen Problems 
gibt gar keine Vorstellung von der mathematischen Kompliziertheit des hier (so allgemein hof- 
fentlich zum ersten Male) betrachteten dreidimensionalen Falles. 


I. Allgemeines Problem 


1. Es seien o die spezifische Masse und u die Poissonsche Zahl des homogenen und iso- 
tropen Materials, X = X(z, y, 2), Y=Y(«, y, 2), Z= Z(z, y, z) die auf die Masseneinheit bezogenen 
Komponenten der Volumkräfte. Der dadurch bewirkte Spannungszustand beschreibt sich dann 
durch folgende drei Gleichgewichts- und ebensoviel Kompatibilitätsgleichungen: 
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ar 02 ee De ta 
7 a 00, 

lt hezel ao <zuuszEEE 

0? 0° OT 
Not ufo, Er ne Spalt — HU — 
te] ale ne, + 2 Vol RR M). 
02 0 ru; 

Ay? lo; — (0, + m) ze Fr) lo, — u, + 0)]—2(1l+u) Öyöz =0, 
le tl + tat 
are rare rer alt d: 


— oo <y,y2<+tmw 


Das Gleichungssystem (1) ist natürlich unter geläufigen Voraussetzungen über das Ver- 
schwinden der gesuchten Spannungskomponenten und ihrer Ableitungen im Unendlichen zu 
lösen. Dies geschieht durch Anwendung der dreidimensionalen Fourier- Transformation. 


2. Die Beziehung zwischen einem bekannten Bedingungen unterworfenen Original 


h(&, y, 2) = Ft[h(&,n,&)] und seiner Fourier-Transformierten h(£,n,£) = F[h(x, y, 2)] drückt 
sich durch 


SR Ne ST | 
KEN.) = Flha&y, 2] = ( ar | h(x, y, z) ei Extny-+H2e) de dy dz | 
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aus. Dabei gilt das Faltungstheorem 
Di RE 1, 996, 1,8) en tetwstn dE an ac = SS] Io Br) 9a y—B,2— y) dxdßay 
= IS ey — 2 7)9@, B, 3) dx dB dy el (3) 


und eine Reihe von (leicht zu beweisenden) Rechenregeln, worunter für weitere Zwecke vor 
allem die Formeln 


| — _iEF[h@y, 2], Pi PR] -izFilh&n0] ... 


(und die sich daraus durch zyklische Buchstabenvertauschung ergebenden) von Wichtigkeit 
sind. Es sei noch bemerkt, daß man sich bei der Berechnung der Integrale (2) auf Cauchysche 
Hauptwerte beschränken kann. 


3. Durch Anwendung der Transformation (2) auf (1) ergibt sich unter gegebenen Ver- 


hältnissen folgendes Gleichungssystem für die Transformierten 0,, 0,, ... der gesuchten Spannungs- 
komponenten o,, Oys...} 


EOEENT, 00; =—ioX, 
0, 3 &tr, +ew,=—ioY,; 
9, +87, +0, =—ioZ, Br) 


we —- + —- 2, tr @® HM), +21 +W)iniy=d; 
Beer HE + ur —- dr, + red +2 ren. 
we, + EI), + - 2 rd, +21 +mM)ELT,—0 


Bildet man nun mit &,n die Linearkombination der zwei ersten Gleichungen (5), so ergibt 
sich leicht unter Zuhilfenahme der dritten Beziehung (5) die erste Gleichung des neuen Systems 


Eu +Py— 2, +27, =—ieEX+nY—LZ), 
-E, +2,42, +29, =—iel-EX+nY+Z2), RT 
eu, — 70, +2, +2lE7,=—ioleX—nY+LZ) 


Durch Hinzufügung der drei letzten Gleichungen (5) kommt man so zum äquivalenten Ersatz- 
system, woraus man leicht die Unbekannten 7,,, Ty., 7,, eliminieren kann. So kommt 


SEEN ee ER ENN LZ 

7, +5 —-(l+%)- e+n ,=—(+w)io Re "u, 
EEE SW LAT NITCh 
,+9,+5,—(l+ u) ara“ ,=-—(i+mie- ae a ie 6). 
Dr amt RE Z 
a o,=—(l+nio ara 


Nun sind aber die Gleichungen (6) äquivalent mit dem neuen System 


+ eree EX+nY+E£z 


u 4) ?+& 2 (1 nie — 2 2+& , e 
ö, 6, . [EX +nY—EZ —EX+nY+LZ 
e+r+le. ea)=iel rear; az), 
ER ER pre nn 
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Anlaß geben. Geht man damit in die erste Gleichung (7) ein, so kommt 


’ X % Z 1 e Eee 
a+mien Fr ROuD EN Eee Br 


mare BER tn FHtZ42@+ +9] 


—EX+4nYH4tZ , EX+4nYH+iZ 
=(1 — ai 1— v 
(-+mio re Regen — 1 HE 
Man hat also 

1 b PILZ 


= — — d a 
l—u A+&. | 


und die Formeln (8) erscheinen in der Gestalt 
= dSWE+rFR Rn ? RER re) z 
-- = Fararall L +2), 
Hr Era ze) Rz HER). Kö | 
gene. Fal- el een 2 
+ Bas) EX+n r| 


Unter Zuhilfenahme von (8) und (9) ergibt sich schrittweise aus der zweiten Gleichung (5.1) 


a EX+nY+LZ)+ 2, — 720, — 27, = _-ie(-2X+nY+LZ) 
N E = 
ter ptel—EX 4nYHIZ) PlEX—nY+LZ) 
-PEXRHHNY ZI HER +I-—PE+HD-CE+M1® 
io NINE ERS nEßB+ELY+nZ 
e+n7 al an % 2n? Y—2R28)=—2iont- BIT : 
und analoge Ausdrücke erhält man auch für 7,,, 7,,. Zusammenfassend hat man also 
N art rnZ = ETIEREFE 
e e Arm+a an 
ee nX+EY+£&n® RA : 
DUNSEE 0 
res 


wobei ® durch (9) erklärt ist. Damit bekommt man endlich 
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Die Formeln (10) und (12) bestimmen die Fourier-Transformierten der unbekannten Spannungen 
und unsere Aufgabe besteht nun in der Auffindung der zugehörigen Originalfunktionen. Natür- 
lich ist dies keineswegs leicht. 


4. Um den Übergang von (10) und (12) zur gesuchten Lösung o,, o,, ... vollführen zu 
können, brauchen wir die Originale zu folgenden Fourier-Transformierten; 


RS in 118 L&° un: 

E+P+E BEip+E Binp+e Eim+ Erprmp 
i® in ie? it 

e+T+ E+r+S Ctr+ errıe 13 
iE2 id ind iEmt 


E+N+ Ern+e® E+nt® Erm+ 


Wegen der inneren Struktur der Transformation (2) und der augenscheinlichen Symmetrie der 
Ausdrücke (13) haben wir jedoch nur 


er: 1 TER „[ tent 
Flargrel "ers "lesen: ler] 09 


zu berechnen. 


Zum Ausgangspunkt wählen wir den Wert des uneigentlichen Integrals 


e-i(Satny+Zz) Ind : v 
er ad en (14), 


bei dessen Berechnung wir uns (wie schon bemerkt) auf die Cauchyschen Hauptwerte beschrän- 
ken können. Damit fallen im folgenden manche mühsame Konvergenzbetrachtungen weg. 


Zuerst hat man unter Benutzung von 
= 0 für A<|y] 
[ Joan) cos ydn ae leg 
ö ya2 


und der BER SE Relation 


cosl 2 Ra ee _® (ASIAN) -IVere 
ae war ayaym 2 | y®+ 22 
für den gesuchten Wert von (14) 
j= 3 [esse Ku Tarn 
0 0 
—4n | cos&xde 4 PELIREI ARTEN EN, (An)cosnyd 
mz IVR+ tz z2 oA nyY N 
0 0 


re 4 les tete, 
J V 02 — 99) (A? + 2) 
® lvl 


Führt man nun anstatt A eine neue Integrationsvariable ® durch 


2?+2=(P? +2)cho 


ein, so findet man sofort 


oo 


Nee | et rrne 1, K,(ey® + 2), 


| Ya— = @+ Ei 
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und es ergibt sich endlich 


9 


—_—— 2m: 
een, TEE | de een 
” |Rsters + 2) cos &x.dE nt 
0 
Wir kommen so zum Ergebnis 
j2 EINEN qua 

ee ‚r=y# N 14.1), 

Nlerrrete® a Re ze 14.1) 


das sich mit Hilfe von einigen Soninschen Tatsachen aus der Theorie der Zylinderfunktionen 
sehr schön auch durch Verwendung von Polarkoordinaten gewinnen läßt. 


Nach (2) kann man auch schreiben 


Aral: EUER = =, 2 rap ran (15) 
& + n° 4 62 p) r E i ’ 
woraus sich unter Anwendung von (4) die Ausdrücke (13.1) leicht berechnen lassen. 
Für den ersten haben wir 


% 5 N 1 RE 
er re ee “ 


während die zwei nächstfolgenden durch Verwendung der Hilfsrelation 
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zu finden sind. Es ergibt sich ‚dafür 
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5. Die bisherigen Ergebnisse machen es nun ganz leicht, die fertige Lösung unseres Pro- 
blems aufzuschreiben. Nach (10), (12), (3) und (16) erhalten wir sofort 
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BIFER = — [ Pr eng 1312 per 
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II. Der Spezialfall einer Einzellast 

1. Der jetzt zu betrachtende einfache Spezialfall ist wahrscheinlich im wesentlichen alles, 
was aus dem nämlichen Fragenkomplex berechnet worden ist. Mitgemeint sind natürlich auch 
die daraus direkt ableitbaren Tatsachen von der Rolle, welche die Singularitäten in der mathe- 
matischen Theorie des elastischen Raumes spielen, worüber man sehr schöne Ausführungen in 
der Literatur finden kann [2]. 5 

2. Wir wenden unsere allgemeinen Formeln (17) auf den Fall an, wo die Belastung lediglich 
aus einer im Ursprung und zwar in der positiven Oz-Richtung wirkenden Kraft Q besteht. Man 
hat nur in (17) 


X, y,2) = Y(,y,2)=0, Zi, y,2) = 0) I 
einzusetzen und die durch | 
III I By) 8 @—) 8 y—P)d@— y) da dB dy = f(x, y, 2) 


ausgedrückte Eigenschaft der Diracschen ö-Funktion zu benutzen, um auf das bekannte Er- 
gebnis [2] 
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zu kommen. 
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Ein Iterationsverfahren zur näherungsweisen Lösung 
von Matrizendifferentialgleichungen 
Von P. Bajcsay und V. Lovass-Nagy 


Die Formel für die Lösung einer linearen M atrizen-Differentialgleichung ist für die numerische Rech- 
nung besonders bequem, wenn die kanonische Darstellung der Koeffizientenmatrix explizit bekannt ist. 
Diese Formel ist jedoch schwierig anzuwenden, wenn sich die Matrix nicht auf Diagonalform bringen 
läßt oder ihre kanonische Gestalt unbekannt ist. Aufgabe der vorliegenden Arbeit ist es, ein konvergentes 
Iterationsverfahren anzugeben, das bei einer beliebigen Koeffizientenmatrix anwendbar ist. Dabei wird 
die gegebene Matrix in eine Summe aus einer Matrix X mit explizit bekannter kanonischer Darstellung 
und einer Matrix ® von möglichst kleiner Norm zerlegt. Die auf diese Weise erhaltene Lösung erscheint 
in einer nach den Rechtseigenvektoren von X entwickelten Form. 


The formula for the solution of a system of linear differential equations is particularly swited to numerical 
computation if the canonical representation of Ihe matrix of coefficients is known explieitly. If this is not 
the case, or if the matrix cannot be transformed to diagonal form, difficulties arise in the attempt to apply 
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the formula. The purpose of the present paper is to construct a convergent method of iteration that may be 
used with any matrix of coefficients. To this end ihe given matrix is presented as the sum of two matrices, 
U and B, where U is a matrix the camonical representation of which is knoum explieitly while B is a square 
matrix with a norm as small as possible. T’'he solution obtained in this way is a series expansion in terms 
of the right-handed eigenvectors of X. 


La formule de solution d’une lineaire &quation differentielle de matrice est specialement convenable & la 
calculation num£rique si la representation canonique de la matrice des coefficients est connue de maniere 
explicite. Si pourtant la matrice des coefficients n’est pas & transformer a une forme diagonale ou que la represen- 
tation canonique de la matrice des coefficients est inconnue, l’application de la formule de solution connue est 
diffieile. Le but de l’article present est la construction d’un proced£ iteratif et convergent qui est utilisable avec 
chaque maitrice de coefficient. Le procede presente la matrice de coefficienis donnee comme somme d’une 
matrice X, dont la representation canonique est connue de maniere explicıte, et d’une matrice carree d’une 
norme la plus petite possible. De cette fagon on obtient une solution qui apparait sous une forme developpee 
apres les vecteurs propres d’ X de cöte droite. 


®opmyıa AA peImeHun IHHeÄHOTO MATPHYHOTO AußhepeHlMAasbHOTO YPaBHeHUA OCOOeHHO 
YAO6HA AA YHCJIEHHOH 06PA60TKH, eCJIH KAHOHMYECKOEIPENCTABJIEHHE MATPHLBI KOSPHHNHEHTOB 
M3BECTHO B ABHOM BH]IO. ÖONHAaKO IIPHMeHeHne 3Toü hOPMyalbl 3ATPYAHHTEJIBHO, eCJIM IIPUBe- 
MeHMe yIOMAHYTOÜ MaTpHlbI K KAHOHNHYeCKOMy BHAy He yAaeTcaı MAIM KAHOHNyYecKmü BUN 
MAaTpHIIbI HEeH3BECTHBIH. Samayueii HacTommei PA6oTEI ABJIAETCH TIOCTPoeHHEe CXONAMETOCH 
HTEPaLMOHHOTO MeTONa, IPHMEHHMOTO AIA AO00TO BuNa MaTpHuupI KoabpummeHTtoB. JIanHası 
MaTpnıa Ko9ßPHIHMEeHTOB IIPH 3TOM NPencTaBllfeTcA B BUNEe CyMMbI MaTpuupI Ü c n3BecTHoä 
ABHOÜ KaAHOHHYeCKoü POPMOH H MaTpHusl ®.c HMO BOSKMOCTH Malnof HOopMoü. Ilosıyuaemoe 
3Necb pemieHme NpencTaBilnerca B BHNe Pa3lo;KeHNUA OTHOCHTENBHO IIPABbIX COÖCTBEHHBIX 
BEKTOPOB MAaTpupl W. 


Bekanntlich ist die Lösung der Differentialgleichung mit der konstanten Koeffizienten- 
matrix 


d 
ee Qi: 


wenn sie die Anfangsbedingung (0) = 4, erfüllt 
a ae EEE. Burn EN 


ferner ist die Lösung der Differentialgleichung 


d £ 
TEN) er nn ee 


wenn sie dieselben Anfangsbedingungen erfüllt: 


t 
a a a OL et ee en 


(Siehe beispielsweise [1], [2]). Die Formeln (1.1), bzw. (2.1) ergeben aber für das praktische 
Rechnen nur dann entsprechende Resultate, wenn die Transformation der Koeffizientenmatrix 
in eine Diagonalmatrix bekannt ist. Zweck der vorliegenden Arbeit ist die Ausarbeitung eines 
solchen Iterationsverfahrens, mit dem die Lösung der Differentialgleichungen (1) und (2) dadurch 
ermöglicht wird, daß sie in einer für das praktische Rechnen zweckmäßigen Gestalt dargestellt 
wird und zwar für eine beliebige Koeffizientenmatrix fl. Denn im allgemeinen kann die Trans- . 
formation der Lösungen der Formeln (1.1), bzw. (2.1) — im Falle einer beliebigen Koeffizienten- 
matrix $} — nicht so durchgeführt werden, daß an Stelle von Exponentialfunktionen mit Matrix- 
Argumenten nur solche mit Skalar-Argumenten stehen. Das folgende Iterationsverfahren da- 
gegen stellt eine Reihe von solchen Formeln auf, mit denen die exakte Lösung in gewünschtem 
Maße angenähert werden kann, die zudem ausschließlich aus Funktionen mit Skalar-Argumenten 
aufgebaut sind, ohne jedoch die Kenntnis der Diagonaldarstellung der Koeffizientenmatrix S 


vorauszusetzen. 


1. Teil. Lösung von Gleichung (1) durch Iteration 


Offensichtlich kann die Koeffizientenmatrix $ stets in der Form {= U + © dargestellt 
werden und zwar so, daß die diagonale Transformation von X leicht durchführbar sei (z.B. 
je nachdem, welches technische Problem die zu lösende Matrizen-Differentialgleichung lieferte, 
kann X als Diagonal-, zyklische, gleichmäßige Kontinuantenmatrix, usw. gewählt werden), 
während 9 eine beliebige Matrix ist. Dadurch kann Gl. (1) umgewandelt werden: 
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Nun ist die folgende Reihe von Differentialgleichungen zu bilden: 
d 
dt y, = Ad > YO) = Yo» 
d 
a 9 Ay + Bdyı ; (0) =, 


+ BH 5 Hu) = 5 


d 
15 Im = Am + B Ym-15 Ym(O) = do 
Es kann bewiesen werden, daß lim 4, = e*' 4. Sowohl beim Konvergenzbeweis des Verfahrens, 
MIX 


als auch bei seiner Anwendung muß selbstverständlich darauf geachtet werden, daß X und ® 
im allgemeinen nicht vertauschbare Matrizen sind! RE, 

Von den in (4) zusammengefaßten Differentialgleichungen kann unter Berücksichtigung 
der Zusammenhänge (2) und (2.1) die m-te Lösung in folgender Form dargestellt werden: 


De do + FEN B pm ale) BE... ee ee 
Oder aber bei Einführung der Bezeichnung 
Met + HM, ld. ee 
ist : 
ESCHE DE. ce Sun = aa De Be 


Nehmen wir nun an, daß W n-ter Ordnung ist und in der Form 


n 
=. 
v=1 
darstellbar ist, worin A, die Eigenwerte von W, u, bzw. vf dagegen die rechts- bzw. linksseitigen 


Eigenvektoren von X bedeuten; v7 u, = ö,, [3]. Dann kann der gemäß Formel (6) bestimmte 
Wert M,, folgendermaßen ermittelt werden): 


N N t 
no Arb AL —Ar. 
Den — 2 ft un dh, + Lob, Bun) Fe" dr u, DE, + 
man! 1 fi) 


n N t T 

E En Be 

> = 2 (or, ® ur.) (0%, ® u) [e TEE TE (a DE +... + 
ka= „= 0 


Ar n (6.1). 
jr > = PR 2. (0, B U) (vr, B 7) u 4 (On B Ur.) x 
klkel Km 
t T u 
% ; eg Fi Fuel PA Sans ne em Akm-ı) ? (dr)r—. Ur, O8, 
Mit (6.1) ist dann das Endresultat von (5.1) 
4 er Er 
Um iY zZ e ut (df, Yo) +2 (vi B Ur) (0X, Yo) fe a) KT dr + 
1 ky= 0 
n n t T 
Rs iS 5 (of, BY u,.) (of, DR Ur) (of, Yo) Sr er A)? (dr)® +... + 
k=1k=1 0 0 (5 2) 
N n n di. 
“1 = 95 ee 3 (0x, B U,) (vr, B U) a8: (bR, -1 ® Ur) (of, Yo) x 
Kkelhml Ay 
t c 
x [et f RAT... Jen em)? Col Ur 
M 0 0 2 
Diese Darstellung des Resultats in (5.2) bietet offenbar zwei Vorteile: 
!) Im Sonderfall, daß X und ® vertauschbar sind, d.h. AB = BW also (siehe [4], Seite 96) u, v* 9 
= Bu, of, nimmt (6.1) folgende Form an: 
Mn LE Bfar + Bf Slim gm fe... f (dem 
m kı Ok, SAH SP ON Ne — 
Nn Ant * Bt 9:2 Ym-1 gm-1 Bt Pr Ym-1 gm-ı 
— Bo 1 b & _ ——uons — u 0 —— —— ... 
ie A an er pn m! Re A. a tr in 
wobei zu beachten ist, daß 


* * * 
un, (dk, BuR,) (0%, Bun) (08, Bur,) dt, = Blur, (dr, um) (dheum)-- (vr, _ı Ur,) vi, = Poly bp. 


RN 
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Erstens: da die in der Formel vorkommenden Produkte (br, B u,,) und (br, 4,) skalare 
Werte sind, kann auch in dem Falle, daß ® nicht in eine diagonale Form transformierbar ist, 
eine solche Lösungsformel gefunden werden, in der nur Funktionen mit skalaren Argumenten 
vorhanden sind. 

Zweitens: die Spaltenmatrix ),, kann in die Komponenten der rechtsseitigen Eigenvektoren 
u, der Koeffizientenmatrix A aufgelöst und so dargestellt werden; die Darstellung in dieser 
Form kann auch als Verallgemeinerung der in der theoretischen Elektrotechnik allgemein ge- 
bräuchlichen ‚,Methode der symmetrischen Komponenten“ aufgefaßt werden. 

Die durch Formel (5.2) gewonnene Spaltenmatrix )„ kann als die Näherung (m—1)-ter 
Ordnung der Lösung (1.1) der Matrizen-Differentialgleichung (1) betrachtet werden. Auf Grund 
des im III. Teil abgeleiteten Konvergenzbeweises ist festzustellen, daß die ersten m-Glieder 
der unendlichen Potenzreihe der in (6.1) dargestellten Matrix M,(d) eigentlich identisch sind 
mit dem Taylor-Polynom (m—1)-ten Grades der Matrizenfunktion e**. 

Zur Abschätzung des Fehlers zwischen der genauen Lösung der Differentialgleichung (1) 
und der in (5.2) angegebenen Näherung kann das folgende Verfahren abgeleitet werden: 

Es sei 

br I Bu,| sa 
und 
bot, ol SB 


ferner im abgeschlossenen IntervalO sr stsa 
er) t <s b 


und 
uie®|<c, 


dann ist?) : 
| R% n N n . 
A 
y— Ym| = RE 4 Ur, 124 zZ #; 2 (0%, BUr,) (08, B Ynrı) (DE, y) x 


ans MET 


EN LSU I eslen, Aveo REN m Be 


mi V: [2 m! yYm—(naba® 


m 


t T = 
x F ee 3 ee f em m) ? (dr)” En ER | 
0 0 


vorausgesetzt, daß m genügend groß ist, damit m>naba. 


II. Teil. Lösung von Gleichung (2) durch Iteration 


Die inhomogene Gleichung (2) kann ähnlich der Umwandlung der homogenen Differential- 
gleichung (1) in (3) auch so ausgedrückt werden: 


d 
rer rdy+g®; MIO) Sn RN RT ATI. (7), 


worin die Transformation von W in eine Diagonalmatrix bekannt ist, B dagegen eine beliebige 
Matrix bedeutet. Bei Beachtung des im I. Teil Erwähnten ist die Näherung (m—1)-ter Ord- 
nung der Lösung der Differentialgleichung (7): 


Ym = Mn (Yo + ! DEE er Een (7.1) 


Zur Abschätzung des Fehlers zwischen der genauen Lösung der Differentialgleichung (2) 
und der in (7.1) angeführten Näherung kann das folgende Verfahren durchgeführt werden: 

Es seien folgende Bezeichnungen gewählt: 

Die exakte Lösung der Differentialgleichung (2) ist gemäß Punkt 2.1): 9 =Yy,-+ ps 
worin d, die in (1.1) angeführte Lösung der homogenen Gleichung (1) und y, die Partikular- 
lösung der inhomogenen Gleichung (2) bedeuten; 

die näherungsweise Lösung der Differentialgleichung (2) ist gemäß (7.1): 9, = Ymn + Ump» 
worin Y„ die in Punkt (5.2) ausgedrückte näherungsweise Lösung der Differentialgleichung (1) 
und 4,» die durch Iteration gewonnene Partikularlösung von (2) darstellen; 


2) Ist = [x], dann ist |r] = V EN 
i=1 
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es sei 
ok, B u,,| sa 


Ivk,Yl SP» 
ferner im geschlossenen IntervallO <rstsa 
)! <b 


und 


und 


unse" <c, 
außerdem 
Ik, 80) = Y 
und 
u. eu] <q; 
dann ist 


19 — Ym| = 9a — Ymn + 9» — Ymp| S 19m — Ymn| + 195 — Umpl= Rn+ RB, . . (7.2), 
worin eine obere Grenze für R, laut Formel (5.3) folgendermaßen gewonnen wird: 
(n aba)" m 


= ee (m>nabo). 


R,< yücß 


h R, = |9p — Ymp| = 
== 1} e u, | > 2 2: (0x, B U) a (Of ® Upnrn) (Ber g9()) x 


rail km+ıml 


< mar fe fr] 
(naba)” 1 (n ab)” m 


<yndya — —=yYndya —————; (m>nabo). . - . (7.3). 
use) /ı 2; m) age N 


i—t T N 

(ax, — A)? (An 1 m)? 

x e A m+1 ... 
J J e m’ (de” + .--t dr 


Unter Berücksichtigung von (5.3) und (7.3) wird das Resultat für (7.2) schließlich: 


B. -(naboy" m 
ee in m! Ym2—(naba) 


(ß+dya; (m>nabo). 


II. Teil. Konvergenzbeweis des Verfahrens 


Aus den Fehlerabschätzungen (5.3) und (7.4) ist unmittelbar ersichtlich, daß der Absolut- 
wert der Abweichung zwischen den exakten Lösungen der Gleichungen (1) bzw. (2) und den 
Näherungslösungen (5.2) bzw. (7.1) gegen Null konvergiert, wenn die Zahl der Iterationen 
im >:c0.36f, 

Ersetzt man andererseits in der Formel (6) die Matrizenfunktionen e*! und e!@=) durch 
ihre unendliche Taylorreihe (wobei beachtet werden muß, daß im allgemeinen X und ® nicht 
vertauschbar sind)?), so gilt 


Mn) =E+AHBDEH, ABEL AHBPP He 4 


1 
k N Yym-ı ma N u Sr). Se 
| (m = 1) | (A t B) si Rn nr (8), 
worin Rt, die folgende unendliche Matrizen-Potenzreihe bedeutet: 
1. 
Rn = le, 
R 2a D 
hierin bedeutet die Matrix $’ nichts anderes als die Summe derjenigen Glieder in der Entwicklung 
der Potenz (A + 3)”, in denen ® höchstens (m—1)-mal als Koeffizient vorkommt (die Zahl 


m—l 
dieser Glieder ist: Bahn Offensichtlich stimmen die ersten m Glieder des in (8) auftretenden 
k=0 


Ausdruckes M,,(2) mit dem Taylor-Polynom (m—1)-ten Grades der Matrizenfunktion et! — e(!+%) 


2) Also ist beispielsweise (U + B?=-WHABLBAHB., 


ZAMM 39 (1959) Heft 1/2, Seite 13—19 183 


überein, in jedem von ? abhängigen Glied v-ten Grades (v > m) von PR, ist der Matrizenkoeffi- 
zient von 


vr 
v! 


identisch mit den ersten St ) Gliedern der Entwicklung von (X + B)”. Die Richtigkeit dieser 
2 m— 

E Feststellung leuchtet durch vollständige Induktion leicht ein. Überschreitet also die Zahl der ' 
Iterationen alle Grenzen (m— &), gilt tatsächlich 


Mt) > et! = eA+ BL, 


. was eingangs behauptet worden war. 
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- Vibrations of a beam on elastic foundation II” 
By P. M. Mathews **) 


In einer früheren Arbeit!) hat der Verfasser die Schwingungsgleichungen eines durch eine hin- und her- 
gehende Kraft belasteten Balkens auf federnder Grundlage unter Vernachlässigung der Dämpfung für den 
Fall gelöst, daß der Angriffspunkt der hin- und hergehenden Kraft sich längs des Balkens verschiebt. In der 
vorliegenden Arbeit wird nun dem Dämpfungseinfluß Rechnung getragen: explizite für alle Frequenzen und 
Geschwindigkeiten der aufgebrachten Kraft geltende Lösungen werden erreicht. In qualitativer Hinsicht unter- 
scheidet sich der vorliegende Fall vom ungedämpften vor allem dadurch, daß bei einer von Null abweichenden 
Geschwindigkeit die Verschiebungen des Balkens um die momentane Kraftangriffsstelle nicht mehr sym- 
metrisch verlaufen, sondern hinten größer sind als vorn. 


In an earlier paper !), the equations describing the vibrations of a.beam on elastic foundation under the 
action of an alternating force whose point of application moves along the beam was solved by the author, 
neglecting damping. The effect of damping is taken into account in the present work, and explicit solutions, 
valid for all frequencies and velocities of the applied force, are obtained. In qualitative features the most 
obvious difference from the undamped case is that, for non-zero velocity, the displacements of the beam are 
no longer symmetric about the instantaneous position of the applied force, but are greater behind than ahead 
of this position. 


Dans une etude anterieuret) l’auteur a r&solu, sans tenir compte de l’amortissement, les equations concernant 
les vibrations d’une poutre sur fondation Elastique sous Veffet d’une force «va et vient » dont le point d’applica- 
tion se deplace le long de la poutre. On tient compte, dans la presente Etude, de l’effect de V’amortissement et 
on y parvient ü des solutions explicites valables pour toutes les frequences et toutes les vitesses de la force mise 
en jeu. ‚Du point de vue qualitatif ce qui differencie le plus le cas amorti du cas non amorti est que pour une 
vitesse non nulle les deplacements de la poutre ne sont plus sym£triques par rapport 4 la position instantanee 
de la force mise en jew mais sont plus importants derriere cette position que devant elle. 


B npensinyımem cooÖleHun!) ypaBHeust, OINChIBAlWIMe koleÖannn ÖallıKu, Ha YoOpyToM 
OCHOBAHHH, IION MeÜCTBueM TIepeMeHHOÄ CHJIBI, TOYKA IPHIIOSKEHHA KOTOPOH NBIHKETCAH IIO 
Öarıke, PeIleHbI ABTOPOM He YYuTbIBası nemıhupoBaHnnn. B Hacronmei pa6oTe nemnmbnpoBanne 
yYYuUTbIBaeTcA M MOJIy4eHbI ABHBbIe PellieHMA CHPABeNAMBbIe IA BCcexX YACTOT U CKOPOocCTeH. 
UTo Kacaetca KayecTBeHHoÄÜ xapakTepuCTUKu, TO HaiboNee ABHOe OTIAIHYHE OT CIYyar NeM- 
IhHPOBaHHA 3akIloyaeTcHB TOM, YTO NIA CKOPOCTEH OTIMYAIINUXCH OT HYIIA CMeIIeHuA Oalıku 
y’Ke He CHMMETPHYHEI BOKPYT MIHOBEHHOTO HOJIOFREHUA IIPNHIIO}KeHHOÄ CHJIbI, A OKA3bIBAIOTCH 
601bIIMMH 3A 3TOM TOYRKoü YeM IIepen Heil. 


I. The Equation of Motion 
It was shown in Paper I that if 


— Mass per unit length of the beam, 

= Young’s modulus of the material of the beam, 

Moment of inertia of the section of the beam, about the neutral axis, 

Rlastieity of the foundation 

damping force per unit length of beam per unit velocity, and if the beam is taken to 


I 


S muy 
I 


and 


I 


*) Contribution from the Division of Pure Physies, National Research Couneil, Ottawa, Canada. 
Issued as N.R.C. No 4931. 
er) National Research Laboratories Postdoctorate Fellow. 
!t) P.M. Mathews, Z. angew. Math. Mech. 38 (1958), S. 105—115 hereinafter referred to as Paper T. 
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extend along the z-axis, then the function y(z, f) giving its lateral displacement at the point x 
at time Z satisfies the equation 


oy 0°y 0y 
a han. I == u ve 1), 
er a a (1) 


the term on the right hand side representing the externally applied force whose point of appli- 
. cation moves in the positive x direction with velocity v, and whose magnitude varies as F, cos ot. 
It was further shown that the steady state solution of (1) can be written as 


yr,d=yl)eosat+hdsnotl „u. vu va 0) 


EI 


where 
Tran ee ee > Se 


and y,(r) and y,(r) have Fourier transforms p,(s) and p;(s) given by 
ö re. 
a FE It —o®s®2—yvis+k— god? +lRvweis— ro? 


and 


P:(s) = Fy (E It —o® 2 —vyvis+k— oo" + (2vwois— vo) 


— (2vwpis—vo) Be 


The analysis revealed the existence of a resonance frequency ®, for the beam, and a critical 
velocity v, representing the limit beyond which the vibrations would be unstable in the un- 
damped case: 


a 
A & ra 2 EN as ap ENG 
0 2 (6) 
and 
4E Ik\\V4 ” 
Di ur serie Elan ee ER) 


Now it is easily seen, e. g. by considering a mode of free vibrations in which the motion 
isin phase all along the beam, that the damping constant » also has a ‘critical’ value »,, given by 


u EU VRR ET NEE SEE . (8). 
It is convenient to use &, d, and », and the constants 
kr\le 4EI\# F, 
d = 77 and = r SE Teer) 


which were defined in Paper I, to express all the quantities of the problem in non-dimensional 
form. Thus, on defining 


N nr DE Me ER 
N n En S at R=ra, and'wy Fun (10) 
and 
- 00 
Dis —iRS 2k 
‚(S) = YuR)emiEldRn F P;($) 1 il, Daran Dee = 111), 
Ars 0 


equation (2) becomes 
YRd) = YlA)eosat # YılR)sunwt . ern (12), 


Yı(R) and Y;(R) now being Fourier inverses (with respect to S) of P,(S) and P,;(S) where 
8 [—4V?®—8N VIS+A4(— W3] 


Pe - Au. - 
ı( ) [S* A V? Ss? 8 v7 s 4 (1 ee, w»)]? + 64 [V Wi S N w] . (13) 
and 
PS) — —_64(VWIiS— NW) 
[?—4 ?®—8N ViS+4A— WöPR+ 6a[Vwis—Nn we‘: (19. 


‚ Inversion of (13) and (14) can be carried out if they can be split into partial fractions, for 
which it is necessary to factorize the denominator, or equivalently, to solve the equation 


[—4 V?®—8NVIS+A4A— WYP +6 IVWIS-NWE=O... (ib). 


A sun Zur 
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Inspection of (15) shows that its roots must occur in pairs like 
A,+iB, As—1B, —AtiB, —Au-iB (16) 
—A,+iB,, —A,—iB,, A;+iB,, AyelB,., 25 


where A, ..., Ag By..., B, are all positive. This information about the signs of the real 
and imaginary parts of the roots is necessary before the formulae of Paper I (equations (24) 
through (27)) can be used for inversion. It is derived from the fact that the equations (15) cannot 
have any real or pure imaginary roots when N, V, W are all non-zero and therefore the signs 
of the real and imaginary parts of any root must remain constant over the whole range of values 
of N, Vand W. Thus we deduce the pattern of signs in (16) from the known signs in special 
cases, e.g. when V<1,NZ0,W=Z0as ae by analytical continuation from the case 
V=-0,NZ0,WZ, discussed below. 


We now have eight parameters A,,..., Ay B»-.., B, as against three in the case 


| 
en 


N =, and therefore the solutions are extremely unwieldy. It would appear to be simpler in 


any particular case, when the values of N, V and W are given, to find out the roots (16) and then 
proceed to carry out the inversion. We shall, therefore, confine our attention here to the simpler 
cases: 


WILSUANEU MR, 
MER EN EU FENZU, 
and compare the solutions with the corresponding ones in the case of no damping, discussed in 
Paper I. 
II. Vibrations in the damped case 
(i) Alternating load at the origin: V=0, NZ0,W=0. (13) and (14) become 


8 [St + 4(1— W3] 


A) [st +4d— Wp + 64 N2 w2 Fer a (17) 
and 
64 NW 
PS) [rad wWopRteamw 00° a Re) 
and equation (15) reduces to 
[St +41 — WW)? +64N?W?=0, 
its solutions being 
Sem At BB -+LAF Si... lan: (19), 
where 
A=Acos®, BE= MSN DABEI EIER (20), 
with 
Amill10ll = Wer aba N We u a ee (21) 
and 
—2 NW 
= UBS N) rar 2a, 
tan4 ® d_ ws’ (a <Aad®<O) (22) 


On splitting up P,(S) and P,;(S) into partial fractions and carrying out the Fourier inver- 
sion using the well known procedure we obtain, after some simple algebra, 


Y,(R) [e-#IRI[B (3 A? — B2) cos A|R| + A (3 B?— A2) sin A 


2 
= m Ri) 
(A®+ B3) 
+ eAR| {A (3 B?— A?) cos B|R| + B (3 A? — B®) sin B|R|}]. ee N 
and 


2 r 
YıR) = er e-BIR{A (3 B®— A2) cos A|R| — B (3 A? — B2) sin A|R]|} 


+ eAlR| {B (3 A? — B?) cos B|R| 


AßB®—AY)sinB|R|}] .. - .. (24). 


It can be easily seen that when N=0, ® = T so that we have A = B from (20), and 


hence Y, = 0 and Y, reduces to the expression derived in Paper I (equation (33)) if the latter 
also is expressedin nondimensional form. WithA und ©chosen as parametersinstead of N and W the 
variation of Y,und Y, with R ispresented in Figs. 1and2 wherei? Y, and A? Y,are plotted against 
AR for various values of ®. For any given N and W, the le values of A and D can 
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be found from (21) and (22) and then Y, and Y, determined from the appropriate curves in Figs. 1 
N y 
and 2. For instance, in the case of resonance (W = 1) we have ® = 3 while A depends on N; 


if W or N vanishes, ® = — and when W>1,9—0. 


30 


XY(R) 


) 2 3 4 5 7 8 9 [e) 
AR 


NT 
“ - a 
: te 
= 
10 
(0) | 2 3 4 5 6 Ti 8 9 IO 
AR 


Fig.2. PıR)inthe case Y=0,W20,N #0 


It may be noted that while Y, is small compared to Y, for small N or W (9- a ss 


magnitude increases to a value equal to that of Y, if the applied frequency is high (W>1, 
® — 0). However, the absolute magnitudes of both Y, und Y, decrease rapidly as the frequency 
is increased above the resonance frequency. In fact it is easy to see from (23) and (24) that if 
Ww:—-1>2NW, Y,(0) & — Y;(0) 2 — (2 WS) -U2, 

It is also interesting to note the dependence of Y,(0) and Y,(0) at resonance on the damp- 
ing N: 


YO) Ir-ı =2{B1)(cos + sin) m 2... (8), 
[YxO)lıraı = 2(y2 + 1)(cos sin EBEN... ...... 00). 


(ü) Amoving load of constant magnitude: W=0,N 0, V+0. 
In this case (13) and (14) reduce to 
8 


P.(S) = dr 
(5) HAWRIENVISIE terre. 
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and 
ENTE TE Er RETTET ET . (28). 


As usual, as a preliminary to the partial fraction expansion of (27) before inversion, we 
have to solve the equation obtained by setting the denominator equal to zero: 


AV EN VESU-AIS0T ee (29). 
Now, 5’ = iS obviously satisfies the equation 
SA +4 S2?—-8NVSHA-=0, 


whose coefficients are all real, and by the application of Sturm’s theorem?) to this equation it 
is easily ‚seen that its roots are all complex or two real and two complex — and consequently 
those of (29) are all complex or two pure imaginary and two complex — according as 


AT —N?) Ver (8 36 NE I. N) Vaed > on... m. or 


E a ae 


0 | 2 S a 5 
N 


Tig.3. The curve 4 (1—N?) V°— (8 — 36 N® + 27 N‘) V! + 4 = 0. The nature of the vibrations depends on whether V, N values are 
such as to be to the left or to the right of this curve, the appropriate solutions being (35) and (41) respectively in the two cases. 


When the roots of (29) are all complex, they must be of the form 


= AT FTIBE N EBAFI Ber a Tel BRD 
where A,, A,, B are positive quantities satisfying the relations 
AHA: 2 BU=AVE ne et 
BIEANZ AN a ee aa 
ER EEE RE EEE RE Re): 


The values of A,, A, and B for various values of N and V are given in Table I. 
In this case inversion of (27) yields the result 


eBRH(R) 4A, BcosA, R+ (AB®— A: + AN) sin A, R 


Sy erregen 16 A? BE + (ABA? + Ay) 
| ORH-R) 4A, Boos A, R— (AB + Ar Ad sin Ay R > 
A, 6%: BLAB+A— Ay ea 


where H represents the Heaviside unit function, i.e. H(x) = 1for x > 0; H(x) = 0 otherwise. 
The most striking feature of the above expression for Y,(R) is the disappearance of the symmetry 
between the ‘forward’ and ‘backward’ directions (R>0 and R<O respectively) which had 
been present in all the cases discussed so far. The first consequence of this is that the position 
of maximum displacement of the beam is no longer the point of application of the load, but a 
point behind it. Furthermore, the instantaneous shape of the displaced beam which has the 
form of a damped wave, appears as if compressed ahead of the load and stretched out behind. 
This tendency gets accentuated rapidly as either N or V is increased, as is evident from Fig. 4. 
In fact, as will be shown below, if N and V are large enough so that the second of the two alter- 
natives in (30) holds, the shape ahead of the moving load is still that of a damped oscillatory wave, 
while the displacement behind has the form ofa difference oftwo exponentials and after attain- 
ing a maximum, decreases uniformly with increasing distance. 

Let us therefore suppose that N and V are such that two of the roots of (29) are pure 
imaginary and two complex. Then these must be of the form 


EB un 000.0, 86), 


2) See, for example, J. V. Uspensky, Theory of Equations (McGraw-Hill, 1948), Chapter VII. 


5) 


4 
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Table I 
BR 0.5 | 1.0 | | 2.0 | 3.0 Br o |» 
12 (0341) | (0215) | (0.158) (0.126) 
7 1.620 | 1.560 1.542 .534 
B=B=B, | 1.089 | 1.060 1.067 1.070 
SER A 1.874 | 2.106 2.294 2.455 
0 A a ne 1.000 | 1172 | 127 1.358 
2 B, 1.572 | 2.000 2.287 2.514 
= ua RR 5 109 ER 0.5998 | 034 | 0254 | 0.202 
Wi . m. 012625 | 2905 3.108 3.295 
“a TE ek 1.063 | 1.258 | 1.414 1.545 
1.0 ds er Fi B, 1.858 | 2346 | 2701 | 2.989 
B Br B, 0268 | 012 017 | 0.10 
< r ad e EEE EN Penn SER \ . 8 | =» 
A | 425 | 446 | 4639 | 4.818 
a a te ee 0.863 | 1.149 1377 | 1.568 
a0 Fi a ra her 1591 | 2211 | 26900 | 3.085 
B, 0.134 | 0.086 | 0.064 | 0.081 
Re 2 6.093 | 6.199 | 6322 | 6.452 
AS ee | 0641 | 099 | 1164 | 1.381 
an = te ler Se | 1193 | 1781 | 226 | 2.728 
B, | 0.089 | 0.057 0.02 0.034 
I 8.022 | 8.05 | 818 | 8.244 
Zr 2 a Se en | 049 | 0727 0.949 1.156 
z I 090 | 1412 1.865 2.286 
B, 0.067 |, 0.03 | 0.032 | 0.02 
A 10.022 | 10.050 | 10.089 | 10.146 
en RR RE 0.398 | 0.592 0.782 0.965 
| 2 RE OS Met, 0.722 | 1150 | 1.538 1.910 
B, 0.054 004 | 005 | 0.080 


For any N > 1, the veloeity at which the vibrations change character from that represented by equation 
(35) to that given by equation (41) — or vice versa — is denoted by V. in the table above, and its value is 
given in brackets. Y. of course is given by the ordinate of the curve of Fig. 3 corresponding to the specified 
value of N. 


Fig. 4. Y,(R) in the case W = 0 as given by (35), for various values of F and N. Unlike the eurves for the undamped case, the present 
ones are not symmetric about R=0. 'Ihe rapid increase of asymmetry with Y and N is evident 


where A, B,, B,, B, are positive and satisfy the relations 


2 Bm Bun Bade Di ana sen ad ne 
#+B—- B-B—BB av 
B(A+B—B BAND re 
(A+B)BB,=4 un weg ee 


Their values are also given in Table I for various values of N and V. 
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The solution for Y,(R) is now | 
eBEH(R) A(2B,+B,+ B,) cos AR—{[A—(B, +B,)(B, + B;)}sin AR 


Yı(R)=8 
R) A {A? + (Bı + B2)?} {A® + (B, + B,)} 
eB:R H(— R) ı j eb» RH(— R) 41 
AA B+ By} (BB) {A®+(Bı + Bo} (Be—B,) I 


_ The displacements, for R > 0 as given by the first term in square brackets, and for R< 0 
' as given by the second and third terms, have of course the forms described in the preceding 
paragraph. The curves in Fig. 5 illustrate the variation of Y,(R) with R in the present case. 


Fig.5. Y,(R) in the case W=0, as given by (41). 
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1 : 
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Über Schraubentriebe mit parallelen Achsen zur 


kombinierten Grob- und Feineinstellung 
Von W. Noli 


Aus Schraubenmutter und Schraubenspindel mit größerem Spiel oder aus zwei Schraubenspindeln lassen 
sich Schraubentriebe herstellen, mit denen man nach Wunsch groben oder feinen Vortrieb erhalten kann. Die 
Berührungsbedingungen dafür benötigter Schraubflächenpaare mit Linienberührung und Wege zu ihrer 
Formgebung werden aufgezeigt. 

Serew and nut with sufficient play, or two male screws, can be used to construct a gear that gives optionally 
coarse or fine advance. For the two screw surfaces required for this purpose the conditions of line contact as 
well as modes of forming are indicated. 


D’un ecrou et d’une tige filetee avec jeu assez libre ou de deux tiges filetees on peut construire des vissages 
au moyen desquels on peut obtenir sur desir un avancement approximatif ou minutieux. Les conditions de 
contact des paires de surface a visser avec contact des lignes y necessaires et des chemins a leur faconnage sont 
indiques. 

3 rafiku U BHHTOBOTO INIINHNEIA C 6ONIbIINM 3a30POM HM M3 ABYX IIIIMHNeJIea MOHKHO 
COCTABUTb YePpBAYHYI Tepenayy, HPH HOMOIM KOTOPOH IIO FKeJIAHMIO MOFKHO MOAIYYHTb CROPy 
HJIH MenIeHHyIo Nomayuy. NVXKA3bIBAITCH YCAHOBUHA COMPHKOCHOBEHUA YIOTPEÖNAEMBIX IIPM 3 
TOM Hap BHHTOBBIX NHOBePXHOCTef C COIPHKOCHOBEHHEM BAOJIBb KPUBOH N MAlTcA yKaaaHun 
Ha TO, KaROH BUN TPpeöyeTcA UPHAATb 9TUM HOBEPXHOCTAM. 


I. 


Werden, einem aus der Praxis kommenden Vorschlag gemäß, eine Schraubenspindel und 
eine Schraubenmutter mit gleicher Ganghöhe h und größerem Spiel r, — r, wie in Bild 1 ange- 
ordnet, so entsteht ein Schraubentrieb, mit welchem sich je nach Wunsch eine Grob- oder eine 
Feineinstellung erzielen läßt. Hält man nämlich die Mutter fest, dann führt eine lockere Drehung 
der Spindel um den Winkel x, zu einem Vortrieb vom Betrage 


en. en ad 5. m ld) 
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wobei p = h]2 rn den Schraubparameter (die reduzierte Ganghöhe) bezeichnet. Drückt man hin- 
gegen die Spindel gegen die Mutter, so daß ein Abwälzen der einander längs einer Erzeugenden e 
berührenden Zylinderflächen erfolgt und auch die Mutter eine Drehung um den Winkel 0%, -Fı/Tz 
vollführt (was durch Längsriffelung der Zylinder erzwungen werden kann), so beträgt der Vor- 
trieb der Spindel nur mehr i 


k=n—n=Pp u—o)=- nz BE (2), 
2 

Durch Festhalten oder Los'assen der Mutter läßt sich daher wahlweise ein größerer oder kleinerer 
Vortrieb erreichen. 

Für die Formgebung der Schraubflanken ist die Forderung maßgebend, bei der zweiten 
Bewegungsart — die die genauere Passung verlangt — wenigstens Linienberührung zu haben. 
Es liegt mithin folgende geometrische Aufgabenstellung vor: Zwei Schraubflächen ®, und ®, 
mit parallelen Achsen a, und a, und derselben Ganghöhe h sind zu konstruieren, die einander 
längs einer Linie berühren. 

Zur Untersuchung der bei solchen Schraubflächenpaaren obwaltenden Verhältnisse stellen 
wir uns zur Vereinfachung der Ausdrucksweise die Schraubachsen lotrecht vor. Da eine Schraub- 
linie vom Radius r die Steigung tgß = h/2rn = p/r auf- 
weist, können wir die zu einem Punkt P gehörige Schraub- 


GG Do a tangente jeweils durch einen Vektor PT kennzeichnen, 
ern dessen Vertikalkomponente den Betrag p hat, während 
IR = seine Horizontalkomponente durch eine Viertelschwenkung 
a aus dem Achsenabstand Pa = r hervorgeht. (Ich habe 


Im ” solche Vektoren in der früheren Arbeit [1] ‚„Tangenten- 
> 8 
re zeiger‘‘ genannt.) 

a Es sei nun P ein Punkt der Berührungslinie der 

beiden Schraubflächen ®, und ®,.. Die gemeinsame 


Bild 1 


Tangentialebene in P enthält die beiden Schraubtangentenvektoren PT, und PT. daher auch 
den Vektor T,T,. Da dieser, wie aus dem Grundrißbild 2 zu entnehmen ist, zur Verbin- 
dungsebene der Schraubachsen a,, a, senkrecht steht, so gilt der grundlegende 

Satz 1: Berühren einander zwei Schraubflächen mit parallelen Achsen und 
gleicher Ganghöhe längs einer Linie, so läßt sich ihnen längs dieser Linie ein 
zur Verbindungsebene des Achsenpaares normaler Tangentialzylinder anschrei- 
ben. — Die Berührungslinie stellt daher einen gemeinsamen Umriß der Schraubflächen für 
Normalprojektion auf die Achsenebene (Aufriß) dar. 

Zur Ermittlung eines berührenden Schraubflächenpaares ®,, ®, sind verschiedene Wege 
gangbar, von denen drei nachstehend skizziert werden sollen. 

1. Man nimmt etwa die Fläche 2, beliebig an, d. h. man verwendet für die Spindel irgendein 
gebräuchliches Gewinde. Die Flankenfläche ®, der Mutter ergibt sich dann als Einhüllende der 
um a, verschraubten Fläche ®,. Die Berührungslinie k läßt sich dabei nach bekannten Regeln 
der darstellenden Geometrie (vgl. etwa [2]) als Umriß von ®, bei Normalprojektion auf die (Auf- 
riß-)Ebene [a,a,] konstruieren, ®, kann anschließend durch Verschraubung von k um a, erzeugt 
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werden. Man findet insbesondere den Achsenschnitt von Ö, durch Verschraubung der einzelnen 
Punkte von kin die Ebene [a, a;]. . 


Wird zur rechnerischen Behandlung ein Normalkoordinatensystem x, y, zso gewählt, daß a, 


in die z-Achse fällt und a, durch x = —.d, y = 0 festgelegt wird, und ist ®, durch den Meridian- 
‚schnitt j 
En HENTZFAZRINTI NS) 
vorgegeben, dann lautet eine Parameterdarstellung von ®, 
Denen ya sn a ON SET 


Die Tangentialebene eines Punktes von ®, wird von den beiden partiellen Ableitungsvektoren. 


- (cos 9,, sin 9, Zı(r})) und (— Tr, - sin Q, Tı - COS Q,, pP) aufgespannt; sie ist, wie es die Berührungs- 


bedingung (Satz 1) verlangt, parallel zur y-Richtung, wenn 
As DC OEL RE RE ER 


‚Die Berührungslinie k wird daher mit r, als Parameter durch (4) beschrieben, wenn 9, gemäß 


ee FED STE 
eingesetzt wird. Die Meridiankurve von 8, kann dann durch Rückverschraubung der Punkte 
von k in die Ebene y = 0 gefunden werden; Schraubwinkel 9, und Schraubradius r, ergeben 
sich dabei unter Verwendung der eben bestimmten Werte aus (4) nach Bild 2 vermöge 
BR De: 005.050... Toreäll Oe.,. anzae ) 
durch 


gn= 7... RB=(d+P 4 P=-r+2d-r.cosp + d ld): 


Der gesuchte Meridian wird somit dargestellt durch 
u =n—d, y=0; 3=2—-7  =a+pr a —9) -- - - Ö). 
Dieses Verfahren ist beispielsweise angezeigt, wenn die Spindel ein normales Spitz- oder Trapezgewinde 
trägt; 2, ist dann eine lineare Funktion von r, (2] = const). 
2. Eine zweite Möglichkeit besteht in der Vorgabe des Grundrisses der Berührungslinie k, 
etwa durch seine Polargleichung 


T] — r (91) a era DE rer (9). 
Die Berührungsbedingung (5) liefert dann durch Integration 
| dr. 
4-—p [ten Er TREE A) 


wodurch der Meridian z, = z,(r,) von ®, festgelegt ist. Die Bestimmung von , vollzieht sich 
dann wie unter 1. 

In der vorangegangenen Arbeit [1] wurde beispielsweise die Annahme eines zur x-Achse symmetrischen 
kreisförmigen Grundrisses behandelt. Durch geometrische Betrachtung wurde dort gezeigt, daß die zuge- 
hörige Berührungslinie k eine Schraublinie mit derselben Ganghöhe h=2rp ist, was zur Folge hat, daß die 
Flanken &, und ©, zyklische Schraubflächen sind. 

3. Als dritte Möglichkeit käme schließlich auch die Annahme des y-parallelen Tangential- 
zylinders in Betracht, dessen Profil (Aufriß der Berührungslinie k) durch 


RER) TADEL DR ER AEUBAIRLE Ur ERTL EU CLT) 


gegeben sei. Die Schraubflächen 9, bzw. ®, entstehen als Einhüllende bei Verschraubung. des 
Zylinders um a, bzw. a,. Ein Punkt (x, y, z) der Berührungslinie k ist dadurch gekennzeichnet, 
daß sein Schraubtangentenvektor (—y, x, p) in der Tangentialebene des Zylinders liegt. Das 
führt auf die Bedingung 

‚ P P 

(X) = — — EN Te ee RE EEE US RRESER ER 6 71° 

ey 2) 
Aus (4) lassen sich dann die Größen r,, g, und z, berechnen, womit die Fläche ®, — insbesondere 
durch ihren Meridian z, = 2,(r,) — vollkommen bestimmt ist. Analog findet man ©,, wobei die 
Gleichungen (7) und (8) verwendet werden können. 
Wird der Berührungszylinder etwa als Ebene vorgeschrieben (2 = m x), dann sind die Flankenflächen 

®, und ©, Schraubtorsen (abwickelbare Schraubflächen), die einander längs einer geraden Erzeugenden 


k (y=—-p/m) berühren. Solche Schraubflanken sind beispielsweise von den Schrägzahnrädern mit Evol- 
ventenverzahnung her geläufig; sie erscheinen wegen ihres gleichartigen Charakters besonders empfehlenswert. 
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| 37: 

In I. bestand ein Schraubentrieb aus Schraubenmutter und Schraubenspindel. Eine kom- 
binierte Grob- und Feineinstellung läßt sich aber auch durch Verwendung zweier Schrauben- 
spindeln erreichen. 

Wir denken uns, wie dies Bild 3 zeigt, zwei Schraubenspindeln mit parallelen Achsen a,, a; 
von gleicher Ganghöhe Ah, jedoch entgegengesetzter Windung im Eingriff. Wie in I. ergibt 


1 =P:’%-»- ae Wer ARE nal air) 


den groben Vortrieb, wenn man bei lockerer Drehung der ersten Spindel die zweite Spindel fest- 
hält, die die Rolle der Mutter übernimmt, während man den kleinen Vortrieb 
Aa a El 2 ER REINT T 
Tz 
erhält, wenn man beim Drehen die beiden Spindeln so gegeneinander drückt, daß die längs einer 
Erzeugenden e berührenden Zylinderflächen aufeinander abwälzen. 

Eine solche Anordnung hat den Vorteil, daß die Spindeln — z. B. durch Ausschwenken der 
einen — leichter voneinander gelöst werden können wie Spindel und Mutter und dadurch ein 
handlicher Ein- und Ausbau erfolgen kann. Das ist besonders bei Maschinen vorteilhaft, bei 
denen öfter eine Änderung des Vortriebes benötigt wird. Man kann dann zu einer Spindel einen 
Satz dazu passender Spindeln anfertigen und nach 
Bedarf auswechseln. Auch der einseitige Lager- 
druck, der besonders bei der zweiten Bewegungsart 
entsteht, läßt sich für eine Spindel aufheben, wenn 
man zwei oder mehrere gleichartige Spindeln sym- 
metrisch um die erste anordnet; der Druck bei die- 
sen Spindeln verteilt sich dann auf mehrere Lager. 


u 


Bild 3 Bild 4 


Hält man schließlich die eine Spindel fest, so daß sich die andere Spindel abwälzend um 
sie herumbewegt, so führt ein Radienunterschied der Zylinderflächen ebenfalls zu dem Vortrieb (2). 
Eine solche Bewegung — es ist im wesentlichen die zweite Bewegungsart — ist eine Schrotung, 
wie sie H. Horninger in [3] eingehend untersucht hat. So kann der vorliegende Aufsatz auch 
als kleine Ergänzung in dem Sinne angesehen werden, daß auf Möglichkeiten einer praktischen 
Durchführung einer solchen Schrotung hingewiesen wird. Auch eine Arbeit von F. Hohenberg 
ist in diesem Zusammenhang zu nennen [4], in der allgemein die Zusammensetzung zweier 
Schraubungen untersucht wird. 

Zur Bestimmung der Schraubflanken wird man wieder von der Forderung ausgehen, daß 
bei der zweiten Bewegungsart Linienberührung vorliegen soll. Die Aufgabe lautet: Es sind zwei 
Schraubflächen 2, und D, zu bestimmen mit parallelen Achsen a, und a, von derselben Ganghöhe, 
aber entgegengesetzter Windung. 

Unter Verwendung der in I. erklärten Tangentenvektoren folgt aus dem Grundrißbild 4, 
daß für jeden Punkt P 

Aam Pl= ATıTEP., 


nn ., 4,2938 
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Nun liegt T, um p höher, T, um p tiefer als P. Eine gemeinsame Tangentialebene wird auf- 
gespannt von PT, und PT,; sie enthält auch den Vektor TOT, der Steigung = zu II, mit 
d= a,as. Da auch für jeden Punkt P T/T3 = d, folgt der 


Satz 2: Berühren einander zwei Schraubflächen mit parallelen Achsen, glei- 
cher Ganghöhe, aber verschiedener Windung, so läßt sich ihnen längs dieser 
Linie ein Tangentialzylinder anschreiben; die Erzeugenden dieses Zylinders 


haben zur Verbindungsebene der Achsen die Neigung tg» = 2 ‚ihre Normal- 
projektion auf diese Ebene ist parallel zu den Achsen. ap 

Zur Ermittlung des berührenden Schraubflankenpaares ®,, ®, sollen die in I. skizzierten 
Möglichkeiten betrachtet werden. 

1. Die Fläche 9, sei beliebig vorgegeben. Dann ist ©, als Einhüllende der um a, verschraub- 
ten Fläche 8, bestimmt. Die allgemeinen Bemerkungen in 1.1. gelten hier entsprechend. . 

Zur rechnerischen Behandlung sei ®, wieder durch den Meridianschnitt 


Be U ON el) ee EEE EN FE (3) 
festgelegt und die Parameterdarstellung 
N er (4) 


für ®, verwendet. Bei Berührung muß die von den partiellen Ableitungsvektoren aufgespannte 
Tangentialebene den Vektor T,T, = (0, d, 2p) enthalten. Es ist also 


cos 9, sing, zur) 
7 emo r con DD EDEN TEE DIANEE ERREILS) 
| 0 d 2p 
oder 
2Ppn de. ns) in 9 —pd.cosee 0. na ne (13a) 


die Berührungsbedingung. Dadurch ist die Berührungslinie k festgelegt. Durch Rückschrauben 
der Punkte von k um a, nach y = 0 erhält man die Meridiankurve von ®,-9, und r, ergeben 
sich aus (4) nach Bild 4 vermöge 

d—xt=T,:0050,, sTänor iR AeanR (14) 
durch 


un=-— ne dem per— 2197 008 rd) 


der gesuchte Meridian ist dargestellt durch 
neo rn, WU, 7,7. Dead tDio.—G), urn. (16). 


Für ein Trapez- oder Spitzgewinde ist 2/(rı) = m konstant, also 


LU re we a (13b). 


rn =; ä ; 
2p — m-d-singı 


2. Gibt man den Grundriß k’ der Berührungslinie k vor, so muß ihre Fortschreitungsrichtung 
(dx, dy, dz) in der Ebene durch den Schraubtangentenvektor (—y, x, p) und die Erzeugenden- 
richtung (0, d, 2p) des Tangentialzylinders liegen, also 


daes..duzı da 
—Y..2..P SD UNE DE RER. VRR AIR DER (17) 
Oundı 2D 
oder 
Bez de page 0 ENT, (17a). 


Diese Berührungsbedingung führt bei gegebenem Grundriß von k (etwa durch x = x{t), y = y(d)) 
zu 


ad 1 EEE 
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Ist z. B. der Grundriß von k eine Gerade «= ay + b, wird 


= Pla +a)y + $ 25a): Igy| + Const. ER ER > (18a) 

eine affine Logarithmuskurve. Für den Sonderfall y = c erhält man 
lm Eareltomt 22er eenne: . (18b) 

Be et 2.5 FT ons 


eine Parabel. che 2 ED 
Da die Spurparallelen s’der Tangentialebenen — wie aus Bild 4 ersichtlich — für alle P durch die Mittel- 


punkte $ der 7,T, gehen, schneiden die Normalen durch P die gleiche Gerade a,, für die aja, = a,a,. Daher 
sind auch Kreise um a, mögliche Berührungslinien. 
3. Schreibt man den Tangentialzylinder der Erzeugendenrichtung (0, d, 2p) vor, so erhält 
man durch seine Verschraubung um a, bzw. a, als Einhüllende die Flankenflächen ®, bzw. ®.. 
In einem Punkt der Berührungslinie k muß die Richtung der Schraubtangente (—y, z, P) 
in der Tangentialebene des Zylinders liegen. Ist dieser durch seine Spur 


ze) De NR 00) 
in y = 0 bestimmt, können wir seine Tangentialebene durch die Erzeugendenrichtung (0, d, 2p) 


und die entsprechende Tangente an (19), etwa (1, 0, z’(x)) festlegen. Für die Berührungslinie 
muß 


a a Re 
di:29 I 7 0 ERS BEREEREE 

bzw. F 
er er NE 
En 27@) (20a) 


Ist z(x) eine Gerade z= mx + 5, wird der Berührungszylinder zur Ebene z=emxz+b+2pyjd. 
®, und 9, sind Schraubtorsen, die sich längs der Geraden vom Grundriß 


_2p Eu 20b 
a a u I Fe 14 Tea ne RE (20b) 
berühren. 
IH. 


Schließlich sei noch erwähnt, daß man einen Schraubentrieb nicht: nur bei gleichen Para- 
metern und gleicher (I.) oder verschiedener (II.) Windung erhalten kann. 

Die Berührungsüberlegungen zu Bild 4 bleiben erhalten, wenn man a, den Parameter p,, 
a, den Parameter p, zuordnet. 7,T, hat dann zu //, die Steigung (p, — p.)/d, so daß man den 
Satz 2 verallgemeinern kann zu 


Satz 3: Berühren einander zwei Schraubflächen mit parallelen Achsen und 
den Parametern p,, p, so läßt sich längs der Berührungslinie ein gemeinsamer 
Tangentialzylinder anschreiben; die Erzeugenden dieses Zylinders haben zur 
Verbindungsebene der Achsen die Neigung tg» = d/(p, — p.), ihre Normalpro- 
jektion auf diese Ebene ist parallel zu den Achsen. 


In Lit, =p, nl. = —.p. 
Praktisch könnte man z. B. durch die Annahme p, = + 2 p, Schraubentriebe aus ein- und zweigängigen 
Schraubenspindeln herstellen. Für eine kombinierte Grob- und Feineinstellung müßte dann zu 9 = —2p, 


r, = 2r, gewählt werden. Die erste Spindel hat dann etwa den doppelten Durchmesser der zweiten, was be- 
sonders vorteilhaft ist wenn man um sie zwei oder mehrere Spindeln symmetrisch anordnet, eine Möglichkeit, 
auf die bereits in II. hingewiesen wurde. — Zwar nicht als Schraubentrieb, doch als Sonderfall mit 2 #0, 
P2 = 0 läßt sich auch ein Gruppenfräser und ein von ihm gefrästes Gewinde ansehen (vgl. Abb. 299 aus [5]). 
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Resistance on a sphere due to a circular vortex filament 


in an uniform flow of a perfect liquid 
By M.K. Mitra 


Die Arbeit gibt einen Ausdruck (nach Kugelfunktionen entwickelt) für den Widerstand an einer Kugel, 
der bei gleichmäßiger Strömung einer idealen Flüssigkeit durch einen kreisförmigen Wirbelfaden hervor- 
gerufen wird. Zuerst wird mit Hilfe des Butlerschen Kugeltheorems die Stromfunktion gefunden; der 
Widerstand an der Kugel wird sodann mittels der Bernoullischen Druckgleichung berechnet. 


This paper gives a formula for resistance (in series of Legendre functions) on a sphere due to a circular 
vortex filament in an uniform flow of a perfect liquid. The stream function of the motion has been found 
by Butler’s sphere T’heorem and then the resistance on the sphere is calculated by means of Bernoulli’s 
pressure equation. 


L’article donne une expression (dans le developpement de fonctions spheriques) pour la resistance sur une 
boule qui, dans un courant uniforme d’un liquide ideal, est provoquee par un fil de tourbillon circulaire. 
D’abord on trouve les fonctions de courant a l’aide du theor&me spherique de Butler; ensuite la resistance sur 
la boule est Evaluee au moyen de l’quation de pression selon Bernoulli. 


B pa6orTe naetca (B BuNe PAa3losKeHus B PAN IIAPOBbIX PyHRKUMÜ) BhIPa3KeHHe MIA CO- 
IIPOTHBJIEHNHA OKOJIO IIapa, BbI3bIBAeMOTO IIPH PaBHOMEPHOM TeyeHuN NTeaslbHONH HKUNKOCTU 
KpyToo6pasHoü BHXpeBoü HuTbw. DB IMepBylo oyepenb IpNu IHoMomm Teopemsi Byriepa 
onpenenAeTcH PyHKRIHSA TOKA; COUPOTHBIIEHHE OKOJIO MAPA BhIYHCHAeTCcH 3aTeM UPu ITOMomM 
ypaBHeHuf NaBıeHun BepHysuım. 


1. Introduetion 


In this paper, a formula for resistance on a sphere due to a circular vortex filament in an 
uniform flow of a perfect liquid is found. This problem arises in connection with the investigation 
of resistance experienced by a sphere due to its wake, when it is moving uniformly through a 
liquid of small viscosity (such as water). 'The wake in three dimensions has been investigated by 
many and it is concluded that the only quasi-steady motion associated with the wake is that of a 
series of vortex rings at equal distances apart and with their planes perpendicular to the direction 
of uniform flow. But that is also proved to be unstable (Ref. Levy and Forsdyke, 1927) and 
it is concluded that the wake is turbulent (Ref. Jefferys, 1930). But initially there is a for- 
mation of a concentrated vortex which derives its vorticity from the free surface layer. At first 
its strength is small. It then gains its strength and size and ultimately breaks away from the 
free surface layer. The corresponding problem in two dimensions (i. e., resistance on a circular 
cylinder due to two symmetrically situated rectilinear vortices with opposite circulation) has 
been solved by Bickley in 1928. But in three dimensions, no such theoretical attempt has ever 
been made, as far as I remember. 


2. Determination of stream funetion 
A, B intersections of the filament with the meridian plane. 
0 — centre of the sphere — origin 
Ox — direction of uniform flow U at infinity. 
(r, 9), (R, x) — polar coordinates of a point P in the liquid and A respectively. 


2°) 


Due to the uniform flow U at infinity past a-sphere (radius a), the stream function is 
1 a? 1 + a? 
Yo -— Zupli— 2)--35 U r?sin? az 2) > 


when there is a circular vortex filament alone in the liquid and the liquid at infinity is at rest, 
the stream function is 


Ti 
gen (r sin O)Y/2 (R sin y)!!?. 9 


(Ref. Lamb, Hydrodynamics, 6th Edn.), 
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2 2 
where 9= (7) F— x E 


2 zi8 


da B= | VRänta., BZ Bi ee 
0 


ea 3 r Reos(O+y)+R® 
0 


where m = strength of the vortex filament, the sense of rotation of the liquid in the filament 
being clockwise in the upper half of the meridian plane. 


When a sphere (radius a) is introduced in the flow due to a circular vortex filament (the 
liquid at infinity is at rest), the perturbed stream function is given by Butler’s sphere Theorem. 
Ref. Milne-Thompson, Theoretical Hydrodynamics, 1955), as 


ig a? m e : ß 
rn le 0) ee (r sin HJ! (R sin y)t!?- 9°, 
here u .— = Kr)" Be E' 
where 7 pP 


Fe ra ers 
yı — k'2sin?«& 
Ö 


rl2 
Br = [Vi—Rrsirac, 
10] 


ra 4a2r Rsin®sinx 
 a@—2@rRcos(+y)+ TR 


EEE) 


when there is a filament in an uniform flow past a sphere (radius a), the resultant stream function 
is given by 


3 
Penn tm zer sin 0(I 5) + 57, sin Pr (Rein ir) 


where 9 and 9’ have been defined earlier. 


Now 
Reale et leder) 
mh ee en than) 
aa Tr 2 EEE 2 ZB RE ee IE ren 
late) gg ara) 
TE n-3.4 n+ 2 me, 


(the n-th term being deduced from the n-th terms of F and E) 


IE. KON 
NER T8 EEE TE 
i6* ı( u 
fdenoting a hypergeometric function. 
Now 


4r Rsin 6 sin x r r pn 
k? = en r u A ....63 37 nn ii 
r—_drRcsOty)+R 1 sin sin | nt 
(where $ = cos (0 + x)) 
N 


r 


P n EEE 73 r® 
pin dsiny/1+2Pz H (4 ß? Imre 1) 


> ER n* A 5 
+ER-RPIDK+RF-2E+6DT, 4 
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expanding in powers of -(r< R) and writing ehe series up to six terms. 
245 6615 22869 
Br 2 a En 1) 8 10 
BR +5 GR + gl PR Ri 
st [r\312 


=; (ey (sin 0 sin y)972 [ 13 Br +1) 55 Ei (8 B®—4 P) ge (16 —1284.1)., 
Merlin same ' 
f + 32 sin Osiny I +28 Tara a pe (8 B®—4 ß) = + (16 BL—12B2°+ ne 


a Da sin? sin? g Ai +25 da, +84) —. 


245 ı r 
z 2 sin? Osin’ Al +25 + a 
6615 r 22869 : 
+ sint 0 sine gi H2PB a 64 sin sing + | 


writing the series in d to such terms as to obtain an expansion of # in powers of 5 correct to six 
terms. 


a/r\2 _ R r 3\r 35 15 r3 
-5(2) at + 
315 15 6: 105 r 
+5 pr nr "pr mt 3 u De 


5 245 r 
1 +37-4+ (667 +2» 3% Beh en 


225 6615 = gr 
3 2 2,2 we. ee Se 
+ (1Rr + er +! ar >) 
22 
+ (48 B* y + 300 8? y? + 490 B? y? + Spy on vr 30p.} “ i 
— 1? 
EEE y® — 3 5 E= a 
(where y = sin Ö sin x) 
re ne | 2 | a ee 
-2(7) Cinosingralı +88 +39 5 + (pr + Dr + RZ) 
105 325 39 245 15 IN 315 
ati RT re Ne: | ont 
2205 = ne Rh 105 105 525 z E> 
a 34 Fr Zen 2 „2,2 ee 4 . Ka = ! 2 „2 
3465 6615 17325 B 693 22869 
u ET eg LES 2 en De ee 
n en 31) ee - 105 riss 2 | 
2 2. 2 3 yak 2 n ehr ws 
ar, on | ch ya ek 
ei (2) (sin 6 sin x)3'? |1 + 35 cosdcosy + : r 16 cos? 4 — 1) (5 cos®?x — 1) 
Tr 
+ SB N) (7 cos? y— 3 cos x) 
n= Zr (21 cos? d — 14 cos?9 + 1) (21 cos x — 14 cos? x + 1) 
64 Rt 
-)- I Rn (33 cos® 6 — 30 cos? 0 + 5 cos 6) (33 cos? y—30 cos? x-+5cosy) +*--|, 


64 R® 
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BR . = ’ Z r 2 , E 2 F 5 r? ; 
-3 Be ee, er > +55 des Ps(A) R +3 Er PSO) 2 


R 
u: “ Pxlw) a ser = Pin) Ps) I Artur " Ps(u) Pro mtr 


N er 


the general term being inferred from the nature of the first six terms, where 1 = cos 0,1 =c0sx 
and P,„ is Legendre’s polynomial of degree n. 
Similarly, 
a? RN 1 
ne Sr, ne au er i N 312 |__ —— 
v ae 0) 5 5) (sin 0 sin x) S+ SR cos# cosx 
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£E r>R, then expanding k? in powers of R/r, we can get a similar expression for y in 
powers of R/r. 


From the value of , gives by (2), the velocity potential ® is given by 


he Sm. P, n azn+ı ’ 
B=— Urcos0ll + 2,)+  sny Zr + en) Am 0<m o, 
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3. Caleulation of resistance 
. ’ 1 . B . 
Since |Pz(u)| < —n. the series in y and ® are uniformly convergent whenr<R. 


Differentiating the series m, ® term by term with regard to time { and remembering that r 
and 0 are constants with ve = time f, we have 


A ap d : 
n 53 Ber u: al nmsnex 2): Ps) 


contribution by the term | 


[7] 
0 = to the thurst along Ox is given by 


o-—| (e 270 cos 0. sin 0.0, 
Ö 


=—ı@a SE sin: ) 
® ar 
all the other terms vanish, because of the orthogonality relations 


+1 
SF Pla) Pla) du = 0, mn 
2 
Bank when m=n. 

On the surface of the sphere, the radial and crossradial velocity components are given by 

(Odr=a = 0, 

3 mer 2nH+1)a- n 
ei; 1 we 2 BE ET 
(ra =—z Usin6+, sin Osin?z a, En (u) Pi(2) . 


On the surface of the sphere, the resultant speed is given by 


9 ß 3% i on 1 a1 r ’ 2 
= (ort = Urin in a I nme Pa Pa 


3m LEN] ie ER 
2 BEIGE 


Pau) Pa(A) . 


Contribution by the term ( — = 0 g*)to the thurst along Ox is given by 


X, = — |(e-5e@)-2#@- c0s0. sin 0.40 
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ade [0000-500 


alle 3mUasin?x-cosx 
=ieo Im - sin? y Ir | Rn Pz(A) Pa+ı(4) Fi 5% 
(Using e results 
’ _2 n (n — 1) 
EAN 2 ZEN 
Fa 1) Pia) - du = 77 


An u(l — 2) Pun(u) Palu) du =0, whenn#=m+1, 
when SE +1, 
Ta —) Pak) Prrıli) du = 2m (m + 1) (m + 2llam + Dem +3). 
(Ref. Wüittaker and Watson, Modern Analysis, 1927 Edn.)) 
Total-resistance along O x is given by 
Da AR nen sine) 


9 r es 1 azn—ı y ; 
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(R, x) is a singularity of stream function Y, to calculate the velocity of the filament, we omit 
from y the part due to filament itself and write the new stream function as 
adn+ı 


8 a 1 ‚ f 
y = — 5 Ursino(t u 2-3 sin’0ein?y I, nn) . Fu ) Pa(}) (r<R) 


the radial and cross-radial velocity components of A due to y’ are given by 


Vr= - nn = ) U ( — =) cos xXx+ 5 sin? ıd . P,() PaQ); 


B Ze 2. es — 
r? sin 6 Aka ou ve 
1 ay' a3 ) R m. } x ; azn+l j p/: F) 
il ++ re Pe. 


Since the velocity of the vortex (regarded as a line filament) due to itself is undefined, we define 
the velocity of the vortex (regarded as a ring of small cross-section) as 


m a 
7 (1 a2) IOx, 


where 7 = radius of the central filament, 
e = radius of the cross-section. 
(Ref. N.N. Sen, 1920—21). 
at A, the radial and cross-radial velocity components are given by 
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Total resistance along Ox is given by 
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Lastly, I wish to thank Dr. B. S. Roy, for his helpful advice in the preparation of the paper. 
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Die zweidimensionale Fourieranalyse spezieller Koppelkurven” 
Von W. Meyer zur Capellen 


Während bereits die Harmonische Analyse der Bewegungsgleichungen von Getriebegliedern wertvolle Er- 
gebnisse für die Beurteilung der Gelenkgetriebe, vor allem in Hinblick auf die Dynamik, gebracht hat, wird 
jetzt auch die Analyse der Bewegungskomponenten von Koppelkurven durchgeführt, die u. a. für die Unter- 
suchung von Schwerpunktsbahnen wichtig ist. Diese Analyse bedeutet im wesentlichen die Aufstellung der 
Fourier-Reihen für den cosinus und den sinus des Koppelwinkels. Die Koeffizienten folgen aus der Form 
und den Abmessungen des Getriebes und stehen in bestimmter Beziehung zueinander. Einige Sonderfälle 
führen zu besonders einfachen Ergebnissen. 


The advantages of Harmonic Analysis in respect of Kinematic Linkages having been already established 
especially from the dynamical point of view, the analysis is now further extended and applied to the com- 
ponents of coupling rod curves, for example to determine the locus of the centre of mass. T'he analysis to be 
carried out consists essentially in establishing the Fourier Series for\ the cosine and sine values of the coupling 
rod angle. These coefficients follow from the type and dimensions of the linkage and are mutually correlated. 
Some special cases lead to extremely simple solutions. 


L’analyse harmonique des equations de mouvement des elöments de mecanismes & manivelle ayant donne 
d’importants r&sultats dans la critique des quadrilateres articules, du point de vue de la dynamique avant tout, 
on est aujourd’hui en mesure de realiser l’analyse des composantes des courbes decrites par les points situes 
dans le plan de la bielle. Cette analyse est importante dans l’ötude des courbes decrites par le centre de gravite. 

L’essentiel de cette analyse est l’etablissement des series de Fourier pour le cosinus et le sinus de l’angle 
forme par la bielle et le plan fixe. Les coefficients resultent de la forme du mecanisme et de ses dimensions 
et # ont une relation determinee entre eux. Quelques cas speciaus donnent des resultats pariculierement 
sımples. 


T'apmoHnyecknüä aHasıu3 ypaBHeHNuÜ NBU5KeHHUA 3BeEHbeB MEeXAHNU3MOB Mall IEeHHBIe Pe3yJIb- 
TAaTbI A1A OOCY>KNeHUA ITAPHHPHEIX MeXAHN3MOB, TIIABHBIM 006PA30M C TOYKHU 3PeHUA IUHAMUKN. 
Tenepb HPoBoNuTcH Tak’ke U aHalIm3 COCTABAIAIOIIUX IBHFReHUA KPUBbIX CIEeIIeHUA, BA;K- 
HbIH M&K]Iy IIPOYUM HIIA UCCHEeNOBAHUA TPAeKTOPHÜ HEeHTPA TAIKECTH. ITOT aHaIU3 B OCHOBHOM 
CBOAUTCH K COCTABJIeEHHUW PANOB Dypbe AA CUHyca H KocHuHyca yrıaa cuenneHun. Koad- 
$UIMEHTBI HOJIYYAIWTCH U3 BuNa MU PasaMmepoOB MexaHusMA U CBABAHbI ONPeMesIeHHBIMMH COOT- 
HOIIeHHUAMN. HekoTOopbIe YacTHbIe CILyYau NAT 0COÖEHHO HPOCTbIE Pe3yJIbTATEI. 


Daß die harmonische Analyse für die Dynamik nicht nur des zentrischen Schubkurbel- 
triebes, sondern auch jeder Gelenkkette wichtig ist, wurde bereits 1940 am Beispiel der Kurbel- 
schwinge erwähnt [7]. Da hierbei die harmonische Analyse der Abtriebsbewegung besonders 
wichtig in Hinblick auch auf die Bewegungscharakteristik und auf die Getriebesynthese erschien, 
wurde für eine Reihe spezieller ebener und räumlicher Getriebe die harmonische Analyse rechne- 
risch durchgeführt [4], [10], [12], [13], [14], [16—21], aber für den allgemeinen Kurbeltrieb 
(Kurbelschwinge oder Doppelkurbel) wegen der langwierigen rechnerischen Auswertung ein 
Analogiegerät gebaut, das beiläufig auch die Bewegungsgesetze selbst liefert [23], ferner ein im 
wesentlichen elektrisch arbeitendes Gerät, das besonders zur Untersuchung von Mehrgelenk- 
ketten dient [3]. 

Bei dynamischen Untersuchungen benötigt man aber auch die Bewegungskomponenten 
und ihre harmonische Analyse [4], [14], darüber hinaus ist diese noch für die Getriebesynthese 
nützlich, wie bereits Wunderlich für zyklische Kurven zeigte [24]}). 


1. Die zyklischen Kurven 


Gemäß Bild 1 und 2 liest man für den Ortsvektor 5 = A,KmitR,=nRund AK =p 
unmittelbar ab: 


= (RB, - Her per (Hypozykoid), y=-—Wm-—1)a.. .... (la), 
bzw. 
— (R, + R) e* + pe‘? (Epizykloide), valnza la ce resrey Hab), 
so daß in Komponenten oder Parameterform 
DE A COS COS TUR A nel Dee 2a), 
I Ehe er PR LS Te) 
geschrieben werden kann, wenn im ersten Fall A,ı = B,ı = RR RA Bn=P 


m = (n— 1) und im zweiten Fall A,ı = Bı =R + R Am = Bm =p m=(n-+]1) gesetzt 
wird. Gemäß diesen Gln. 2 sind dann x und y in einfacher Weise je durch eine Fourier- Reihe 
dargestellt?). 


*) Herrn Prof. Dr. Alwin Walther, Darmstadt, zu seinem 60. Geburtstag. ”n 

1) Vgl.auch [1] und das Referat in VDI-Z. 92 (1950) S. 870, wenn auch bei de Beauclair der Begriff 
„‚mehrdimensionale Fouriersynthese‘ in anderer Weise verstanden wird. > 2 

2) Für n = 2 ergeben sich im ersten Fall die Ellipsen, für n = lim zweiten Fall die Pascalschen Kurven. 
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Wunderlich [24], [25] überlagert nun den GIn. 2 noch eine weitere „Schwingung“, kommt 
damit zu den „höheren‘‘ Radkurven und benutzt diese zur Synthese, und zwar zur Darstellung 
des „‚Vierecks mit abgerundeten Ecken“ [5]: Die für n=4 und = p/R = 1/(n — 1)=19 
gewonnene Epizykloide?) weist zwar die Form eines Vierecks mit abgerundeten Ecken auf, aber 
die „‚geradlinigen‘ Stücke sind relativ kurz. Und um diese „Geradführung“ zu verbessern, d.h. 
um sich nicht nur auf eine vierpunktige Berührung zwischen Kurve und Tangente zu beschrän- 
ken, wird in Gl.2 noch ein weiteres Glied hinzugenommen. 


Bild 1. Gegenläufiger Umlaufrädertrieb Bild 2. Gleichläufiger Umlaufrädertrieb 


Man könnte infolgedessen in Fortführung dieser Überlegungen jede geschlossene Kurve 
durch eine zweidimensionale Fourier-Analyse bzw. Fourier-Reihe darstellen, d.h. für den 
Vektor eines bewegten Punktes K, der die betrachtete Bahn beschreibt, in Fortführung von Gl. 2 
den Ansatz machen 

seat > ment DU Se 
m=l, 2, ml, 2, 
worin 3, einen festen Ortsvektor bedeutet und r„ die Amplitude der m-ten Kreisschwingung 
4m = Im €”* bedeutet, oder in Komponentenform, wobei jetzt m immer positiv, r„ aber ge- 
gebenenfalls auch negativ sei, 


ven Tr Im cosme, 
y=Y+t & "sinme, 


worin dann x„ bzw. y„ die Fourier-Koeffizienten und x, bzw. y, die konstanten Glieder dar- 
stellen. 


m we An a 


2. Die Koppelkurven 


2.1. Der Ortsvektor 
Der Kurbeltrieb sei gemäß Bild 3 in die komplexe Zahlenebene eingebettet, und es sei A,A 
die Kurbel und ferner A,A=a, AB=c, B,B=b, A,B, = d, Antriebswinkel BA =o, 
Abtriebswinkel eB,B = ß und der Koppelwinkel, d.h. der Winkel, den die Koppelmittellinie 
B A mit der positiven Richtung der x-Achse bildet, gleich y gesetzt. Der Koppelpunkt K habe 
in dem mit der bewegten Koppelebene verbundenen Koordinatensystem die Koordinaten u 


und vo [15] bzw. r und y. Zur Vereinfachung wird unten auch noch ß durch x — ß ersetzt. 
Dann liest man für den Vektor A,K = AA ++ AK sofort ab: 

zm ad Lu? + Ve d— ae H wen hiver 2.2: (da) 
3) Vgl. a. [6] und [22]. 


u ee ee 
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oder 
ee ae en . (4b). 


2.2. Funktionsgleichungen für den Koppelwinkel 


In der vorstehenden Gleichung ist y eine durch die Maße des Getriebes besti i 
. . > t P 
des Kurbelwinkels &. Schreibt man die Gleichung des Getriebes in der Form 


Der ce heß=.], 2 7 soul OR ee (5) 


hin [2, 11] und eliminiert aus den beiden Komponentenformen den Winkel ß, so erhäl: 
£ . 5 ält 
einfachen Umformungen für die genannte Funktion F(a, Ye 0die on Be 


ee . (6a), 
worin 
Zw, Veasın, Rp (6b) 
und 
ef 2200 a (6c) 


Fr x 
Il) 


Bild 3. Beziehungen am Kurbeltrieb mit Koppelpunkt 


die veränderliche Strecke / = AB, darstellt, Bild 3. Man kann aber dieser Gleichung auch eine 
ünmittelbare Deutung geben, wenn man die Winkel = <& A,B,A und ö= X B,AB, ein- 
führt. Es ist dann 

el IA Te BR 


worin ö nach dem Kosinus-Satz, d.h. aus 2 
rer 2 TC = 2 E—DR TC neueren) 
leicht berechnet werden kann, evtl. mit Hilfe eines Nomogrammes [9] und schließlich wird 


asın a Asina 9°) 
ip = —— = ee ee ANENL ONCE 
d—-deos& INCH 


Hierbei stellt aber @ nichts anderes als den Schleifenwinkel einer zugeordneten zentrischen Kur- 
belschleife mit A = a/d dar, wenn man beachtet, daß dieser ja durch Gleiten der Geraden 9 
(Bild 3) durch den Punkt B, entsteht [12], [20]. Dreht man beiläufig A,A um A bei Festhaltung 
von Punkt A und der Geraden g, so hat man einen Schubkurbeltrieb, an dessen geradgeführtem 
Punkt B, eine Schubschwinge B,BA angeschlossen ist. 

In Sonderfällen vereinfachen sich die Gln. 6 oder 7, wie unten nochmals benutzt wird. 
Für b=c hat man eine besondere Vereinfachung, ebenso für «= b und c =d, wobei cos ö 
—=coso, d.h. ö= 9, y = 0 (Parallelkurbeltrieb) oder öd = — 9, y = — 29 (Antiparallelkurbel- 
trieb )ist, in Übereinstimmung mit [16]. 


2.3. Die Funktionsgleichung für den Winkel 
Eliminiert man aus den beiden Komponentenformen von Gl. (5) den Winkel y, so erhält 
man unter Ersatz von ß durch m — ß eine ähnliche Gleichung G(a, ß) = 0 wie oben?): 


P cosß + Q sin  — Se ee), 
4) Auch leicht durch ein Nomogramm darzustellen. 
5) Dies sei hier der Vollständigkeit halber angegeben; beiläufig gilt auch noch 


£=@+b:+2bccos(ö-+%). 
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wobei P und Q aus Gl. (6b) folgen, worin ferner 
T=-(P+P— Ab... .:. urn (8b) 


gesetzt ist und f wieder aus Gl. (6c) folgt. 
Ähnlich wie oben kann geschrieben werden 


BED 
worin © den Winkel AB,B bedeutet (Bild 3) und aus 
PFRZIET AEPER 


ee I ze ee re ; - 
A 3fb Ba (9b) 
folgt. 

2.4. Die harmonische Analyse bei der Koppelkurve 


2.4.1. Allgemeiner Ansatz. Mit den Gleichungen (4) wird man hier im allgemeinen nicht aus- 
kommen, denn in der x-Komponente vermißt man die sinus-Glieder, in der y-Komponente die 
cosinus-Glieder. Man muß infolgedessen zu einem allgemeineren Ansatz greifen, d. h. zur Dar- 
stellung 

z= Ay. +t % A, cosma + 3 B,„sinmo, 


y=Ay + 2 A,ncosma 42 D,„ um, 
Im Ansatz nach Gl. (4) wird somit 


ma... 0 


To 7 Azo» w= Ayo; Im = Ar Ym =, Re und Bin == Ay —0. 


Um nun aus Gl. (10) wieder zu einer vektoriellen Form zu kommen, werde in Anlehnung 
an [14] der x-Achse aus Bild 3 der Einheitsvektor i, der y-Achse der Einheitsvektor j zugeordnet, 
so daß ein Punkt K durch den Vektor 3=1x + jy dargestellt ist. Bildet man diesen aus den 
Gln. (10), so ergibt sich 


TEE 

= ut m—l a... 2 
Dabei stellt 3, einen festen Punkt dar, 3.m cos m« eine lineare Schwingung von der Amplitude 
Zm = |äcm| und der Richtung des Vektors 3.„, entsprechend 3, sin mx eine lineare, um /2 
phasenverschobene Schwingung von der Amplitude 2,7 = |3;m| und bilden beide zusammen 
eine elliptische Schwingung, deren konjugierte Durchmesser 3. und 3m Sind. Setztmana= wtf, 
so sind die Frequenzen der Schwingungen durch m w,/2 x gegeben. Über die Ellipse der Grund- 
schwingung vom Mittelpunkt 3, überlagern sich also die elliptischen Schwingungen höherer 
Frequenz und approximieren die gegebene Kurve®). 


Die Vektoren Z.m» äsm und 3, werden aus den Komponenten gemäß 


3= 50 + 3 dmecosma + 3 3. sinma, 


dcm — i Aym AR } Aym ’ dism = 1 Bm + j Bym ’ jo =1 Ar + Ayo 
gewonnen. 


2.42. Die Komponenten-Koeffizienten. Da hier nur drehfähige Getriebe betrachtet wer- 
den, hat man sich nur auf die Koppelkurven der Kurbelschwinge und der Doppelkurbel zu be- 
schränken, denn nach dem Satz von Roberts [2], [11] kann die Koppelkurve der drehfähigen 

no 2 « ’ » 7 ’@1 ü n ler 1 3 ıry st . a “ u: %L 
Doppelschwinge auch durch zwei Kurbelschwingen erzeugt werden. Schreiben wir Gl. (4a) in 
Komponentenform, d.h. 


T= AC08% +Ucosy—vsiny, 


y=asinx + usiny-+ Dcosy 


ati 


su zeigt der Vergleich mit Gl. (10), daß es genügt, die Werte cos y und sin y harmonisch zu analy- 
sieren, ‚d.h. zu setzen [4], [14]: 


> UA) — Bl J ke = Su » 4 > . 
csy=&=Ast 2% Am c08ma + 3% Bmsinma, (13 
7 I -l Br Su » N E7 7 « « h 3 5, 
siny=n=Ayn+t & A,mcosma+ 3 D,.„sin ma, EN - 
6) Diese Überle a BR rs : R ne j 
) se Überlegungen erinnern an die Bahnen der Himmelskörper. — Da beiläutig eine Ellipse auch mit 


einem Kardankreispaar erzeugt werden kann, vgl. Gl.1 fürn = 2 i i 
Ü kann, . Gl. = 2, d.h. durch Kreisschwingungen i 
Frequenz mit entgegengesetztem Drehsinn, so können auch die einzelnen Ellipsen aus BER IE ee 
schwingungen mit gewissen Phasenverschiebungen entstanden gedacht werden. — Im übrigen wurden Ansätze 
ähnlicher Art auch von O. Tolle in einem Vortrag gemacht, erschienen in VDI-Berichte Bd. 29 (1958), S. 45 
A ‚8.45. 


ih u u 
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Damit schreiben sich die in Gl. (10) angegebenen Koeffizienten 


Ay =UAg—V Ay: Apmu A, dr 0.A;s; 
Au = HU dA Am N, HOT DB 
A A, EA. ml, .. (14), 
Ba UB.nzV Bun: ER Er: Bı=a-tu B.ı+ vB.; 
Bm  UBmt DV Ben Eon 
und es lassen sich damit die Vektoren 39, 3.9 sm 2. B. wieder in komplexer Darstellung angeben: 


de = [Auer A, ee] FL En 
em En Am er 4 A, eva] al en 
m? Ben EB, er En ee 
m=1,2,3,.... Hierbei sind iv, und mw, gleich Null für allem > 1, und fürm =1 gilt 
a ae a ee ai‘ 


2.4.3 Beziehungen zwischen gewissen Koeffizienten. Integriert man die Gl. (6a) in den Gren- 
zena=0unda=2n, d.h. für eine Periode, so ergibt sich, da bekanntlich 


2n 2n 2n 

1 
|. An=— | &eosadı und Bı,= . | sin «as 
[0 0 0 


IT 


1 
9 


De 


ist mit ald =}: 


Au= 


a2 1 da 2 p2 
cd 


Multipliziert man Gl. (6a) weiterhin mit cosm«a und sinma, m=1,2,3,..., und integriert 
in den Grenzen «= 0 und a = 2n, so folgt auch 


BA (Aue DB) SERIEN 


2A Acmaı) a Acmaıl MBamıy Bm) = 0. . 72 (178), 
a 
ee EN ee re EB ai 
DE 0 A nn ana. a ae 
2 Bem— A [Bem+1) + Bem nl + A [dymıy —Anm+y]l=0 . . . (17d), 
a ne 


Es sei noch bemerkt, daß gewisse dieser Koeffizienten auch bei Berechnung des Flächen- 
inhaltes von Koppelkurven eine Rolle spielen [15]. 


2.4.4. Ermittlung der Koeffizienten. Die Fourier-Koeffizienten für & und n können in be- 
sonderen Fällen, wie unten gezeigt wird, leicht rechnerisch gefunden werden.’) In komplizierteren 
Fällen — und das ist bei den Koppelkurven im allgemeinen so — genügt es, die x- und die y-Kom- 
ponente eines Punktes der Koppelmittellinie (u. a. auch des Punktes B) als Funktionen von « 
entweder instrumentell mit einem Analysator oder rechnerisch nach einem Rechenschema zu 
analysieren, da jafür einen solchen Punkt doch = a c0osa-+ u cos y, bzw. y= asina + usin y 
gilt. Ebenso kann man auch die x-Komponente eines Punktes der u-Achse und eines Punk- 
tes der v-Achse benutzen, da ja für den letzteren doch x = a cosa — v sin y wird. Diese Kom- 
ponenten lassen sich instrumentell aber leicht mit einem von E. Lenk [3] entworfenen Gerät 
ermitteln, also auch bei beliebigen Getrieben. 


25. Sonderfälle 

2.5.1. Die Schubkurbel. Hier gehen B, und damit b und d nach unendlich, aber es bleibt 
d—b endlich. Schreibt man in Gl. (7b) ®— b®=(d—b)(d+b)=e(d + b), so erhält man 
beim Grenzübergang in bekannter Weise [14] cosö= e/e — (ac) - cos» und 9=0, daA=ajd 
nach Null geht,. vgl. auch Bild 4, so daß y=6 wird. Damit gilt Ay= eje, Ay = — ale, 
während alle A;„ für m> 1 und alle B;,„ verschwinden — in Übereinstimmung mit GLEZ) 
Die für die Entwicklung von sin d = sin y notwendigen Werte sind an anderer Stelle angegeben 
[4], [14], so daß sich eine Wiederholung erübrigt. Dort [14] wird aber der W inkel & von der 
Parallelen zur Geradführung durch A, aus gerechnet, vgl. Bild 4 sowie Bild 5 mit der Zusammen- 
setzung der einzelnen Vektoren aus [14]. Bei der kinematischen Umkehrung, der geschränkten 


?) Auch bei gleichschenkligen Getrieben treten, wie eine noch unveröffentlichte Arbeit zeigt, merkbare 


Vereinfachungen auf. 
3 * 
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Kurbelschleife [21], geht B und damit c und b nach unendlich und bleibt ce — b= e endlich, 
so daß cos ö = e/d wird und Gl. (16) z. B. auf 2 A,, —A (Ası + Baı) = 2ejd führt. Da hierüber 
an anderer Stelle berichtet wird [21], sei hier, zumal sich langwierige rechnerische Entwicklungen 
ergeben, darauf verzichtet und dafür auf zwei Gruppen von Koppelkurven mit besonders ein- 
fachen Entwicklungen eingegangen. 


Bo @ 


x x 
(R) (RR) 


Bild 4. Schubkurbel, Almessungen Bild 5. Schubkurbel, Koppelkurve 

2.5.2. Der gegenläufige Antiparallelkurbeltrieb [8]. Hier sind wie bei der Schubkurbel die 
Koppelkurven vierter Ordnung und gemäß Bild 6 liest man in ähnlicher Weise, wie für den Win- 
kel £ des Abtriebes gezeigt wurde [16], für den Koppelwinkel unmittelbar oder auch aus Gl. (7b), 
wie oben bereits erwähnt, y=—29 ab. Damit hat man cosy = cos? = co®?p — sin? 
und da, vgl. auch Bild 3, fsinp = asina, fcosp = d—acos«a wird, folgt schließlich 

ET 1—R—21cosa +24? cos?a 
1 1+2—22c0sa H 


Bild 6. Gegenläufiger Antiparallelkurbeltrieb, Abmessungen 


Diese in cos«& = » unecht gebrochene rationale Funktion läßt sich aber auch 
Tue j (1 — /2)2 
| a 
schreiben. Nun wurde für das Übersetzungsverhältnis i, — do/da einer zentri j 

c \ i s S B Be er zentrischen, scl 
Kurbelschleife, d.h. bei A < 1 die Entwicklung [12], 18] R u 


Als — A) ans 
a ee 2 Acosmm . ee 


m=1,2,.. 


1.7: a 
008 = FÜR) —A00sa + 


ui #0 
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gegeben, woraus auch durch Umformung 
1—22 ea 
1+2—21Icosa 


folgt. Setzt man Gl. 20 und die Fourier-Reihe gemäß Gl. 19 in Gl. 18 ein, so ergeben sich für 
die harmonische Entwicklung von cos y die einfachen Koeffizienten 


Au=l—R, Au=—R, Am=(—Ma, m=23, 
und sämtliche B,, =. 
Für siny=—sin29 = —2sing cos kann man 


Ben tern __ 2Asina(l—Acoso) 
7 Er 1+42—21cosa 


schreiben, für (1—Acosa) :(1 +42 —2%Acos«o) aber 1-+i, so daß 
-7=--- hrs227malliit) Sr a (21) 


TERN a (20) 


wird; und setzt man hierin die Entwicklung für i, aus Gl. 19 ein unter Beachtung von 


2sinacosma— sin(m + 1)&—sin(m—1)«, 


so folgt 

een.  ssferner. Bi M22 1) 
und 

Don A uns RE 

Damit schreiben sich die Koeffizienten für die Entwicklungen von x und y: 

An, zu ll 2), Au=vr(l—M). Aı=a—uf, Au=—vR, 

Bel 3) WA ever (Ve), ME. (22)8 

urn onmu.iB, a. Biuuua sn, ) 

DR PA, Sr, wc 


Damit hat der Punkt Omit den Koordinaten 2,—= A,, und yy= A,, die Lage „= r(1— AP) e'?, 
d.h. er liegt auf dem gleichen durch den Winkel y gegenüber dem Steg festgelegten Strahl, auf 
dem auch das dritte Fokalzentrum der Koppelkurve oder der dritte Festpunkt C, bei der drei- 


fachen Erzeugung der Koppelkurve [11] liegt, und es ist AO =r(l—A)=AK-(1—P2). 
Die Grundschwingung liefert, da die Vektoren 


31 =tla— u) —jv®=ia— (tu+jv)? 


und 
31 =11 (2 —72)v + jfa—I (2 — 22) u} 
= (iv—jWIR—R)+ja 
aufeinander senkrecht stehen und nicht die gleichen Beträge haben, eine Ellipse. Fürv = 0, 
Koppelpunkt liegt auf der Koppelmittellinie, hat man eine Ellipse, deren Achsen parallel den 
Koordinatenachsen verlaufen (vgl. a. u.). 
Für die Oberschwingungen ergeben sich aber Kreisschwingungen. Denn es ist 


dcm =" (1 —M)(iu+jo) 
m = A" —M)tv—ju) 


d.h. die Amplituden beider Schwingungen sind für jedes m einander gleich und die Vektoren 
stehen aufeinander senkrecht. Danach gilt auch 


Em = Im cos (ma — y), Ym = Im sin(ma& — y), 
mit 

Im = ir (1 — dr ;U:: u 
Vgl. hierzu Bild 7 für eine beliebige Koppelkurve?°) [8]. 


8) Eine Kontrolle liefert x(0)=u-.a, v(0)=0 und z(r) =u—.a, v(z) = v, wie unmittelbar folgt. 
9) Die Koppelkurve selbst wurde einer früheren Arbeit [8] entnommen. 
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Für die Koppelkurve des ee der Koppel AB gilt mitv = O0 und u=c[2 = d/2 
noch zusätzlich 


a 
Au=s@-M=a—z# und Ba=z#< As 


als Halbachsen der Ellipse für die Grundschwingung'®). 
Bei der Lemniskate [8] wird A = 1/y2 und gilt somit A,, = 3 a/4, B,, = a/4. 


[2} 
Bild 7. Gegenläufiger Antiparallelkurbeltrieb, Koppelkurve 


2.5.3. Der gleichläufige Antiparallelkurbeltrieb"). Auch hier gilt, vgl. Bild 8, wieder 
cosy=c0os2p und siny=—sin2o 
wie soeben. Nur hat man es jetzt bei der zugeordneten Kurbelschleife mit einer umlaufenden 


zu tun, es ist jetztA = a/d größer als eins. Es gelten dann zunächst die gleichen Entwicklungen 
wie oben, aber für i, muß jetzt angesetzt werden [12], [18] 


= —1 3 17 cama=— 1.3. A%008mal ale Be 
m=1,2, m=1,2 
a N 
Für die Entwicklung von & = cos y hat man somit 
Aso — () ; De = - 1 A Aa = (22 - 7 1)/Am u (1 — 12) Am—2 , n= 2‘ 3, ... 


und sämtliche B;„ = 0. Aber fürn = sin y gilt jetzt eine gleichartige Anordnung, d.h. A „= U 
Byı = Az» Bm = u Setzt man diese Koeffizienten in die Entwicklung für x und y ein, so 
folgt schließlich 3, = 0 und 

de1 I (a u/A) -jv/A, 

dsı = iv/A + ja ul); 
diese Vektoren stehen aufeinander senkrecht (das skalare Produkt verschwindet), und ihre Am- 


plituden haben den gleichen Betrag, d.h. die GrundschwinEuge stellt für jeden Koppelpunkt 
eine Kreisschwingung dar, vgl. Bild‘ 9a, deren Radius r, aus 7 = (a — u/A)® e (v/A)? folgt. 


10) Hierbei ist beiläufig r, der Halbparameter der elliptischen Relativpolbahn. 
") Vgl. [8]. 


m 
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Ebenso stellen aber die Oberschwingungen genau wie bei dem gegenläufigen Getriebe 
E Kreisschwingungen dar, Bild 9b, so daß ähnlich wie dort 


Bild 8. Gleichläufiger Antiparallelkurbeltrieb, Abmessungen 


357 


Ja-z) 


Bild 9. Vektordiagramm 


= Tm cos (m +Y), 

Ym = Im 608 (m& + y) 
mit 

Im = I (A? — 1)/Am 
gilt, vgl. a. Bild 1012) für eine beliebige Koppel- 
kurve. Wird v=0 gewählt, so wird un =a—.u/A 
und diese Amplitude verschwindet, wenn u=/a 
— a?/d ist, d.h. es verschwindet dann die Grund- 
schwingung völlig. 

Bei diesen Koppelkurven des gleichläufigen 
Antiparallelkurbeltriebes liegen also, da sich nur 
Kreisschwingungen überlagern, höhere Radkurven 
vor. Wenn man zwei Schwingungen, z.B. nur die 
ersten beiden Glieder der Entwicklung benutzt, so 
liefern diese addiert eine zyklische Kurve als mög- 
liche Näherung. 


2.6. Geschwindigkeit und Beschleunigung 


Setzt man = wtund differentiiert die Rei- 
hen für x und y ein- bzw. zweimal, so erhält man 
die bezogenen Geschwindigkeits- und Beschleu- 
nigungskomponenten [4], [14] 


v,/o = dıldx, v,/o = dyjd« , 
b,/o® = d?x/do?, b,/w? = d’y/da? 


A=9%=:3 
Ya =2 
Wa =15 


Bild 10. Gleichläufiger Antiparallelkurbeltrieb, Koppelkurve 


und hat damit auch die Harmonischen Reihen für diese Komponenten, sofern die Reihen noch 
konvergieren (was bei Getrieben, die nicht in der Verzweigungslage durchschlagen, immer der 
Fall ist). Dabei gilt dann in der Geschwindigkeit A,» = Mm Bm, B,m = — mA, und in der 


2) Vgl. [8]. 


40 ZAMM 39 (1959) Heft 1/2, Seite 40—49 
N N a 


Beschleunigung entsprechend Aym = — m? Ay» Bym = — m? B,, worin die Indices v und b auf 
Geschwindigkeit und Beschleunigung hinweisen sollen. Es können Harmonische, welche in der 
Entwicklung von x und y klein sind, in den Beschleunigungskomponenten groß werden, „Diffe- 
rentiieren rauht auf.“ 
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Shock Wave Propagation in an Infinitely Electrically Conductive 


“ . . - - * 
Gas with Transverse Magnetic Field and Gravitation ” 
by S. I. Pai!) 
Das Verhalten eines eindimensionalen starken Verdichtung 
elektrischer Leitfähigkeit bei Anwesenheit eines transversalen Magnetfeldes und eines Gravitationsfeldes 
wird mittels der Lagrangeschen Methode untersucht Das magnetische Feld bewirkt eine Erhöhung der 
Ausbreitungsgeschwindigkeit sowie eine Verbreiterung des vom Verdichtungsstoß beeinflußten Gebietes. 


} A wesentlichen Punkte des gleichen Problems werden auch vom Bulerschen Standpunkt aus kurz be- 
vandelt. 


sstoßes in einem Gas mit unendlich großer 


The behavior of a one dimensional strong shock wave propagating itself in an infinitely electrically 
conductive gas in the presence of a transverse magnetic field and a field of gravity has been analyzed by 


means of the Lagrangian method. The effect of the magnetic field is to increase the speed of the shock wave 
and the domain influenced by the shock wave. 


*) This research was supported by the United States Air Force through the Air Force Offi ientifi 
Research, Air Research and Development Command Under Contract AF 18 (600) 993. 


!) On sabbatical leave from Institute for Fluid Dynamics and Applied Mathematies, Universi 7 
land. Visiting Professor in Institut für Theoretische Gasdynamik, DVL e.V. for oT 5 ne 


wem. 
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‚ A brief description of the chracteristics of Ihe same problem from the Bulerian point of view is also 
given. 


Au moyen de la methode de Lagrange on analyse Vallure d’une forte onde de choc & une dimension dans 
un gaz @ conductibilite electrique infiniment forte dans la presence d’un champ magnetique transversal et 


d’un champ de gravitation. Le champ magnetique provoque une augmentation de la vitesse de la propagation 
ainsi qu’un elargissement du domaine influence par l’onde de choc. 


En outre les points principaux du m&me problöme sont bri&vement traites aussi du point de vue d’ Euler 
Meronom Jlarpanska u3yyaeTca HoBeneHue ONHOMePHOÜ CHAIBHOÜ YAapHoü BOIHLI B Taze 
c ÖecKOHEeYHO Ö0BINOH IANEKTPOIPOBONUMOCTLE IPH HAJIHYMM TONePEeYHOTO MATHHTHOTO HOAAH 
U NOJA TATOTeHUA. MATHUTHOE HONIEe BbI3bIBAeT YBeiIHyeHHe CKOPOCTH PacıpocTpaHenun u 


pacumpeHue 00JlacTu BIIHAHHA yMapHof BOoAIHbI. ÖCHOBHBIE IYHKTBI TOÜ ;ke 3anaym KpoMe 
TOTO PaCCMaTPUBAWTCH KPaTKo C TOYKH 3peHuA Dürepa. 


1. Introduetion 


We investigate the one dimensional flow produced by a strong shock wave in a gas of in- 
finite conductivity for electricity subjected to a transverse magnetic field and to gravity. The 
shock wave is produced by the motion of a body (e. g. by a piston acting on a column of gas in 
a tube). We shall study only the flow outside the shock region, so that viscous stress and heat 
conduction may be neglected and the electrical conductivity may be considered as infinite; but 
the change of entropy must be taken into account. We also assume that the usual magneto-gas 
dynamic approximations are applicable in our problem [1]. 

There are two ways to investigate this problem: (1) the Eulerian method and (2) the 
Lagrangian method. In the Eulerian method, we may use the method of characteristics. 
In the present case, we have three sets of characteristics [see reference [1], eq. (39)]; i. e., 

(a) the lines of flow 


Fine AR NER 2 0 


where zis the spatial coordinate; {, the time and u, the flow veloecity. 
(b) the characteristiec lines 


eh eh er. ee +, 2): 
(c) the characteristic lines 
dx 
ir ern er, 53 
where a, is the effective sound speed in magneto-gas dynamic, as defined by 
| GSP RT In, HYos dar a ee ENTE: 


Here y is the ratio of specific heats, R is the gas constant, T is the absolute temperature of the 
medium, wu, is the magnetic permeability, H is the transverse magnetic field strength and o is the 
density of the gas. 

The characteristics (1) are double roots of the characteristic equation. Hence along the 
flow lines, there are two relations, viz, (i) entropy S = constant and (ii) H/o = constant. Along 
the characteristic lines (2) and (3), respectively, the following relations hold [2]: 


Jöp 

öt et 
In the case of isentropic flow or in the flow with weak shock waves across which the change 

in S and in (H/o) are negligible, only the two sets of characteristic curves (2) and (3) are needed 
to determine the flow field in the (x, f) plane with the help of relations (5). But as soon as S 
and (H/o) are not constant in the whole plane, we have to use all three sets of lines (1) and (3) 
to determine the flow field. Usually it is difficult to work with three sets of lines in a plane. 
Hence in the case of strong shock waves with a large change in entropy S and in (H/o), it seems 


cu cul 


op! = 
an (u + a,) Fr (u + a,) ern Zu (2). 


| 
rF 
1 


oH oH) 
Me H n Hr (u + a,) ax) 0 dl, | Or) 


ö 


more appropriate to give main attention to lines of flow — u along which the entropy $ 
and (H/o) are constant (except across the shock wave front), instead of using the characteristic 
curves - —=u-+a,. Thus it is more convenient to apply the Lagrangian method of des- 
cription of the flow in the case of strong shock waves. 
2. Lagrangian Equations of Motion [3] 
Let x be the coordinate of an element of volume of gas at time {. We write 


A 2 er er | |: 
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where s is a parameter, characterizing that particular element of volume of gas and thus its path 
in the x-i plane. In general the choice of s is immaterial, but in certain problems, a special 
choice of s will be more convenient than others. We will choose s such that 2=s so long as the 
element of volume was at rest. 

By definition 


de © 
= dt 2 nm a Ei ae Ga BE (7) 
The equation of motion is 
nr 2 ee . (8). 
PT a a nn ; 


The partial derivatives of @ with respect to fin which s remai constant, refer to the changes 
with time observed when following the motion of the elemen‘ of volume of the gas. The left 
hand side term of equation (8) represents the acceleration of the element considered. The three 
right hand side terms represent the accelerations due to the gas pressure gradient, the magnetic 
pressure gradient and the gravitational acceleration g, respectively. 


If we write 


equation (8) becomes 


De +9)=—2@P+pm) - + (10), 


shock wave where pg = —- z 


In the x-t plane (Fig. 1) we shall draw the 
path of the shock wave as x —= X(f). Let o, be the 
density of the gas at rest and o, the density of the 
gas after crossing the shock wave. By the equation 
of continuity, we have 


PIE = Pi Sr ee 
op i 
Ko=S KordXo x Now here ddr = —- ds=®ds and since ds = dr, 
Fig.1. Path of element of hock Poben 
ıg.1. ath of element O as acT0oss a SNOCK wave 
s oD=M-, 2.0. rate, 


where o, may or may not be a constant. In general, one must assume o, to be function of s. 

If we consider only the case where the flow behind the shock wave is continuous so that 
the entropy is a constant along each line of flow, the pressure may be expressed in the following 
form: 


/ 


\Y ®,? 
p=-al)=n(g) en EEE 


where the subscript 1 refers to the state immediately behind the shock front. 
Similarly, we have 


. D, 
H=H = H,| 2): (14). 
Hence the equation of motion (10) becomes 
0p 1) D,\r D,\: 
0; +)-—2 (g) ar 2 (2) | 1) Ve 
. Hi - 
where py, = E Z " = magnetic pressure immediately behind the shock front. 
The velocity of the shock wave is 
= = ni 16 

dt rd Bath nee ee 


At the path of the shock wave we have 
Ash) =s un ne Da 
hence, taking the time derivative along this path, we have 


E—u 
an ee 
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3. Boundary Conditions 

When we consider the problem of a 
shock wave produced by the motion of a body 
(e.g. by a piston acting on a column of gas 
in a tube), we can take s = O0 for thz element 
of volume of the gas immediately adjacent 
to the body. If the motion of the body or 
piston is fully given, (0, t) will be completely 
known (Fig. 2). In other cases, a condition 
oran equation may be prescribed, determining 
the motion of this body. In either case, this 
fixes a boundary condition for the function 
%(s, t) referring to s =. 

A second boundary condition refers to 
the path of the shock wave. It must be 
obtained from the condition (17) together 
with (18) and the generalized Rankine- 
Hugoniot conditions [4] in magneto-gas dy- 
namics. The generalized Rankine-Hugoniot 
relations in the coordinate system in which the 
velocity in front of the shock wave is zero, Fig2. Boundary Condition in the Flow of Shock Ware 
are as follows: EL 


gE=a(l— u) 
_Pı 1 ah (1 ei Po Als Pr, a, 


y® € o& 0, CH 306 & 

HE=H, (—U,) (19), 
£ Pass! Be 1  _Po Pr ee! BR: -. 
EN A an a TER Te 

Da, HN 
—. Pal 7 


where the subscript 0 refers to value immediately in front of the shock front and the subscript 1 
refers to value immediately behind the shock front. 
From equations (19), we haye 


2 


22 Ne > Dre Ep 
„-| 1 743% % Bm one +8 ram + ne | 


y + 1 2 £ € 2 & 08 0 
(29). 
Furthermore we have h 
BP FE Fr Ur - DE RI REIT IT 
and 
ER 27 a EN En Eh RAS ie 


4. Solution by Series Expansion 
For small values of s and f, we may use the method of series expansion to find o(s, £). 
We write 


oshHd=t+ zit t+bhts+b, 2 + +ots+tat? cs +--- (23). 
The path of the shock wave may be written as 
Des hs RR 2A): 
Here a,, az; b,, Da, ba; €; CC» 04; --„ Ip Ta Tz... are unknown constant coefficient to be deter- 
mined. 
If the motion of the body is known, we have a condition at s = 0, which takes the form 
DIE ee P-Lı . 220.2. 2200.2. (5) 
from which «a, b,, C,, . .. can be found. 


For instance, if we take the case of uniform motion with velocity V, so that (0, t) = Vt, 


then 
di ==iV. , EA Dee Eee re N. 


These values will be accepted in the following lines. 
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Along the path of the shock wave, we have (5,1) = s. Thus with equation (26), we have 
se V(ns+r82+.)+,s+bn24+. +4 ve SR 
Hence 
„el—-Vn, betdhrhrn) u - —-(V/n+bn+6n + GN). 120 
From equation (24), we have 
Pe 


= en+t+2ns+31n2+° EA es a a I a a 29), 
€ ds 
EIEELT ) 
n(1+ an (30). 
ber. 2r, n__gls (12 12 gEls 2 
= 1- sta 3 32) L ei + 2) er ) 


Substituting equations (29) and (30) into equation (20) and expanding the expression into a 
power series of s, we may write formally the expression of u, as follows: 


j 
u u PART RETTET 


where A,, A,, etc. are constant coefficients depending on the exact expressions for both py, 
and o,. A special case will be worked out in section 6. 
By definition 
u- 8 _ VınsIBtsuRten Mn Mo 
and thus 
„u=V/+bs+2,(nstn2?+:'.)s+5g2+*-» 
=V+bs+Q@con+c)®?+:::- ee 
= VYrbshetsr er 
where * =2cr, + tz. 
Equating coefficients in equations (31) and (33), we have 
A, 3 A 
Bi — te 
a » a 2 he EEE: 
Since the expressions A,, A,, etc. contain the coefficients r,, r, etc., we may find relations for 
Tj, T, etc. in terms of b,, C, C, ... and the parameters y, a, A, V ete. 
Finally we determine the coefficients b,, c,, €, etc. from the equation of motion. 
We are going to calculate the coefficients for the case of a shock wave in a gas of con- 
stant temperature, subjected to a constant magnetic field strength H,, and having a density 
decreasing in consequence of gravity. The corresponding case without magnetic field has been 


worked out by Burgers in reference [4.] We shall compare our results with Burgers results 
in order to show the effect of the magnetic field on this problem. 


5. Shock Wave in an Ionized Gas of Constant Temperature and Subjeeted to Gravity 
and under Constant Transverse Magnetie Field 


We consider a vertical column of gas of constant temperature and take x as the vertical 
distance, positive upward. We have the following conditions for the gas at rest. 


pressure . DH) =ARRT,: ) 
density 2 al 2)" 
temperature = T, = constant, 
magnetic field = H, = constant, ee 
N * / m 
sound speed = m, =yYyRT,= constant, 
dp, 
= —0o( 
dx 003 
The equation of motion (8) in the present case reduces to 
dpo 


a ee n, 


A a u = 
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Substituting from equations (35) into equation (36) and noting that s = x for the gas at rest, 
we have 


d y 
n = a Se er (37) 
Integration of equation (37) gives 
gs 
RU EN EEE ER a 1a 00 


where 0% is a constant. 
Substituting equations (18), (21), (22), (35) and (38) into equation (15), we have 


Ze. je 2 (« )z Dre oe 

za t9 ds ae ram a Den 
u [7] 

where d, = Be and =. 


Equation (39) is the equation of motion for the flow produced by a strong shock in a gas 
of constant temperature and constant magnetic field, subjected to gravity. 


The boundary conditions are: 
(1) a given condition for the path s=0e.g. 9(0,) = Vt, 


(2) g=s along the trace of the shock wave together with the generalized Rankine- 
Hugoniot relations across the shock. 


6. Numerical Example 


We are going to determine the flow field behind the shock wave in an infinitely conductive 
gas of constant temperature and constant magnetic field subjected to gravity with the following 
data?) 


Tr h) Px: 
v=03, dsl, =—, — =01. 
; 5 3 000 
From equation (20) with the above numerical values, we have 
7 rei Ki 11% 2 ST mil: ven sine 
3[7, "1 ef 2]. a 
Hi 2: er Eis newer > 
Also 
Bet 5 23 125 
Ber | eo Se ed EIER 4 
e IrZgshVs Ele et tet: 


Substituting equations (29), (30), (33) and (41) into equation (40), we have 
0.30 +5b,s +c* +... 


317 ae De 
-. = S 2K, + KK = 42 
(ER RRl— 5 Kıt2Ke+z Kıkı) + 5 Rı) «2 
with 
2 AT ar, 2 05... 
Eee a Bu DE 


R,=[n +2ns-+37,S°+---, 
> 25 

Bear 

Equating the coefficients of equal powers of s on both sides of equation (42), we have for the 
zeroth power of s 


Rn 


ee 50T 7 \2 
97 0 Auge 47. 6 0.8) 
or 
ae ER =0O...2.. rl). 


=) The numerical data, except PHo/Qoo, are the same as those given in reference [4] by Professor Bur- 
gers. In this way, we may compare the results with and without magnetic field effects. 
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Equation (43) has three roots which are 
ner 7 00088 
= 10001 em re 
Tj3 = 3.4260 
We must discard the root rj, because, according to equation (24) would give a negative value 
of t for small s which is outside the region { > 0 in which we are interested. Similarly we must 
discard the root r), because according to equation (28) it would give a negative of a,, and thus 


in equation (23) would make x = 9(s, !) negative for small sand {. Hence the only root which 
we can use is 7, — 7, — 0.7588. We shall use this value for all the following calculations. Thus 


1 


r, = 0.7588, a 1317. um Er Ben. nn 
1 _ 
From the first equation of (28) we then obtain 
4 el—03n = 0.7722, er (46). 


From the coefficients of the terms with the first power of s in equation (42) and the value 
of r, given in equation (45), we have 
1 = — 0.2307 55 — 044310 2 WE a EEE (47). 
From the coefficients of the terms with s? in equation (42) and the values of r, and r, given 
in equations (45) and (47), respectively, we have 
T3 = 0.05650 b + 0.03545 9 d; — 0.01910 E — 0.1538 c* . .. . .. - (48). 
Substituting the values of r,, r,, and r, into equation (30), we arrive at 
E = 1.3179 + (0.8013 b, + 0.1557 g) s 
+ (0.1929 53 + 0.0045 g 5, + 0.1179 2 + 0.8013 cH) + --- . . . . (49). 
We further have 
u, = 0ISW ER Ss erst Fun: re (50) 
and consequently 
Eu, = 0.3954 + (1.5583 db, + 0.046719) s 


+ (1.5583 c* + 0.8592 b3 + 0.1571 gb; + 0.0357 9)2+--- ... GH). 
Next from (23) and (24) we find: 
ir ) 


=— o Pr = 2 S - 0... 
28 4 +b,t+2b,5 4 


= 07724 +b,s+2b,s+:--- 
D, = 0.7724 + (0.758856, + 25b)S + . 2... 2 20 (53). 


Substituting the expressions (23) and (51) to (53) into the equation of motion (39) and 
equating coefficients of the same powers of s and tf, the following results are obtained: 
For the zeroth powers of s and f, we derive 


3.4067 bz —= 0.6123 q 


or 


bb = O1707 0. 0 une Re Re 
Substituting equation (54) into equation (47) we find 
a BE. re rn a En 


From equation (28) we have 
b; Vrtbr) = 0.1150 Se 
For the terms with the first power of Zi, we derive 
— 0.2551 9? + 4.8719 c, = 0 
or 


Cs = V.OBBBE RE en ne 
“rom the terms with the first power of s, with the help of (57), we obtain 
= — 0.012503 GERT ee 
Hence 
et = 27,6; +0; = 0.033341 GE (59). 


Wr 
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Equation (48) then gives 
T3, = —.0.01603 9 (60) 
Equation (28) gives 
4a=-  (Vn Fb Ion Fan) - 00208 2... .. (61). 


The following table gives the numerical values of the coefficients for the flow field with 


and without magnetic field (Burgers values) 


nr Without With 
Coefficienta Magnetic Field Magnetic Field 

FE 0.3 0.3 
8 0.7541 0.7724 
E I 0 0 
6: 0.1900 g 0.1797 g 
b, sr 0.1105 9 
& 0 0 
> 0.0118 g2 _0.01253 g? 
& 0.06297 g? 0.05236 g2 
& _0.03243 92 —0.01220 g? 
c* 0.0437 g? 0.03334 92 
a 0.8198 0.1588 

5 0.0509 g 0.086830 9 
= 0.005829 | —-0.01603 g 


Summar:zing we obtain for the case with magnetic field (Fe — 0.1) 
000 


9 = 0.3000 £ + 0.7724 s + 0.1797 gs t— 0.1105 9 
— 0.01253 9 s ? + 0.059236 2 21— 0.0120 2s°-+---. 
Along the path of the shock wave, we have 
t = 0.7588 s — 0.08630 g s® — 0.01603 95? +» - - 
u, = 0.3000 + 0.1797 gs + 0.033341 9° + - 
& = 1.3178 + 0.2998 gs + 0.1517 8° + --- . . 
We further have 


& = 0.7724 + 0.1797 gt — 0.2210 gs — 0.01253 9? 2 + 0.1047 9? st — 0.03660 ! 5? -+... 


d z 
PO 
.G # 
50 OA 
„e 
er 
NIS nee 

40 RZ m mag &: 
8 Y w'! -nock wave Di 
S 7” 2 ir 
& / 2Z 
2 30 7 ./ 
Z ER 
= / ? 
» De ---- Portion of series development 
& L RE 
Z > which is not accurafe 

UEORER—.n. — Portion obtained by numerical 


integration 
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Distance (5) Unit length =350. m. 


0 20 +0 


Fig. 3. Flowfield of our numerical example 
v— 5/8, pP a,[%s = 0,1, a, = 350 meters/sec., g = 9,81 meters/sec?. 


ED 


le, 
> (68), 
....(65). 


(66). 


We have plotted the path of the shock wave from equation (63), with three terms on, in 
Fig. 3. The path of motion of the body (s = 0.31) and that of the shock wave without magnetic 
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field are shown also in Fig. 3. The flow field between the path of the body and that of the shock 
wave are what we are interested in. It is evident that the results of our series expansion are not 
valid for large s, because at large values of s, the time appears to go back ward. 

A way to determine the range of validity of our series expansion up to three terms is as 
follows: 

We calculate the flow velocity u, immediately behind the shock wave front by two different 
methods: one is obtained from equation (64) which may be denoted by u,,; and the other, which 
may be denoted by u. is calculated from equation (20) with the value of £ given by (65), if the 
results of series expansion are accurate enough, these two values of ız,, i. e., u,, and u,, should 
be equal. On the other hand, if the results of the series expansion are not accurate enough, 
discrepancies between ı,, and u,, will occur. We may use the expression 


—u \ 
Asa, ee 
U. 
as a criterion to determine the range of validity of the results of the series of expansion. 


We found that for small values of { und s, say below f < 10, A is negligibly small. As s and f 
increase, / increases too. For instance 


t=15 s=2198 20 

= 20,897, Ba 0E 

T=80, = rel 
Probably we may conclude that for £<20 and s< 30, the results of the series expansion are 
accurate up to 10%. Above {= 20 and s = 30, the results of series expansion is poor, and we 


should use a numerical integration method to calculate the flow field, starting from the results 
of the series expansion. 


7. Numerieal Continuation of a Solution 


We will assume that the solution is known by the series expansion method until the line AB 
along which ? = constant (Fig. 4). We then know the values of, u, p, o and H at all points on the 


XS 


ig, 4. Method of numerical continuation 


line AB, and thus we can calculate the derivatives of p and H with respect to s at all interior 
points of the line AB. By equation (10), we can calculate du/dt in these points. As a result, we 
can determine u in all points from B’ to A’, e. 9. a: 


ou 
u = u + (55), At Bann 8, 0% ln Far uele wo re EEE 


The end conditions remain to be considered. 
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From equation (13) we have 
erg Esel Mi 
ot ot ® öl 8 \® ER : 

From equation (14) we have 


os 8% 


a Fe ee a] a wahr 
At the end A, we suppose that (0, £) is given. This will fix u at B.. | 
The conditions at point B, may be determined approximately by the following manner: 
From the value of & at B and the interval At, we may determine the location of B, which 
is approximately the location of shock on the line £ = t,. From the values of u and du/ös at B’, 
we may determine the value of u at B, which is approximately u, on the line {= t,. From this 
value of ız, and equation (20) we determine the value £ at B,. After the value of & at B, is known, 
we may proceed further in a similar way. 


We have applied this numerical integration method and found that when { < 10, the results 
of the numerical integration are practically the same as those of series expansion. As fincreases, 
the errors of the results of series expansion increases. Usually the error in the values of pressure 
is larger than the corresponding error in velocity. When equation (67) is used to check the results 
of numerical integration, the error is negligibly small. 

We started from {= 20 and found that the results of the numerical integration differ 
considerably from those of the series expansion. The new results are shown in Fig. 3 too. We 
find that the shock wave velocity obtained by the numerical integration is smaller than that 
found by series expansion. Since the larger value of the shock wave speed given by series ex- 
pansion leads to inaccurate results, our numerical integration must be considered to be more 
reliable. 


8. Conelusions 


The propagation of the one dimensional flow produced by a strong shock in an infinitely 
conductive gas subjected to gravity and to a transverse magnetic field has been analyzed by means 
of the Lagrangian method. At first the method of series expansion has been worked outin detail. 
This method is applied to the case of shock wave in gas of constant temperature and under 
constant transverse magnetic field. The effect of the magnetic field is to increase the speed of 
shock wave in comparison with the corresponding case without magnetic field. Hence for the 
same motion ofthe body, the range of the f!ow field influenced by the shock wave with a magnetic 
field, is larger than that in the corresponding case without magnetic field. 

A simple criterion is suggested to determine the range of validity of the results of the series 
expansion (67). It was found that the series expansion is valid in a rather large range. For 
t = 20, the error is about 10%. Beyond the range of the validity a method of numerical inte- 
gration is given to calculate the flow field. The same criterion (67) is applied to the results of 
numerical integration. It is found that the error is negligibly small. 

The characteristic of this problem in Eulerian method of description is also briefly described. 


Acknowledgement 


The author wishes to express his thanks to Professor J. M. Burgers for his inspiring dis- 
cussions and his appreciation to Mr. A. Speth for his work in numerical caleulations of this report. 


References 


[1] 8. I. Pai, One dimensional unsteady flow of magneto-gas dynamics, Tech. Note BN-86, AFOSR-TN-56-537, 
AD 110356, Institute for Fluid Dynamics & Applied Mathematics, University of Maryland, 1956. 

[2] R. E. Longhhead, Solution of problems involving in hydromagnetic flow of compressible ionized fluids, 
Phys. Rev. 99, No. 6, September 15 (1955). ' 

[3] The method applied in section 2, the series expansion in section 4 were given for the case without a magnetic 
field by J. M. Burgers in ‚‚Non-Uniform propagation of Shock Waves“, Lecture series No. 10, prepared 
by 8. I. Pai, Institute for Fluid Dynamics and Applied Mathematics, University of Maryland, 1951. They 
have been extended here to the case where a magnetic field is present, with some alterations at a few place. 
A different method for numerical continuation is used in this report from that of J. M. Burgers. 

4] 8. I. Pai, On exact solution of one dimensional flow equations of magneto-gas dynamics, Tech. Note BN-82, 
AFOSR-TN-56-486, AD 97370, Institute for Fluid Dynamics & Applied Mathematics, University of Mary- 


land (1956). 
Manuskripteingang: 10. 3. 1958 


Anschrift: 8. I. Pai, Institute for Fluid Dynamics and Applied Mathematics, University of Maryland, Maryland 


(USA 
4 


50 . ZAMM 39 (1959) Heft 1/2, Seite 50—67 


Asymptotische Gesetze der turbulenten Ausbreitung 


von Heißluftstrahlen in bewegter und ruhender Außenluft 
| Von W. Szablewski 


Auf der Grundlage des Prandtlschen Schubspannungsansatzes für freie Turbulenz werden die @e- 
schoin tigkeit und ee runder und ebener Heißluftstrahlen in bewegter und ruhender Außen- 
luft für große axiale Abstände von der Düsenmündung explizit mittels der Ähnlichkeitshypothese berechnet. 
Der Vergleich mit Messungen runder Heißluftstrahlen ergibt befriedigende Übereinstimmung. Besonders 
hervorzuheben ist dabei das Ergebnis, daß sich der empirische Koeffizient x, der Theorie für das Modell 
eines runden Düsenstrahls von der Mündung bis zur Asymptote praktisch als eine Konstante erweist. 

The Prandtl expression for the moved shearing stress for free turbulence is used to compute according 
to similarity, the velocity and temperature fields of hot air jets both for circular and for slit-like cross-section 
for large axial distances from the nozzle. The results obtained for circular eross-section agree quite well with 
the experimental results. Special mention deserves the fact that the empirical coefficient #1, appearing in 
the theory, for the jet with circular cross-section is practically constant along the line from nozzle to infinity. 


Sur la base de la disposition de frottement turbulent selon Prandtl pour la turbulence libre les champs 
de vitesse et de temperature de jets circulaires ou plans d’air chaud dans lair exterieur agite ou calme sont 
calcul&s de fagon explicite pour de grandes distances axiales de la bouche d’une tuyere au moyen de U’'hypo- 
these analogique. La comparaison avec des mesures de jet circulaires d’air chaud a pour resultat un accord 
satisfaisant. Il y faut souligner avant tout le resultat disant que pour le modele d’un jet circulaire de tuyere 
de la bouche jusqu’a l’asymptote le coefficient #1 empirique de la th£orie se trouve ätre pratiquement une 
constante. 

Ha ocHoBaHuu Ipmema Ilpanıtıa, Ucxonnmmero U3 HAUPA;KEeHHUA CABHTA, UIH CBOOOAHON 
TYPÖYJIeHTHOCTU BbIYUCHAIOTCA MOAA CKOPOCTel U TeMHepaTypbI KpyT!IkIX H ILIOCKUX CTpyü 
ropAuero BO3AyXa B ABMSKYINelcH HM IIOKOMINeÄcH BHeIIHeÜ BO3AyIIHOÄ cpene MIA 6OABIUHX 
OCeBbIX PACCTOAHMÜ OT BbIXONA COIMIA; 3TO OCYINECTBIAETCH IPH IIOMOIIM TeOPHH IONOÖHA. 
CpaBHeHnHe C U3MepeHUAMH HA KPyTlIbIX CTPYAX TOPAUyeTO BO3AyX AAaeT YAOBJIETBOPHTEJIBHLIE 
Pe3yJIbTaTbI OTHOCHTEJIBHO COTAJACHA C ONBITOM. OTYacTu CJIeNyeT HOoNYePKHyTb, 4TO B pe- 
3yIIBTaTe HCcCcHeNOBAHHÜ 3MIIHPHYecKUh KoePPHlMeHT xı TeOPHH MIA MoNelm Kpyraofi Co- 
MA0BOH CTPyM OT BbIX0Na COMIA MO ACHMITOTBI OKA3bIBAeTCHA IPAKTHYeCKH NOCTOAHHOÄ Be- 
AHYuHoü. 


I. Einleitung 


Untersucht werden die durch turbulente Mischung mit der umgebenden Luft erzeugten 
Geschwindigkeits- und Temperaturfelder von Heißluftstrahlen in sehr großem axialem Abstand 
von der Düsenmündung. Den Ausgang der Untersuchung bilden die vom Verf.!) aufgestellten 
Gleichungen zur Berechnung turbulenter Strömungen von Gasen stark veränderlicher Dichte, 
die sich hier — wenigstens für die Bewegungsgleichungen in der Richtung der Strahlachse — 
auf die bekannten Gleichungen für Felder konstanter Dichte reduzieren, für welche die asympto- 
tischen Gesetze im wesentlichen schon ermittelt worden sind?)®)*). Demgemäß werden wir hier 
die schon bekannten funktionalen Formen erhalten, welche das axiale Verhalten und die Profile 
des Geschwindigkeits- und Temperaturfeldes beschreiben. Dennoch ist es nicht so, daß man 
sich mit dieser Feststellung begnügen könnte, denn die in diesen Funktionen auftretenden Bei- 
werte, die mittels Heranziehung der Erhaltungssätze des Impulses und der Energie bestimmbar 
sind, zeigen eine wesentliche Abhängigkeit von der Temperaturdifferenz zwischen dem aus- 
tretenden Heißluftstrahl und der umgebenden Luft. Dieses Ergebnis war insofern von vornherein 
zu erwarten, als bei der Vernachlässigung molekularer Wirkungen das Feld auch in der Asymptote 
als eine Mischung von individuellen heißeren und kälteren bzw. leichteren und schwereren Par- 
tikeln erscheint. 


Die weitere Untersuchung fußt dann auf der Ähnlichkeitshypothese, d.h. auf der An- 
nahme, daß die Profile der Geschwindigkeits- bzw. Temperaturverteilung über den Schnitten 
senkrecht zur Strahlachse zueinander affin sind. Diese Hypothese wird dadurch gerechtfertigt, 
daß für große axiale Abstände der Einfluß der charakteristischen Länge des Düsendurchmessers 
oder der Düsenspaltbreite als vernachlässigbar anzusehen ist: in der Asymptote erscheinen die 
Düsenmündungen als Punkt oder Linie. 


Der in Abschnitt V für die Asymptote der turbulenten Ausbreitung rotationssymmetrischer 
Heißluftstrahlen durchgeführte Vergleich mit dem Experiment zeigt gute Übereinstimmung. 
Von grundsätzlicher Bedeutung ist dabei, daß der sich für den empirischen Koeffizienten x 
der Theorie ergebende Zahlenwert nur wenig verschieden ist von dem Wert, den die Unter- 
suchung?) der Mischzone turbulenter ebener Heißluftstrahlen, die zunächst durch eine Scheide- 


1) W.Szablewski, Ing.-Arch. 20 (1952), S. 67. 

°®) W. Tollmien, Z. angew. Math. Mech. 6 (1926), S. 468. 
®) H. Görtler, Z. angew. Math. Mech. 22 (1942), S. 244. 
4) H.Schlichting, Grenzschichttheorie. Karlsruhe 1951. 
5) W.Szablewski, Ing.-Arch. 25 (1957), S. 10. 
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wand getrennt sind und sich dann mischen, ergab, wobei zu bemerken ist, daß das letztere Modell 
das Feld an der Düsenmündung beschreibt. Damit wird erhärtet (vgl. Abschnitt IV), daß der 
Prandtlsche Ansatz für den Austausch von Strömungen der freien Turbulenz als sachgemäß 
anzusehen ist. 


II. Bewegungsgleichungen 
Wir zeigen zunächst, daß die Bewegungsgleichungen in der Hauptströmungsrichtung sich 
für sehr große axiale Abstände von der Düsenmündung auf die Gleichungen für Felder konstanter 
Dichte reduzieren®). 
Im rotationssymmetrischen Fall (runder Düsenstrahl) wird das turbulente, im zeitlichen 
Mittel stationäre, Mischungsfeld bei Vernachlässsigung von Gliedern kleinerer Größenordnung”) 
durch folgende Gleichungen?) beschrieben: 


Impuls 
EEE ner 2 kA OR: on _.00 
tr ala tEra,,) — 
Energie 2e7 
oeelu croc os 2 Le 
Er + 2r — E eo) Ba ee (II), 
Masse 
orou  OroDv do 
a Te 2 Kb 


——_—_ Ei 2 
3) 


2 
— 
Ug | 9 
Er 


| 
| 
a en Se ex 
| ucz) ca) 
| 
I 


Bild 1 


Hierbei bedeuten x, r Zylinderkoordinaten, u bzw. v die Geschwindigkeitskomponenten in der 
x- bzw. r-Richtung, o die Dichte, T die absolute Temperatur, c, die spezifische Wärme bei kon- 
stantem Druck; E ist eine empirische Zahl, die das Übertragungsverhältnis von substantiellen 
Eigenschaften und der (Druckschwankungen während des Transports ausgesetzten) Geschwindig- 
keit beim turbulenten Austausch angibt. Das Überstreichen bedeutet zeitliche Mittelwerte; im 
folgenden lassen wir die Querstriche fort. 

Die scheinbare kinematische Zähigkeit wird hier durch 

&(2) = x, b1(®) (u(&) — U) RE A ee N) 

beschrieben, wo x, ein empirischer Koeffizient, b,(x) die (halbe) Mischbreite des Geschwindig- 
keitsfeldes, u,(x) die Geschwindigkeit des Düsenstrahls längs der Strahlachse bedeuten; u, bzw. 
u, bezeichnen die Geschwindigkeit und T, bzw. T, die Temperatur des austretenden Düsenstrahls 
bzw. der umgebenden Luft. 

Die Größe c, wird als Konstante angesehen. n 

Beachtet man (III), so können wir bei Einführung der Übertemperatur 9= T—T, 
das Gleichungssystem (D—(III) auch schreiben: 
Impuls 


‚oo du 


e ou a ou @| ö [ 2) 
Bo.U el 
ee Sr är\ dr 
2. 6) Diese Feststellung ist nicht etwa trivial. Zwar wird in der Asymptote per definitionem die Differenz 
der Dichten von Strahl und Umgebung klein. Dasselbe gilt aber auch für die Differenz der Geschwindigkeiten 
von Strahl und Umgebung, so daß es von vornherein nicht feststeht, ob der Temperatureffekt in der Asymptote 
vernachlässigbar ist. 

?) Dabei setzen wir voraus, daß die äußere Luftströmung unter keinem Druckgradienten steht. 

8) Vgl. z.B. W.Szablewski, Ing.-Arch. 26 (1958), S. 358. 
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Energie 

0 od o[ 9 
rgun. +reu -Ee@)lei[r%) 
Masse | 


Oxzu Orv 
— Eee... Wa Et WERE VIEL 
OX a or y | ur, 

Gegenüber dem Feld konstanter Dichte treten die unterstrichenen Glieder zusätzlich auf. 


Für sehr große axiale Abstände von der Düsenmündung haben wir nun 0 X 0, (Dichte 
der umgebenden Luft), u x u, und damit sind öo/ör sowie du/ör als kleine Größen anzusehen. 
In der Asymptote sind demgemäß die unterstrichenen Glieder in Gl. (V) und Gl. (VI), die den 
Temperatureffekt enthalten, offenbar vernachlässigbar. - 

In entsprechender Weise ergibt sich der gleiche Sachverhalt für ebene Strömungen. 


III. Asymptotische Felder der Geschwindigkeit und Temperatur 
Wir behandeln zunächst den Fall runder Düsenstrahlen. 


A. Rotationssymmetrische Strömung 
Nach Abschn. II lauten die Bewegungsgleichungen: 


Impuls 
ou ou o/ ou 
(u in -H D är) Pe &(X) dr [ =) 2% ee er er un (IL, 
Energie 
[70 0» 0/ 
| reale) un (IX), 
Masse 
u 0 
BR FÜLLEN NE (X). 


Wir werden hier die Fälle bewegter und ruhender Außenluft zu unterscheiden haben. 


1. Bewegte Außenluft (u, # 0) 
Wir wollen die anzuwendende Methode an Hand dieses Falles, der auch für Felder kon- 
stanter Dichte noch nicht behandelt worden zu sein scheint, ausführlich darstellen. 


Geschwindigkeitsfeld: 


Unter Einführung der dimensionslosen Größen 
sr Hrn), fkehluriee D 0% 
ie + —, De und = me n—7Z 


machen wir den Ähnlichkeitsansatz: 


u—u, ua) —ıu T/ro v 'x\P 

——— — N), = r, Düsenradius),; —— —[ : 

UV— U, w—Uu em) l (z/r,)" (aD ipeirainn); W— U \r ) von); 

y an / m» \m rn (1.1). 
wir setzen dabei ein: — u(*) x a = B(&) 
I — U, To In Äro; 
Mit den 'Transformationsformeln 
1) = l [7 [ | eh n ö| 
Or |Tolein,  (&ro)* Onlam OK Tori, a: Oc/roln zT, " nern 


erhalten wir dann mit der asymptotischen Beziehung u x u, für Gl. (VII) 


’ 


a2 fmpenn 2) +(2] m ep 1d 
Te ieh um Fe v2 -(2) But) 


und für Gl. (X) 


Re KR; dop\ ac\® / day‘ 
=) U (m np—nn? an) | Iv +n 2) 0, abet (1.3). 


To To 
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Vergleich der Potenzen von x/r, ergibt 
m—l=p+m—n=2m—n*) und m+n-—1 == pn, 
Eine weitere Gleichung für die Exponenten liefert der Impulsintegralsatz 


12 
[ ouw—u)rd=o,u(w-—u,) 7 RR 7 R (1.4). 
ö 
Für die Asymptote mit 
4 01 % 
RU, NmSaelenn 
5 0% T, 
sowie dem Ähnlichkeitsansatz Gl. (1.1) ergibt er 
20, j ia, 1 
— [ei dn = ee rl a a 1.5), 
5 R a u,1+ (/T;}) m 
also m=—2n. 
Physikalisch sinnvoll ist nur eine Lösung mit n > 0 sowie m<0,p<0. Man erhält eindeutig 
1 2 4 
men le Dar Re ee (1.6) 


Die demgemäß aus Gl. (1.2) folgende Gleichung 


dp —u (n$ dp r 
2 4 1 
2nY een 34, BU” r A 
mit den Randbedingungen (0) = 1, p(®) = 0 hat die 
1 
= en i 2 = ———— ..... 0.0... 1.73: 
o=e mit o a (1.7) 
die aus Gl. (1.3) folgende Gleichung 
U do 2) = : dıy 
= (ene-+w9 ale 
mit der Randbedingung (0) = 0 wird gelöst durch 
ay= = RE a Eule (1.8). 


Für den Streckungsfaktor o gewinnen wir noch eine weitere wichtige Beziehung aus dem Im- 
pulsintegralsatz. Gehen wir mit Gl. (1.7) in Gl. (1.5) ein, so folgt 


ar un(147 rn) RT AIR TE (1.9). 


In dem erhaltenen Formelsystem treten die drei Beiwerte B, U und o auf. Zu ihrer Bestimmung 
haben wir zunächst die in Gl. (1.7) gegebene Definitionsgleichung für o sowie Gl. (1.9). Eine 
dritte Gleichung erhalten wir durch das Ablesen der Breite des berechneten Profils Gl. (1.7): 


r 
(0 N) Breite = De BE ee N A, Rn (1.10) 
Die drei Beiwerte berechnen sich dann zu 
4 a er I ee 
ee 1.11% 
An 
1 (u/u)"® 1 
——— RR 2 ee ee 13 12)8 
een 5 
ug T, 
Ei — R 
1 43 Uy 1 
„B=,(6m0oB) — 5 (1.13). 


ae wur A1/SE 
(u,/uo)?! 1 + Sr 
T, 
*) Dabei ist zuzulassen, daß die eine oder andere Potenz von x/r, gegenüber den anderen vernachlässig- 
bar ist; im vorliegenden Beispiel ist das — siehe unten — die Potenz (</ro) )n-tm—n, 
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Dieses Formelsystem, das in seiner Aussage ursächlich an die Prandtische Formulierung des 
turbulenten Austausches gebunden ist, legt nun zusammen mit Gl. (1.7) und Gl. (1.8) das Strö- 
mungsfeld der Asymptote in allen Einzelheiten fest. Dabei nehmen wir an, daß unter sonst 
gleichen Versuchsbedingungen (vgl. Abschn. IV) x, eine Konstante ist, die unabhängig von den 
Strömungsparametern (u, — u,)/u, und d,/T, ist. 

Die Mischbreite des Geschwindigkeitsfeldes sowie der Abfall der Strahlgeschwindigkeit 
längs der Achse hängen demnach sehr wesentlich davon ab, ob es sich um einen Heißluftstrahl 
handelt oder nicht. ; 

Man kann aus den Gleichungen u.a. ablesen, daß für heißere Düsenstrahlen einerseits 
der Abfall der Geschwindigkeit längs der Achse schneller erfolgt und andererseits die Mischbreiten 
kleiner sind. Der in Abschn. V durchgeführte Vergleich mit dem Experiment wird eine Bestäti- 
gung dieser Aussagen ergeben. 


Temperaturfeld: 
Der Ähnlichkeitsansatz lautet 
du) r/Tg 
== m N= ar 
5 d, x) l (z/r,)"® (1.14) 
wir setzen dabei = Tan)" 
0 
Der Energieintegralsatz 
r = 
eudrd=UÜ zn (1.15) 
0 . 
liefert in der Asymptote und mit dem Ansatz Gl. (1.14) 
ae ik ee 
u E d= . — rer i 
" UT rwm1+ N) 
also 
2 
ne An) rs RER ARE 
Gleichung (IX) ergibt dann für x 
dx wW—u,[ d? d 
28) Be N nl Bl) 2,04 
IKT dn GR TR. dı? + = 
mit den Randbedingungen 
0) -1, x(©) = 0 
und der Lösung 
1 1 
| -rem 
a num: 


Bezeichnen wir mit b,(x) die Mischbreite des Temperaturfeldes, so erhalt 2 5 

unmittelbar für das Verhältnis der Mischbreiten i et 
b, 
b, 


ee. 0: 


Eine wichtige Relation folgt noch aus dem Energieintegr wi EG 1 
| :h au: ‚nergieintegralsatz. ir 
Een oo Korn g gralsatz. Gehen wir mit Gl. (1.18) in 


2 m d 
02 Er | f n) a en 
Der Vergleich von Gl. (1.9) und Gl. (1.20) ergibt dann das Gesetz 
er) 
en, 


Man kann hier schon von einem Gesetz sprechen; 

i ıesetz sprechen; denn Gl. (1.21) folgt aus den hypothesenfrei 
ne des Impulses und der Energie mittels des Ännlichkeits ne A ns Rela- 
un X g'®, Beide letzteren Annahmen sind durch das Experiment mehrfach bestätigt wor- 

en; vgl. hierzu auch Abschn. V. Die Glockenkurve der Funktion p nach Gl. (1.7) schließlicl 
ermöglicht auf jeden Fall eine gute Approximation. von | r 


PO 
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Wir berechnen noch die Fördermenge des Strahls 


[0,0] 
a NR, BR UT (1.22) 
bzw. 
F eU—uT Tr 
Ze a dee 
F% ; WU Tg To 


wo F, = % (U, — u,)rr, die Treibmenge®) des Strahls ist. In der Asymptote ergibt die Inte- 


5 } d et 
gration mit S = ee = und bei Einsetzen von Gl. (1.7) gemäß Gl. (1.1) 
0 (1) 1 

Fe 2) 
a. 

bzw. nach Gl. (1.9) 
F Uy 
N. (1.23) 


Die Fördermenge des Düsenstrahls nähert sich also asymptotisch einem konstantem Wert. 


2. Ruhende Außenluft (u, = 0) 
Hier liegt für das Geschwindigkeitsfeld konstanter Dichte die Lösung von H. Schlich- 
ting!®) vor. 
Mit dem Ähnlichkeitsansatz 
0 FREIE, v 1 


= 0) 
Be alt” 777 


werden die Profile durch die Gleichung 


d d dy 1d 
Ute) Hr meh: En 


mit den Randbedingungen (0) = 1, p(®) = 0 sowie durch die Gleichung 


N 


d d 
tn =Unle+ng}) 


mit der Randbedingung y(0) = 0 bestimmt, welche die Lösungen 


1 : 1 U 1— (on)? 
=. 0 2 — —— ; — —— 0... (2.2 
ee ne EEE 
haben. 
Der Impulsintegralsatz 
oo 72 
lewrar-au. ee NDR) 


0 


liefert uns wieder die Kennzeichnung des Heißluftstrahls mittels der asymptotischen Relation 


U2 | - 5 
eh ne One EN Picyn BE REEN ERTL, E00 ra oe 2.4 . 
Wir können dann wieder bei Ablesung der Breite o B des berechneten Profils @ die Beiwerte 
berechnen: 


0% 


ıl 

ee a a Dr a a REN RS 
5 8x,0B ) 
ERS ee 2.6 
ee 2) BE ri 

Ilm: 

T, 
4 B= 5 @4,0 B} ir . N % & . 5 . « s hr . . . (2.7). 


9) Bezeichnung nach W. Wuest, V.D.I. 92 (1950). S. 1000. 
16, H.Schlichting, a. a. O, 
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i i i i Luft, vgl. Gl. (1.13) — die 
Es fällt auf, daß hier — im Gegensatz zur Ausbreitung in bewegter 
u oischen Mischbreiten nach Gl. (2.7) keine Abhängigkeit vom Temperaturparameter 
9,/T, aufweisen. Für das Temperaturfeld erhält man 


ER! BE : _ 
d =) nr 
Das Profil wird durch die Gleichung 
a2) % -Erhre dy 
Uplern u u! (+ 
mit den Randbedingungen x(0) = 1, x(©) = 0 bestimmt und lautet 
1 
ee ne TELE 0 
ar OTTHEnPe 
Der Energieintegralsatz 
2 
|eudrar- win. ee 
Ö 
liefert die asymptotische Relation 
UT. a 
a _ — (1 N RAN 
"mtr | 
Vergleicht man die beiden Relationen Gl. (2.4) und Gl. (2.11), so folgt 
LIE RT = 
U 3 


(Anmerkung: Sollte sich — wie derzeit nur der Vergleich mit dem Experiment zeigen kann - 
hier die Funktion 9 = e”(°”* als eine bessere Approximation erweisen, so würde durch Eingang 
mit dieser in die Integralsätze 


folgen.) 
Für die Fördermenge 


F= an [BUurde .. 2 VRR 
0 


erhält man mit der Treibmenge F, = 0, uy x r} die asymptotische Relation 


bzw. nach Gl. (2.5) und Gl. (2.6) 


' = j ) 1/2 E 
2 en v3 (8 x [03 B) 1 - H- do 5 . - » . - . . . . . . (2.1 ». 
F) T, To 


Bemerkung 


Der in Erscheinung tretende fundamentale Unterschied zwischen den Fällen bewegter und 
ruhender Außenluft, der einen Widerspruch gegen die zu fordernde stetige Änderung physi- 
kalischer Zustände anzuzeigen scheint, ist nur durch die Beschränkung auf die Asymptote be- 
dingt. Der Widerspruch löst sich auf, wenn wir den Vorgang der turbulenten Ausbreitung eines 
Düsenstrahls als Ganzes betrachten: 

Bei ruhender Außenluft wird der asymptotische Zustand bekanntlich bereits in verhältnis- 
mäßig kleinem Abstand von der Düsenmündung realisiert; vgl. Abschn. V. Lassen wir nun eine 
kleine Außengeschwindigkeit u, zu, so handelt es sich um eine kleine Störung jenes Zustandes, 
die erst in sehr großem Abstand — wo u, von gleicher Größenordnung wie die Strahlgeschwin- 
digkeit ist — wirksam wird. Bis zu größeren Abständen von der Mündung haben wir also hier 
im wesentlichen den Status wie bei ruhender Außenluft zu erwarten; der asymptotische Zustand 
wird erst in sehr großem Abstand anzutreffen sein und verschwindet in der Grenze u, —0. Mit 
größer werdender Außengeschwindigkeit u, ist die Realisierung des asymptotischen Zustandes 
dann in kleineren Abständen zu erwarten; vgl. Abschn. V. 


ar be ” 
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B. Ebene Strömung 
Wir können uns im Folgenden kurz fassen. 


Die asymptotischen Bewegungsgleichungen in der Hauptströmungsrichtung reduzieren 
sich, wie in Abschn. II dargelegt, auf die für ein Feld konstanter Dichte: 


Impuls 
ou ou Au 
or az @, ro or och Sy nr Va oe (XD), 
Energie 
0» od 0» 
u 2x +D y E £e(X) ap N ee Bun (XI), 
Masse 
cu 
Zr + 24 NET EN BRETT NIE (XII). 


x bedeutet dabei die Koordinate in Richtung der Strahlachse, yin der dazu senkrechten Richtung. 


3. Bewegte Außenluft (u, # 0) 
Wir machen wieder den Ähnlichkeitsansatz 


el n.g). so a ae 


1 — u, 1 — U, 
0) x \P 
N 4) von; Dt). 
dabei setzen wir 
u,(x) Ad u > SF b, b* Er 
I — U lg j De ) 


Der Vergleich der Exponenten von x/d an Hand der Bewegungsgleichungen (XI) und (XIII) 
sowie des Impulsintegralsatzes ergibt 


n= 5; m-—5; De Bar 
Für das symmetrische Profil p(n) folgt aus (XI) die Gleichung 


1 


= h ‚U, — u, d’p 
ale tn) mBt ( 


mit den Randbedingungen 
0) =1, pa)=0; 


die Lösung lautet 
I 


m en: i I = —  ——————.. 1.2.2... 33 
nt a 4x, BU (w — u,)/u, oo 
Für das Profil y(n) erhält man die Gleichung 
dy a do 
arten) 
mit der Randbedingung y (0) = 0 und der Lösung 
U 
a RE ELER N ee 
Der Impulsintegralsatz 
Beni eduiin U)dı.... 20 00.) 
0 
ergibt dann die asymptotische Relation 
_ IR do 3.6 
o-Fuilıı?) FR . (3.6). 
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Die Beiwerte ergeben sich zu 


a yr je u,/ug ai Oo \W2 Be 
en y P x oB E = ie 1 + ) > Ra ( )s 
Ug 

en ee e 
U yr (5 % oB Be = A FE ae Be: ( yM 

Ug I, 

En 

N - 2x, 0. B\?? Ip: ne 1 e$ 
B= ja| | Ta rau en 

T, 


und zeigen die Abhängigkeit des asymptotischen Strömungsfeldes von dem den Heißluftstrahl 
charakterisierenden Parameter %/T;. 


Für das Temperaturfeld erhält man mittels des Ansatzes 


DO Bu 
= Mn) 
aus dem Energieintegralsatz (3.13) 
Ra 
» =D a 0 en = N EEE 
Gl. (XII) ergibt für das symmetrische Profil y die Gleichung 
1 dx ee: 
let E Ban u d? 
mit den Randbedingungen 4(0) = 1, xy(©) = 0 und der Lösung 
= — (on) 

Ve N si I TR Si ae 
aus der für das Verhältnis der Mischbreiten des Temperatur- und Geschwindigkeitsfeldes folgt 
b 4. 

SuVE.. 2 A 
b, 
Der Energieintegralsatz 
eu dy— MRS 2 
0 
liefert dann asymptotisch 
u DI 
o- Euer! “(14 2) a 
Aus dem Vergleich von Gl. (3.6) und (3.14) folgt die Beziehung 
Tal 
nn VE - 13.18); 


der — wie der entsprechenden Formel (1.21) des rotationssymmetrischen Falles — wieder eine 
gewisse Allgemeingültigkeit zukommt. 


Für die Fördermenge 


Pi zn U UN. 
erhält ınan mit der Treibmenge F, — £, (u, — u) 2 d asymptotisch 
F U Ü 
F_@U( 9) 
F, 2 Tui 
bzw. nach Gl. (3.6) 
I Uy 


Fe a 
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4. Ruhende Außenluft (u, = 0) 
Hier liegt für das Geschwindigkeitsfeld konstanter Dichte die Lösung von H. Görtler!!) vor. 


URRANT Sl a! 
San Ben. Se 
Die symmetrischen Profile werden durch die Gleichungen 
dp 
= u.(2 Hug) tr = re 
Randbedingungen (0) = 1, g(®) = 0; 
und 
en a ah Randbedingung (0) = 0. 
dn dn!’ 
bestimmt. 
Die Lösungen lauten 
N 2 1 
+ —2 + ee er Bee ‘ 
= 0p) "(or) mit 20 = 2.0 oy u (on &of (on) 3 TIg(lon)) (4.2). 
Der Impulsintegralsatz 
old. 0, 1 dis se. Re (4.3) 
ergibt dann asymptotisch 
2 | d 
=— De ne 7; t 
ll, | (4.4) 
Es folgen die Formeln 
1 
NE RETTIEENENETÜENEEEHEE 4.5), 
ee 70 B (45) 
1 1 
= l\/- Re (4.6), 
” E (4x, 0 BJ? (1 + 8) T"? z 
uB- a ENT go 20 a Re a (4.7). 
Temperaturfeld: 
v T 
» =. (Jay? 4m) Se a Er a en (4.8) 
mit der Gleichung 
U 2, = Em BO 
= „(& ng +Y- ale 
und den Randbedingungen x(0) = 1, x(o) = 0. 
Nach H. Reichardt!2) haben wir die Lösung 
VE N u U CE) re er re N, 
Der Energieintegralsatz | 
fouddy = U, d, d a ar er A = (4.10) 
0 
ergibt asymptotisch 
f DR in es 
= ur(i a Cof | #2) (ae BER Kerr 
T, 
0 
und der Vergleich von Gl. (4.4) und (4.11) liefert die Beziehung 
oo rs B: SEE, 
= = j Cof +2) gas) Were NR), 


0 


11,7; @ärtler, a. a0). € 
12) H. Reichardt, Z. angew. Math. Mech. 24 (1944), S. 268. 
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i i j i Heißluft- 
Nach vorliegenden Messungen, vgl. auch Abschn. V, haben wir bei der Mischung von 
strahlen für das Übertragungsverhältnis den Zahlenwert E = 2. Aus Gl. (4.12) folgt dann 


re 
Udo 
(Anmerkung: Approximiert man in den Integralsätzen durch die Funktion @ = e”), 

so erhält man mitE=2 


(208 ea ee 


Te | 
_—— — x 0,866) , = 
j IL@Z ( } 
also fast denselben Wert.) 
Für die Fördermenge 
F= foudy..... un 


ergibt sich mit der Treibmenge F, = 0 uy2d die asymptotische Formel 


Es 17 a “ . .(4.15). 
5 5402 +2 Auen : 


IV. Der Prandtlsche Ansatz für den turbulenten Austausch bei freier Turbulenz 


In voller Allgemeinheit kann man im Fall konstanter Dichte (bei ebener oder rotations- 
symmetrischer Strömung) offenbar für die turbulente Schubspannung setzen 


2) - 
Terre ren 


dabei soll bei Betrachtung von Strömungen der freien Turbulenz x die Hauptströmungsrichtung, 
y die dazu senkrechte Richtung bezeichnen. Das Symbol e(z, y) bedeutet die sog. scheinbare 
kinematische Zähigkeit. 


1. Die erste wesentliche Aussage des Prandtlschen Ansatzes ist, daß man für Schnitte 
senkrecht zur x-Achse sich mit Mittelwerten des turbulenten Austausches begnügen kann: 


x, y) = al)? II am Me 


Dieser Ansatz hat seine Tragfähigkeit an Hand der bekannten guten Übereinstimmung der so 
berechneten theoretischen Profile der Geschwindigkeitsverteilung mit dem Experiment erwiesen. 


2. Hat man den Mittelwert e(x) akzeptiert, so verbleiben als bestimmende physikalische 
Faktoren offenbar nur die Mischbreite b,(x) und die Geschwindigkeitsspanne Uyar — Umin- Aus 
dimensionsanalytischen Gründen folgt dann sofort 


et) = able) lu mt en 
also der bekannte Prandtlsche'?) Ansatz mit dem dimensionslosen empirischen Koeffizienten x,. 


3. Da es sich bei obigem Ansatz um eine Mittelwertbildung handelt, ist von vornherein 
nicht zu erwarten, daß der Koeffizient x, für die Gesamtheit eines Strömungsvorganges, also z. B. 
für einen Düsenstrahl von der Mündung bis zur Asymptote, gleich einer Konstanten wäre. Doch 
wird man, wenn andererseits dem Ansatz eine Bedeutung zugemessen werden soll, erwarten 
dürfen, daß dieser Umstand nur geringe Schwankungen des Zahlenwertes von x, zur Folge hat. 

Weiter ist zu bemerken, daß nach (y) Schwankungen von x, dadurch bedingt werden 
dürften, daß die Festlegung eines Maßes für die Mischbreite der längs des Strömungsvorganges 
variierenden Profile, die mit gegen -- © wachsender Querkoordinate y nach Null abklingen, 
einer gewissen Willkür unterliegt. Aber auch hier dürfte zu erwarten sein, daß die Schwan- 
kungen sich in engen Grenzen halten. 

Insgesamt gesehen, sollte x, den Charakter einer Modellkonstanten haben. 


Dagegen ist mcht zu erwarten, daß x, den Charakter einer universellen Konstanten hat. 
Die Verhältnisse dürften hier ähnlich liegen wie bei dem Koeffizienten x des logarithmischen 
Geschwindigkeitsgesetzes turbulenter Grenzschichten, für den die ermittelten Zahlenwerte eine 


merkbare Streuung zeigen, wenn sie auch in der Mehrzahl den Wert x x 0,4 ergeben. Doch läßt 
sich zur Zeit wenig dazu sagen. 


13) L. Prandtl, Z. angew. Math. Mech. 22 (1942), 8. 241. 
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V. Vergleich mit Messungen 

Messungen von Heißluftstrahlen sind bisher anscheinend nur für runde Düsenstrahlen 
durchgeführt worden, auf die sich der Vergleich beschränkt!#). 

Die Untersuchung der turbulenten Vermischung ebener Heißluftstrahlen!5) — also des 
Modells, das durch das Mischungsfeld in unmittelbarer Nähe der Düsenkante 
realisiert wird — hatte bei Ablesung der Breite o B zwischen den Grenzen 0,95 und 0,05 den 
Zahlenwert x, = 0,0082 ergeben. Das genannte Maß der Breite o B behalten wir im Folgenden 
bei. Der hier gezogene Vergleich mit Messungen der turbulenten Ausbreitung rotationssymme- 
trischer Heißluftstrahlen in größeren axialen Abständen bei sowohl ruhender wie bewegter Außen- 
luft ergibt dann das bemerkenswerte Resultat, daß in der Asymptote nur ein geringfügig von der 
genannten Zahl 0,0082 abweichender Wert des Koeffizienten x, anzunehmen ist, was nach den. 
Erörterungen von Abschn. IV den Prandilschen Ansatz als sachgemäß erscheinen läßt. 

Wir wollen noch hinsichtlich des Vergleiches mit dem Experiment bemerken, daß allerdings 
gerade die für die asymptotische Theorie interessanten kleinen Werte der Geschwindigkeits- und 
Temperaturdifferenz des Strahls gegenüber dem umgebenden Medium in größeren axialen Ab- 
ständen nur schwer meßbar und daher oft mit Unsicherheit behaftet sind. 


A. Ruhende Außenluft (m, = 0) 


Hier liegen die Messungen von S. Corrsin und M.S. Uberoi) an Heißluftstrahlen mit 
der Übertemperatur 9, = 15°, 170° und 300° C gegenüber der Außenluft bzw. 


do 
T 0,05 0,57 1,0 
vor. 
Der Vergleich ergibt hier 
a WOOSDR EEE ME a re) 
17 
du) 2. u) 
PDS Ra 
u (X) 
ed 


09-08 -07 -06 -05 -06 03 -@2 07 0 
| 09 
% |08 

07 


06 


x =0 Messung 5. Corrsin and M.S. -02 
2 Uberoi,Naca Rep. 998 (1950) 


05 


Theorie 
04 


03 


02 


ıg al) 
% 
- -70 
01 02 03 0 05 06 07 08 0970 
Ua (IT) 
Vo 
Bild 2 


14) Für ebene Düsenstrahlen liegt für den Fall konstanter Dichte die Untersuchung von E. Förthmann 
[Ing.-Arch. 5 (1934), S. 42] vor, der die turbulente Ausbreitung in ruhender Luft für größere axiale Abstände 
gemessen hat. Der Vergleich mit der Theorie nach Abschn. III, B, 4 würde hier einen wesentlich größeren 
Zahlenwert für x, als bei rotationssymmetrischen Strahlen ergeben. Wir gehen an dieser Stelle nicht weiter 
darauf ein. 

15) W.Szablewski, a.a. O., S. 2. 

16) S. Corrsin and M. S. Uberoi, Naca Rep. 998 (1950). 
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Diese aus Gl. (2.12) mit E — 2 (vgl. 4.) sich ergebende Relation, die in der dort geschilderten 
Umschränkung als hypothesenfrei anzusehen ist, wird nach Bild 2 durch das Experiment be- 
stätigt. Im Fall 9,/T, = 0,05 liegen für die kleineren Werte u,(x)/u, die zugehörigen experimen- 
tellen Werte 9,(2)/9, zu niedrig; vgl. hierzu auch Bild 5. Der Fall ®,/T; = 0,57 wurde nicht 
berücksichtigt, da die experimentelle Temperaturkurve infolge Überschneidens anderer experi- 
menteller Temperaturkurven als unsicher anzusehen ist. 


2. Temperatureinfluß auf den Verlauf der Achsenfunktionen 
Nach Formel (2.6) ergibt sich die — von x, freie — Relation 


uue) ee 

Up Klar Pe 1+ (W/Ty)zı 

u4() 1+(@/Ty)ı 
Up KlTa)yı 


und die analoge Formel für 9,(@)/do- 


v4 (X) 
25 "=005 


-0 -09 -08 -07 -06° -05 -06 -03 -02 -U 0 


elg 


ı0 
Baer 09 
vo || ge 
Tr, | For 
07 
06 . = = 057 Messung 5.Corrsin and M.S. -02 
ö . Uberoi, Naca Rep. 998 (1950) 
05 Z10 -03 
Theorie 
04 -04 
03 Ss 
-06 
02 -07 
-08 
09 I a) 
g Up % 
-10 1) 
01 02 03 0% 05 06 07080910 


ualX) 


u ae 
Ya =005 
ı 


Bild 3 


Aus Bild 3 kann man ablesen, inwieweit die Messungen diesen Effekt zeigen. Für 9/7, =1,0 
ist die Übereinstimmung recht befriedigend, weniger für d/Tı = 0,57. Auf die analoge Dar- 
stellung für die Achsenfunktionen der Temperatur | 

se B wurde aus den unter 1. genannten Grü 
verzichtet. i wa 


3. Verlauf der Achsenfunklionen der Geschwindigkeit und der Temperatur 
Nach Gl. (2.6) erhalten wir mit (*) und dem aus Gl. (2.2) folgenden Maß (co Be — 
0,05 x 


in der Formel u,@)/u, = U(z/r,)! für 9/T, = 0,05 bzw. 0,57 bzw. 1,0 U ) 

’ Et = 0,0% RN .1,0 U = 14,62 bzw. = 11,96 
bzw. = 10,59. Aus Bild 4 ist ersichtlich, daß die Theorie befriedigend wiedergegeben wird, wenn 
auch die theoretische Kurve für 9,/T, = 0,05 etwas zu hoch liegt. ; 


Für 9,@)/d, = T(e/r,) ! erhält man nach Gl. (2.12)mit E = 2 für 
e _erhä an nac Ale —=2 für 9/T, = 0,05 bzw. 1,0 
die Werte T = 9,74 bzw. = 7,06. Der Fall 9,/T; = 0,57 wurde aus ne 1. erwähnten 
Grunde außer acht gelassen. Der in Bild 5 gezogene Vergleich zeigt eine weniger gute Überein- 


stimmung, insbes. für den Fall 9,/T, = 0,05. Allerdings ist zu beachten, d i i 
* » a i = \ ® i i z aß 
nach 1. die kleinen 9,(z)/d,-Werte zu niedrig liegen. ö ne 
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® = = 0 05 
ne h 057 Messung 5. Corrsin and M.S. 
J En Uberoi,Naca Rep, 398 (1950) 


. 


Theorie 


01 - -/0 
STEIN TO) 20 30 +44 50 60 70 8090100 


Bild 4 


Messung 5. Corrsin and M.S.Uberoi, 
Noca Rep. 398 (1950) 


Theorie 


01 
ee ET IT, 20 30 44 50 60 70 80 90100 


Bild 5 


4. Profile der Geschwindigkeit und Temperatur über den Mischbreiten 

Meßwerte sind hier leider nur für 9,/T, X 0 von S. Corrsin"”) publiziert worden. 

Nach Gl. (2.2) haben wir für das Geschwindigkeitsprofil 

1 EI, 
PTnronm am’ 

wobei nach GI. (2.5) mit (o B).« — 1,7 und mit (*) o = 8,65 ist. Der Vergleich in Bild 6 ergibt 
ausreichende Übereinstimmung sowohl hinsichtlich der Profilform als auch hinsichtlich der Misch- 
breite. 


2) S. Corrsin, NACA - Wartime Rep. W - 94 (1946). 


“ 
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Für die Temperaturverteilung (2.9) 


ip" 
ist aus Bild 6 für die experimentellen Werte 
E=2 


ablesbar. Das ergibt y = p!?. Es zeigen sich die bekannten stärkeren Abweichungen im Rand- 
gebiet, wenn auch die experimentellen Punkte hier ziemlich stark streuen. 


Messung S.Corrsin, 
Noca Wortime Report W- 94-1946) 


2,0 °® 
60 » Geschwindigkeit 
80 ° 
o 
° Temperatur 


Theorie 


02 ar 0 ) 02 03 0% 
Bild 6 


B. Bewegte Außenluft (u, #0) 


Wie wir in Abschn. III, 2 (Bemerkung) dargelegt haben, erscheint unsere Theorie der 
Asymptote nur bei größeren Außengeschwindigkeiten im üblichen Meßbereich (der sich bis 
ungefähr (u,(x) — u,)/(ug— u,) bzw. 9,(x)/d, X 0,2 erstreckt) realisierbar. Wir werden das 
weiter unten auch experimentell belegen. 

Zum Vergleich mit der Theorie ziehen wir zunächst eine Messung von O. Pabst!®) mit den 
Daten 

W— U, Ö 


05 = 1,0, 
Ug T, 
sowie eine Messung von F. Landis und A. H. Shapiro!P) mit 
U, — 
a 08: vo _03 
U, T, 
heran. 
Der Vergleich ergibt dann hier 
1.3 U,DOBR 0.000 as na A 


dx) " 1u@—u, 
u 2 nen 


Diescs Gesetz, das sich aus Gl. (1.21) mit E = 2 ergibt, wird durch die genannten Experimente 
auch nicht annähernd wiedergegeben. Ja, das Experiment zeigt sogar Überschneidungen der 
Achsenfunktionen der Geschwindigkeit und Temperatur. Da andererseits sich aus dem Experi- 
ment recht deutlich E = 2 ergibt (vgl. 4), so muß man nach den Darlegungen im Anschluß an 
Gl. (1.21) annehmen, daß diese Diskrepanz zu Lasten des Experiments geht. 


Wir führen im Folgenden den Vergleich nur für die Achsenfunktion der Geschwindigkeit 
weiter. 


1) Q,Pabst, U. u. M. 8007 (1944). 
») F. Landis and A. H. Shapiro, Heat Transfer and Fluid Mechanics Institute 1951. 
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2. Temperatureinfluß auf den Verlauf der Achsenfunktion der Geschwindigkeit 
Nach GI. (1.12) ergibt sich die — von x, freie — Relation 


ul) — U. U Er 
U—u, karmdr _ / + BolTdır 

ul) — u, A+/T)ı 
U — U, BulTı)yı 


In Anwendung auf die genannten Messungen folgt 


ul) — U, ud 
W—U, ym=10 _ ee 
u > 
U — U, [8/7 = 0,3 
wenn wir den Einfluß der Differenz zwischen den Geschwindigkeiten en 0,53 und 0,50 
vernachlässigen. Us 


In Bild 7 haben wir einige zugehörige Punkte aufgetragen; dabei wurden die den Meß- 
punkten von O.Pabst entsprechenden Werte der amerikanischen Messung einem Isotachen- 
diagramm entnommen. 


05 
0. Pabst, UV.u.M. 8007 (1944).- 
Messung «\F. Landis and A.H.Shapiro, 
Heat Transfer and Fivid Mech. 
Institute 1951 
Theorie 
—A=- 03 
0 05 10 


Bild 7 


3. Verlauf der Achsenfunktion der Geschwindigkeit 
Wir haben hier nach Abschn. III, 1 die Formel 


(ua) — u)) — U(x/ry 23 . 


(u, — U) 
0,95 


U berechnet sich nach Gl. (1.12) mit (**) und mit dem aus Gl. (1.7) folgenden Maß (oB) _ =1,5. 


0,05 


a a Een m ER 053.0 10 
Ug T, U, T, 
U = 6,55 bzw. 5,345. 
Bild 8 zeigt, daß die theoretischen Funktionen in der Größenordnung richtig liegen. Aller- 
dings könnte man aus den Meßpunkten unschwer auch eine andere Neigung ablesen, als die 
Theorie angibt. Da jedoch die Profile — vgl. Bild 92%) — das Ähnlichkeitsgesetz mit dem Misch- 


20) Auf die Wiedergabe der aus dem Isotachen- und Isothermendiagramm von F. Landis und A.H. 
Shapiro ablesbaren Profile haben wir verzichtet. 


D 
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1/3 r 
breitenverlauf A bzw. dr m (=) annähernd befolgen, folgt aus dem Impuls- bzw. Energieintegral- 


To () 0 3 ar 
satz zwangsläufig er bzw. 94) m & i 


U, f) To 
Unl&)-U, ig = ice 
Kozyı 05 10 15 
10 0 
09 


04 e Messung FE Landis and A.M.Shapiro, 
z Heat Transfer and Fluid Mechan. Institute 1951 


« Messung 0.Pabst, U.u.M. 8007 (1944) 


Theorie 


en NT DT Er Bu 2) 0 0 01” 


os 


Bild 8 


4. Profile der Geschwindigkeit und Temperatur über den Mischbreiten 
Nach Gl. (1.7) haben wir das Geschwindigkeitsprofil 


2 T/Tg 
= — (an) u 
N 
co berechnet sich nach Gl. (1.11) mit (**) und dem Maß (co B) a — 1,50. Man erhält für 
Ti, nr ee o= 2,4. 
Ug T, 


A »- 
SSPE 1 


Uo - 
Messung 0. Pobst, U.u.M. 8007 (1944) 
; 
a 
A SE 
= = . 
o (fo n, = 24 
a s 32 2 ee 
Temperatur Geschwindigkeit 
o 40 m 
v 48 v 
Theorie 
asymptotisches { oSsymptotisches 
Profil a h Profil 
5 
& 
wo 
 B:; \ 
vo SNME e r/rg 
v bg 1/3 
°‘. an =) 
eree je (EEE RTBeR l 4 ae . 
1 a 


Bild 9 
Mit E=2 wird das Temperaturprofil nach Gl. (1.18) durch 


A 
beschrieben. 


- 
3 


# 
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Bild 9 erhält die Messung von O. Pabst. Sowohl in bezug auf die Profilform als auch die 
Mischbreite ist die Übereinstimmung als ausreichend anzusehen?!). Insbesondere kann man aus 
der Messung von O.Pabst E=2 ablesen. Daß die asymptotischen Profile als „äußere Be- 
grenzung‘ der Schar der experimentellen Profile erscheinen, ist als sachgemäß anzusehen, wie 
wir noch auseinandersetzen werden. 


Im Kernbereich einer Düsenströmung, der durch das Vorhandensein des Strahlkerns defi- 


niert ist, haben wir bekanntlich?) in sehr guter Näherung - bzw. 2 m X, infolgedessen 
Be 


0 0 0 
müssen auf dem Weg zur Asymptote die Mischbreiten den Übergang zur asymptotischen Formel 


b b a : 3 
= bzw. n m a vollziehen, worauf bereits O. Pabst??) hingewiesen hat. Dieser Übergang 
{) {) 0 

muß sich bei einer Auftragung der gemessenen Profile über der Ähnlichkeitskoordinate 
Ye (r/ro 
 (@frg)t® 

Hüllkurve‘ erscheint. Diese Fächerung wird, entsprechend unserer Bemerkung in Abschn. III, 2, 

um so markanter in Erscheinung treten, je kleiner die relative Außengeschwindigkeit u,/u, ist.) 


als Fächerung der Profile darstellen, wobei das asymptotische Profil als „äußere 


Ug-Ly [2 .% 5 
277 = 0% ; 7, I. 
Messung 0. Pabst, U.u.M. 8007 (1944) 
o er = 24 ® 
I) 
SLR 
Temperatur Geschwindigkeit 
o 40 0 
v 48 v 
Theorie 
asymptotisches a osymptofisches 


Profil 


Bild 10 


Bild 10, in das wir für (u, — u,)/ug = 0,74; 9,/T, = 1,0 eine Messung von O. Pabst?°®) und 
die theoretischen Profile mit o = 1,35 eingetragen haben, zeigt diese Fächerung deutlich. Auch 
dieser Messung kann man übrigens das Verhältnis E = 2 entnehmen. 


21) Nach (1.11) zeigt die Theorie als Temperatureffekt eine Verkleinerung des asymptotischen Aus- 
breitungswinkels mit wachsendem Parameter 9/7, an: Für (u, — u,)/üg = 0,5; %/Tı = 0,3 (F. Landis und 
A.H.Shapiro): o = 2,06; für (u — u)/w = 0,53; 9%/T, = 1,0 (O. Pabst): o = 2,24. Dieser Effekt findet 
sich experimentell anscheinend bestätigt. Dem Isotachendiagramm der zitierten Messung von F. Landis und 
A.H.Shapiro kann man folgende Mischbreiten, gemessen in n, bei &/r, = 48 entnehmen: 9 = 0,30; 0,25; 
0,15; 0,10 mit 7 x 0,50; 0,54; 0,64; 0,70. Diese Werte vergleiche man mit dem Verlauf der entsprechenden 
experimentellen Kurve in Abb. 9. 

22) Vgl. z.B. W. Szablewski, Ing.-Arch. 26 (1958), S. 358. 

28) O. Pabst, a.a. O0. 

24) Für kleine relative Außengeschwindigkeiten u,/uy # 0 z. B. decken sich bis zu großen x/r, hin die 
Profile bei Auftragung über der Koordinate y = (r/r,)/(z/r,). Trägt man nun die Profile über der Koordinate 
n = (r/r,)/(z/r)1/3= y(z/r,)2/3 auf, so erhält man demgemäß eine starke Fächerung. 


Manuskripteingang: 23. 4. 1958 
Anschrift: Dr. W. Szablewski, Berlin- Adlershof, Büchnerweg 43 
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Theorie und Experiment zum Reibungswiderstand it it 

der glatten rotierenden Scheibe bei turbulenter Strömung 
Von E.Broeckeer nel aaa 


Der Reibungswiderstand rotierender Scheiben in Gehäusen ist Iheoretisch und. EREEWER ntell 
die a als auch für die laminare und turbulente Grenzschicht zwischen der Keı 
und der rotierenden Scheibe bzw. der gegemüberliegenden Wand untersucht worden; entsprechen 
nisse sind für die frei rotierende Scheibe bei laminarer und turbulenter Grenzschicht bekannt. 2% er 
Theorie und Experiment sind für die verschiedenen Fälle in guter Überreinstimmung mit A PER 
. Falles der Scheibe im Gehäuse bei turbulenter Grenzschichtströmung, was bisher mit den Annahmen zur. 2 

Lösung der Grenzschicht-Impulsgleichungen begründet wurde. Die folgende Untersuchung zeigt, daß ds 
theoretische Ergebnis experimentell nicht nachweisbar ist und behandelt zwei Fälle unterschiedlicher Wand- 

begrenzung, die eine allmähliche bzw. schnelle Änderung des Drehmomentenbeiwertes vom unteren Grenzwert 
der Scheibe im engen Gehäuse zum oberen Grenzwert der freien Scheibe ergeben. . 


E The frictional resistance of rotating discs enclosed in a casing has been investigated theoretically and 
P experimentally both for luminar flow and for the laminar as well as for the turbulent boundary layer between 
5 


= 


core flow and rotating disc, or casing wall, respectively. Corresponding results are known for the freely 
. rotating disc both for laminar and turbulent boundary layer. 23 19 ae 
In all these cases the theory agrees quite well with the experiments, with the exception of the case of the 
enclosed disc with turbulent boundary layer. The explanation for this disagreement has always been sought 
in the assumptions made in solving the momentum equations for the boundary layer. The present investigation 
r shows that the theoretical result is not demonstrable by experiment; it treats two cases of different boundary 
conditions corresponding to slow resp. fast change of the torque coefficient with increasing casing dimensions, 
from the lower limit for a narrow casing to the upper limit for free rotation. 


La resistance & la friction de disques en rotation en boites a &i€ examinee de fagon theorique et experimentale 
pour le courant lamellaires ainsi que pour la couche limite lamellaire et turbulente entre le courant du 
noyau central et le disque en rotation, respectivement le mur oppose; des r&sultats correspondants sont connus 
pour le disque en rotation libre avec une couche limite lamellaire et turbulente. 

La theorie et l’experience se trouvent bien en accord pour les cas differents, a l’exception du cas du disque 
en boite avec courant turbulent de la couche limite, ce qui a &tE explique jusqu’ici par les hy 
employees pour la solution des equations de la force impulsive pour la couche limite. Dans l’examen 
suivant il est demontre que le resultat theorique ne peut pas ätre prouve par des experiences, et deux cas de 
murs limites differents sont traites quiont pour resultatun changement lent respectivement rapide du coejficient 
de mouvement de rotationde la valeur marginale inferieure du disque dans la boite Etroite jusqu’ä la valeur 
marginale superieure du disque sans boite. 


| 
e 
4 
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CoIpOTUBJIEHHE TPeHHA BPAmamImmxcH B KApTepax MMCKOB MIA NAMHHAPHOTO TeyeHun A 
TaK>Ke JUIH NAMHHAPHOTO HU TYPÖYJIeHTHOTO HOTPAHHYHOTO CIIOA MEeik]Ly OCHOBHEIM TeYEHHEM U 
BpammalmımuMcAH NUCKOM COOTBETCBEHHO IPOTHBONOJIOFKHOU CTEHKOH ÖBINO HCCHENOBAHO Kak 
TeOpeTuyecKH TaK H IKCHEPHMEHTAAIBHO. ÄHAJIOTHYHLIE Pe3yJIbTaTbI H3BECTHEI H IH CBOOOAHO 
BpammamImeroch ucKka IIpu JaMUHAPHOM U TYPÖYJIeHTHOM HOTPAHHYHOM CIIOW. 

‚ln pasaımyHBIXx CIIyyaeB MeKAy TeopHeü MH OIBITOM HMEeTCA XOPolIMe CoTIacHe € HC- 
KJIbIYeEHMEM TOTO CAYYAA, KoTNa AMCK HAXONUTCA B KaprTepe Upm Haumyum TYPÖYACHTHOTO 
TeyeHNA B IIOTPAHHYHOM CIIOIO. ITO NO CHX IIOP OÖAICHSIIOC HONYINCHHAMN, IPHHSATBEIMH IPH 
pemeHhnu yPaBHeHuUA UMIIyIbCaA IIPOTpaHnyHoro canon. B pa6oTe NOKA3bIBAeTcH, YTO Teo- 
peTruyeckuf pe3y.IBTaT HeBO3MO3KHO HOTBePAHTB IKCHEPHUMEHTANBHEIM IyTeM. Paccmarpn- 
BAIOTCH BA CIIyYası C PA3HBIMH BUNAMNM CTEHKU, NAIOINME IOCTENEHHO COOTBETCBEHHO CKOPOE 
U3MEHEHHE KO9ßhHLMEHTA KPYTUJIBHOTO MOMEHTA OT HNKHETO NPenena, Koria NuUCK IOTPysKeH 
B TeCHOM KapTepe, NO BEPXHETO IIPeNeta, COOTBETCBYIOIIETO CBOOOAHOMY NUCKY. 


1. Einleitung 


Außer für den laminaren Fall der frei rotierenden Scheibe [1], [2] und der Drehströmung 
in der Nähe einer festen Wand [5] ist neuerdings auch die allgemeinere Strömung zwischen zwei 
mit verschiedener Winkelgeschwindigkeit umlaufenden Scheiben als exakte Lösung der Navier- 
Stokesschen Gleichungen ermittelt worden [9]; wird insbesondere die eine Scheibe zur festen 
Wand, so ergibt sich bei verschwindendem s/a oder Re = a? w/v (Bild 1) eine lineare Verteilung 
der Tangentialgeschwindigkeit zwischen Scheibe und Wand, mit zunehmendem s/a oder Re 
beschränkt sich die Anderung der Tangentialgeschwindigkeit grenzschichtartig auf Scheiben- 
und Wandnähe unter Ausbildung einer Kernströmung, deren W inkelgeschwindigkeit im Grenz- 
fall zu 5 X 0,3 @ ermittelt wird [13]. Die reale Strömung führt jedoch für größere s/a oder Re 
zu in der negativen Radialrichtung sich ausbreitenden turbulenten Grenzschichten, worauf 
die Übersicht zur vorliegenden Theorie im folgenden Abschnitt Bezug nimmt. 

Außer der endlichen axialen Begrenzung wurde theoretisch auch der laminare Fall bei 
gleichzeitig radialer Begrenzung des Flüssigkeitsgebiets behandelt [13]; für den entsprechenden 
turbulenten Fall liegen experimentelle Ergebnisse vor [3], [4], [7], [11], [12], während die Theorie 
den Einfluß der zylindrischen Begrenzung vernachlässigt. 

Die vorliegenden Ergebnisse können für die Berechnung der Einzelverluste von Turbo- 
maschinen nicht ohne Einschränkung zugrunde gelegt werden. Dies gilt insbesondere für solche 
Ausführungen, bei denen der Drall der geförderten Flüssigkeit die Strömung an den Laufrad- 
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scheiben beeinflußt. Die Grundform dieser Anordnung erhält man nach Bild 1 aus der allseitig 
begrenzten Scheibe bei Entfernung der Zylinderwand in der nur seitlich begrenzten Scheibe mit 
im Außenraum ruhender Strömung, während bei der realen Laufradströmung im Außenraum 
eine Wirbelquellströmung vorhanden ist. 


Freie Scheibe Scheibe ım Gehäuse 
Wand Scheibe Wand Scheibe 
5 & 
2a 
w w w 
allsettigbegrenzt seitlich begrenzt tangential radıal 
Bild1. Scheibenanordnungen Bild 2. Geschwindigkeitsprofile, turbulente 


Grenzschicht 


2. Theoretische Ergebnisse für die freie Scheibe [1] und für die allseitig begrenzte 
Scheibe [3] bei turbulenter Grenzschicht 


Das v. Karmänsche Näherungsverfahren der Impulsgleichungen für die Grenzschicht an 
der rotierenden Scheibe liefert in 


ö ö 
Radialrichtung ar | Des ) e— j ne = Pro RBor. en elle), 
b 
i ö ö 
d 
Tangentialrichtung 1(r | Dru* “) —r2ß Er ( N DE «) re li), 
ö ö, 
sowie für die Grenzschicht an der festen Wand in 
9 ö 
d 
Radialrichtung TE ( | v2 «) — H Pdlserroe Dorn (2a), 
ö ö 
8 9 
d d 
Tangentialrichtung SE [" | vu «) — r2B Ar | v “) INT DOES (2b). 
ö ö 
Die angenommenen Grenzschichtprofile (Bild 2) 
z\1/7 z z\1/7 
„= 62) Be u A alı (2) la 
1/7 1/7 
v =, (5) a u 185) ER N a ne on (3b) 
führen mit dem Wandschubspannungsgesetz 
T y 1/4 
FT a 0,0225 ee ee ee ee sen (4 a), 
oe (u* + v*) Yu* HEY, 
1/4 
N DE I lie ne Ferien Map: (4b) 
e (u? + v?) ee) 


zu den Impulsgleichungen in r, ß, ® und ö, vü bzw. ®, v,. 
Die Lösung für die rotierende Scheibe lautet 


dead, De ae) au ER. (8) 
mit 

„0,0278 + 0,1944 Blo — 0,222 (Bo) e 

ee ee nn 
und 


y -( 02% = (1 + 9310 (1 — Bo) Wo). 2.2. (M. 


0,313 + 0,504 Bo 
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Damit wird, wenn ie 
’ M — 2 JT I T r2 dr 
{) 


das Drenmenent einer Scheibenseite und 


a. Ge te 


den Momentenbeiwert sowie 


die Reynolds-Zahl bedeuten, der reduzierte Momentenbeiwert (turbulente Strömung) 


2M 


für die rotierende Scheibe zu 
Cy=4nr = (0,0225)*/5 [x (0,313 + 0,504 B/o) ]/? (1 + &2)310 (1 — ß/o)%® (rotierende Scheibe) (11). 


Für $ = 0 folgt daraus die Lösung für die freie Scheibe mit & = a, = 0,162 zu 
Cir» = 0,0728: 7 (frei rotierende Scheibe) 7. 1.777 Fer (12). 


Der entsprechende Ansatz nach (5) für # und o, liefert wegen des positiven Vorzeichens 
von r,in (2a) eine imaginäre Konstante; Schultz-Grunow setzt deshalb ® = v, # und ermittelt 
durch Potenzreihenentwicklung 


Cy = 4n: 0,0225 - 0,367 (B/o)® (feste Wand)... ... (13). 


Durch Gleichsetzen von (11) und (13) folgt für das Verhältnis der Winkelgeschwindigkeiten von 
Kernströmung und Scheibe 


ßl® = 1/1,954 und damit Cy = 0,0311 (allseitig begr. Scheibe) .. . . . (14), 
wobei die Reibung an der zylindrischen Begrenzungswand vernachlässigt ist. 
u 


0526 


05 % 
fe) = falB/o) 


0 05 - .3/w 7 


Bild 3. Kennzeichnende Größen für die rotierende Scheibe bei 
turbulenter Grenzschicht 


3. Die kennzeichnenden Größen für die rotierende Scheibe bei turbulenter Grenzsehieht in 
Abhängigkeit vom Winkelgeschwindigkeitsverhältnis B/o, Bild 3 
Aus (11) und (12) folgt für das Verhältnis der Momentenbeiwerte 


Eule AuCarn R + 1,606 B/o) ir ( a Ka (1 — Bo) = f(Blw) . (15); 


X 1+03 


die Grenzschichtdicke läßt sich wegen (5) mit y = (v/oy15 f(B/w) in der Form 


ö [ro\l5 5 /a\88 0,072/x  \85 | 
(2) = 5) Eu \ + 1,606 ai) ni DS 


u rue 


. 
2 
= 
. 
3- 
s 
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schreiben, ferner wird die auf die Tangentialgeschwindigkeit bezogene radiale Wandschub- 
spannungsgeschwindigkeit 
vo 


le) = fBlo) ee Ber done, 


4. Allgemeine theoretische und experimentelle Ergebnisse 


Bild 4 zeigt außer bekannten Versuchsergebnissen einige Meßwerte [12] im turbulenten 
Bereich, auf welche noch näher eingegangen wird. Mit Ausnahme des berechneten Momenten- 
beiwertes für die Scheibe im Gehäuse bei turbulenter Grenzschichtströmung (*) sind die theore- 
tischen Ergebnisse für die Bereiche 


Laminarströmung zwischen Scheibe und Wand [3] Cy Ra —n.dis, 

laminare Grenzschicht an Scheibe und Wand [3] Cy Rei? ='1,334, 
dgl. an der freien Scheibe [1], [2] 133% 

turbulente Grenzschicht an Scheibe und Wand [3] cy Re! = 0,0311 (*), 
dgl. an der freien Scheibe [1] = 0,0728 


gut bestätigt. Der quantitative Unterschied zwischen Theorie und Experiment für den er- 
wähnten turbulenten Fall wird in [3]-und [8] mit den zugrunde gelegten Annahmen zur Aus- 
wertung der Grenzschicht-Impulsgleichungen erklärt. Dies befriedigt jedoch insofern nicht, 
als die übrigen mit vergleichbarer Näherungsrechnung erhaltenen Ergebnisse durch zuverlässige 
Messungen belegt sind; so liefert bei laminarer Grenzschicht die Näherung nach [3] für dıe freie 
Scheibe einen nur 2% größeren Momentenbeiwert als die exakte Lösung [2], während die Meß- 
werte für die Scheibe im Gehäuse zwischen der Näherung und der exakten Lösung [9], [13] 
liegen. 


Freie Scheibe: Theorie [1] [2]—-— Experiment [6] o o o 


allseitig begrenzte Scheibe: Theorie [3] - - —\ Experiment [3] 
005 (Zylinderwend .\s/a = 0,02 ® 
{ reibungsfrei) =00 *+ 


seitlich begrenzte Scheibe: Theorie (lamınar) [913] 


Experiment (turbulent) [12] 
ne > 
= : r | 
>= ae 
& 


[4 


s/@ = 0015 x 


001 


Bild 4. Drehmomentenbeiwert der rotierenden Scheibe 


5. Winkelgeschwindigkeitsverhältnis ß/» und Drehmomentenbeiwert CO, für die allseitig 
begrenzte und für die seitlich begrenzte Scheibe in Abhängigkeit vom Abstandsverhältnis s/a 


Nach (2) ist der Momentenbeiwert cy der rotierenden Scheibe eine Funktion von Re und 
Blo. Da cy > Re-U5 nach Bild 4 genügend bestätigt ist, kann sich die Betrachtung auf den 
reduzierten Momentenbeiwert Cy beschränken. Den Ausgang hierzu liefern experimentelle 
Feststellungen der Winkelgeschwindigkeit ß der Kernströmung: Für die allseitig begrenzte 
Scheibe liegt der Meßwert ß/® = 0,357 bei s/a = 0,30 vor [3], ein weiterer Wert ergibt sich 
nach [7, dort. Bild 20] zu $/o® = 0,43 bei s/a — 0,135, dessen Ermittlung Bild 5 zeigt, dabei 
ist (ra)? = E gesetzt, so daß 


1 
.__..@ 1/2 
Blo = —; (2 Apıle) 


leicht aus der Figur folgt. Entsprechende Druckmessungen an einer die rotierende Scheibe 
seitlich begrenzenden festen Scheibe [12] sind in Bild 6 dargestellt; dabei waren die Druckmeß- 
stellen linear in & angeordnet, um bereits im Experiment die Abweichungen gegenüber dß/dr =0 
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zu veranschaulichen. Wie später begründet wird, wurde der Druckverlauf durch AP, reduziert 
wobei mit (18) 
| AP, = o(a o]2)j2 
definiert ist, somit kann unmittelbar an der Abszisse abgelesen werden 
1/4 Bo)" 
ge 2 


Für größeren Scheibenabstand s beginnt nach Bild 6 von der Drehachse aus eine Abwei- 
chung von der starren Rotation der Kernströmung im Sinne der Potentialwirbelströmung. 


S=b=5mm 
Re= 2.075-10° 


05 
20 2633 mm 
! 
300mm Hg 200 700 [0] 
3 ae) 7 4p VARe— 05 [4 
Bild 5. Ermittlung von ß/w® bei allseitig begrenzter Bild 6. Ermittlung von ß/® bei seitlich begrenzter Scheibe nach 
Scheibe nach Messung [7] mit Wasser Messungen [12] in Luft 


Zur Umrechnung der gemessenen Momentenbeiwerte der Versuchsanordnung nach Bild 6 
(s,; =konstanter, s—= variabler Scheibenabstand) auf eine symmetrische Anordnung gleicher 
Scheibenabstände s mit beliebig dünner rotierender Scheibe wurden gesetzt: 


Mes = Gesamtmoment der Scheibe nach Bild 6, 


M’ = Moment der rechten Scheibenseite, 
M' = Moment der linken Scheibenseite, 
M„ = Moment am Scheibenumfang der Breite b, 


so daß 
M ges = M' + M’” + M,„ und für die symmetrische Anordnung M = M'’ folgt. 
Es können 


ad a 
M'=2rfurd, M"’=2r[uUrd 
0 0 


wegen 
u = ualrja)8l5 , u = a (rja)8l® 
in der reduzierten Form 
MAMI b 
REES 


geschrieben werden; zur vereinfachten Berücksichtigung von 
M,;=?2n@br,, 
sei 


1 


N 5 (Ta ie Tu) 


genommen, so daß 
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wird, wobei fürs=s, M’= M’” gilt. Mit (8) wird dann 

Bw 32m, 

En 

Die Meßwerte Jo und Cy = ty Rell5 sind in Bild 7 und 8 sowohl für die seitlich begrenzte 
Scheibe als auch für die allseitig begrenzte Scheibe eingetragen, wobei die Werte cy nach [7] 


Cy 


: Ir i ZT, 7 
der dortigen empirischen Gl. (9) in cy = 5 R,eRe— En Rey; umgerechnet wurden. 
Rotierende S, k 
5 ch 
Bo eisen bean 
De Te Bi (20) 
Ala rn Ri 
N XAu— 
Ns 
Se 
Sn, Messergebnisse: [7] + 
Fed 2 [3] x 
RN [12] 
en 
17) 01 02 —— 5/8 03 


Bild 7. Winkelgeschwindigkeitsverhältnis ß/® bei allseitig begrenzter bzw. seitlich begrenzter 
rotierender Scheibe in Abhängigkeit von s/a 


Man erkennt für den ersteren Fall ein ausgeprägtes Minimum von Cy und eine schnellere 
Zunahme an den Wert der freien Scheibe als für den letzteren Fall, wie auch aus dem unter- 
schiedlichen Verlauf von ß/® hervorgeht. Insbesondere liefert das Experiment 


ee) I 
as I a Er el ORTES 
sja 0 ® 2 ( ) 
gegenüber dem theoretischen Ergebnis 1/1,954 nach [3]. Die Gültigkeit des theoretischen Wertes 
wird trotzdem nicht widerlegt, weil die Theorie von s/a unabhängige Grenzschichten annimmt 
und wegen der verschiedenen Vorzeichen von r, in den Impulsgleichungen (1a) und (2a) den 


Rotierende Scheibe 
seitlich begrenzt 


0.0728 (s/2 =) 


Gin. (I1u.23) 


Messung [127 


alıseitig begrenz!_ _ ____— 


004 


Re = 2,075: 10° 
002 


[) 07 02 —n je 03 


Bild 8. Reduzierter Momentenbeiwert der allseitig begrenzten bzw. seitlich begrenzten rotierenden 
Scheibe in Abhängigkeit von s/a 


Wert ß/® = 1/2 ausschließen muß, während beim Experiment bei abnehmendem s/a die Kern- 
strömung allmählich verschwindet und die Strömung dem laminaren Grenzfall $/»o = 1 — z/s 
oder 


entsprechend (s/a) Re—>0 [13, Bild 5] zustrebt. Das theoretische Ergebnis nach [3] ist damit 
experimentell nicht nachweisbar, zumal auch für solche s/a, bei denen entsprechend dem späteren 
Bild 10 nachweislich die Grenzschicht an der ganzen Scheibe turbulent ist, die Reibung am 
Umfang der Kernströmung von nicht mehr vernachlässigbarem Einfluß ist. (Die Eliminierung 
dieses Einflusses durch eine mit der Winkelgeschwindigkeit ß rotierende Zylinderwand würde 


74 E. Broecker, Reibungswiderstand der glatten rotierenden Scheibe _ 


nicht den betrachteten Fall der im festen Gehäuse rotierenden Scheibe realisieren). Im folgenden 
seien noch die Verläufe von ß/» und Cy für die in Bild 7 und 8 gegenübergestellten Fälle theore- 
tisch untersucht: j 

Bei der allseitigbegrenzten Scheibe wird außer an der rotierenden und der festen Scheibe 
auch eine Grenzschicht e <a zwischen Kernströmung und Zylinderwand mit der Geschwindig- 
keitsverteilung 


angenommen (Bild 9), welche unter Vernachlässigung der axialen Bewegung 


7, y \1 
= 0,0225 () 

cap aßel 

ergibt. Da sich die Betrachtung auf kleinere s/a be- 
schränken soll, wird in erster Näherung die Grenz- 
schichtdicke e mit der Grenzschichtdicke ö, am Um- 
fang der rotierenden Scheibe in Beziehung gesetzt, 
welche nach (16) zu 


°e Rei = [,(B]o) 


ermittelt wurde, wobei für die Größe von e die Win- 
kelgeschwindigkeit £ statt der Differenz (o—P) be- 
stimmend ist, so daß 


Bild 9. Geschwindigkeitsprofile am Umfang der Kern- 


= ae 2 € ar 
t bei allseitig b ter bzw. seitlich b t sr 2 =: 
strömung bei allsei ee _ ch begrenzter . Re l Ta (1 ß /®) fa ( B /o) 


gesetzt wird; f, ist danach die an der durch ß/» = 1/2 gehenden Senkrechten gespiegelte Funk- 
tion f, in Bild3. Mit M,=2n@sr, als Reibungsmoment der Kernströmung an der zylin- 
drischen Begrenzung wird 


Cr = 47-0025 (Bo) (fe) Mesa 2». ou ee . (19). 


Unterscheiden weiter Cy, bzw. Cy, die reduzierten Momentenbeiwerte der rotierenden bzw. 
festen Scheibe, so ergibt sich aus Cy, = Cyr — Cyr 


sja=3,4A(oı (FR — CH) >. 2 en... 


worin nach (11) bzw. (12) u. (15) Cyr = 0,0728 - f(P/o) und ferner Cy; = 0,116 - (8/®)95 sind; 
die Konstante vor dem letzteren Ausdruck ist dabei unter Berücksichtigung von (18) eingesetzt, 
während (13) den kleineren Wert 0,1038 ergibt. 


Bei der nur seitlich begrenzten Scheibe wird nach Bild 9 von einer radialen Ver- 
mischungsbreite 7 zwischen Kernströmung und Außenraum ausgegangen, in welcher die Ge- 
schwindigkeitsverteilung 

u, Ka 
Euer ae RN 
ap n 


herrsche und für welche das Schubspannungsgesetz (y = 0) 
T; „[ » \4 

v_ 0,0205 ( 5 ) 

o(a B)° aß 


gelte. Wegen der endlichen axialen Breite der freien Mischströmung zwischen Kern und Außen- 
raum wird in den Ansatz für 7 noch das Abstandsverhältnis s/a mit einbezogen entsprechend 


Ha 
ja aß - F(s/a) ; 


eine befriedigende Übereinstimmung mit den Meßwerten wurde erhalten mit F-14 — 15 . 16,85 s/a 
9 


wonach der reduzierte Momentenbeiwert der Vermischung der Kernströmung mit dem Außen- 
raum zu 


Cy, = 47: 0,0225 (Bw) 15. edle sa nn (21a) 
und danach 


sja » 8510 = 0,236 (@JB)8 (Cr — Cup) = 0,236 [Cyr,(o/P) — 0,116]. . . (2b) 
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werden. In einem zweiten Rechnungsgang mit dem Ansatz 
y \15 
na=(5) -F*6le) 
konnten die Meßwerte ebenfalls, da Bl®o bei gegbener Anordnung nur eine Funktion von s/a 
ist, dargestellt werden, wenn F*1/® — 0,164 . (B/w)!3/10 gesetzt wurde, was anderseits 


Az - 0,0225 
Cyv = — sea  Ploy" - sja See ca DE: (22a) 


sowie f 
SIE VS OP) Ca OR ER (22b) 
liefert. Aus (21a) und (22a) folgt auch 


1 1 
— __ 0712,96 8fa m __ 5713 8/a 
Do —= ze un ze 


Nach Bild 7 erfolgt die Abnahme von ß/® mit zunehmendem s/a für die nur seitlich begrenzte 
Scheibe erheblich schneller als bei der allseitig begrenzten Scheibe, der Vergleich von (19) mit 
(21a) bzw. (22a) ergibt Ja—0 Cu, 84 Cyr- 

Die berechnete Winkelgeschwindigkeit der Kernströmung für die-allseitig begrenzte Scheibe 
(20) wird durch die Messungen nach [3] und [7] sehr gut bestätigt; die Messungen von [12] 
für die seitlich begrenzte Scheibe folgen befriedigend einer der angeführten Gln. (21b), (22b) 
und (23), so daß die Extrapolation auf größere s/a berechtigt erscheint. 

Der berechnete Verlauf der reduzierten Momentenbeiwerte für die betrachteten Fälle 
unterscheidet sich nach Bild 8 derart, daß für die nur seitlich begrenzte Scheibe bereits bei 
kleinem s/a eine Übereinstimmung mit der Messung erreicht wird, während für die allseitig 
begrenzte Scheibe erst für größere s/a Messung und Rechnung zusammentreffen. Der Grund 
hierfür kann darin gesehen werden, daß im ersteren Fall der in der Grenzschicht geförderte Volum- 
strom am Umfang der rotierenden Scheibe 


B 
V, _2nalo: dz= 2r.a 0,409 08,0, = 27: 0,409. f, fa @® @ - Renil® 
ö 


in den Außenraum austritt und ein gleichgroßer Volumstrom am Umfang der festen Scheibe 
eintritt, wogegen beim letzteren Fall ein Volumumlauf im Bereich der Grenzschichten vorhanden 
ist und insbesondere bei kleinem s/a die Eckenströmungen nicht mehr gegen die Grenzschicht- 
strömungen an der rotierenden und festen Scheibe sowie am Zylindermantel vernachlässigbar 
sind. Für die nur seitlich begrenzte Scheibe läßt sich denn auch der Übergang von der turbulenten 
Grenzschichtströmung in die Laminarströmung an der Abweichung der Messung von der Rechnung 
besser erkennen: Bezeichnet (s/a); das kritische Abstandsverhältnis, für welches die turbulenten 
Grenzschichten an der rotierenden und festen Scheibe am Scheibenumfang zusammentreffen, so 
liefert die Abschätzung 
(sd) 7 2 6./a = 2 Fz(Blo) Ret® 


für die Messung nach Bild 8 mit der dort angegebenen Reynoldszahl 
2 
(s/a); X 18,35. Plo) 3 


aus Bild 10, das den Bereich kleiner s/a von Bild 8 näher darstellt, erkennt man die beginnende 
Abweichung zwischen Messung und Rechnung für s/a X 0,015, wofür nach Bild 7 ß/o = 0,42 
und nach Bild 3 f,($/») = 0,14 ermittelt werden, was (s/a); X 0,0153 übereinstimmend ergibt. 

Für kleinere s/a nimmt Cy erheblich zu und erst für sehr kleine s/a bei über die ganze 
Scheibe sich erstreckender Laminarströmung nähertsich C,, dem theoretischen Verlauf. Letzterer 
ist sowohl nach der approximativen Lösung nach [3] als auch nach der exakten Lösung nach 
[9] und [13] eingetragen; für sehr kleine s/a verschwindet der dargestellte Unterschied der 
Lösungen, wobei der in [13] angegebene Beiwert (10.2) 


ee 


v 


log A > log| 


eht und mit cy Re=naj/s als 0. Näherung von [3] identisch wird, für größere s/a geht nach 
[13, Bild 11] cy Re? — x n/2 (vgl. Bild 4) bzw. Cy — % 0,02 in Bild 10. Schließlich finden 
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die geringen Abweichungen der gemessenen von den berechneten C “ für die nur seitlich begrenzte 
Scheibe bei größeren n in Bild 8 darin ihre Erklärung, daß der in Bild 7 dargestellte Verlauf 
von ß/®, welcher der Rechnung in Bild 8 zugrunde liegt, aus dem (gestrichelt verlängerten) 
linearen Druckverlauf nach Bild 6 bestimmt ist, wobei die Änderung des Strömungszustandes 
in der Nähe der Drehachse bei größeren s/a nicht berücksichtigt wird; bei Beachtung dieser 
Abweichung ergäben sich hierfür im Mittel größere f/o und entsprechend kleinere Cy im Sinne 
des experimentell festgestellten Verlaufs. 


Re = 2,075:10® 


Messung [4] [7] 
| 
g ' “ 
M 
| Messung [12] 
| Laminarstromung- = |_„ /unbulente Grenzschicht 


Übergang yynbulente Grenzschicht 


. ‚heibe 
H Rotierende Sc, 
| seitlich begrenzt 


y ale begrenzt 
.(1.u.23)__ _ = z— 


0 007 002 — 5s/a 003 


Bild 10. Reduzierte Momentenbeiwerte im Übergangsbereich laminar-turbulent 


1: Laminare Grenzschicht [3] 
2a bzw. 2b: Laminarströmung [3] (0. bzw. 2. Näherung) 
3: Laminarströmung [9] [13] mit Übergang in lam. Grenzschicht 
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KLEINE MITTEILUNGEN 


Zur Berechnung des spannungsvermindernden 
Einflusses einer Nabe bei der rotierenden ko- 
nischen Scheibe 


Mittels der von Honegger, Zürich!), angegebenen 
Festigkeitsberechnung rotierender konischer Scheiben 
läßt sich ein spannungsvermindernder Einfluß einer 
Nabe zunächst nur so ermitteln, daß man diese 


1) Z. angew. Math. Mech. 7 (1927), 8. 120—123 u. Biezeno-Gram- 
mel, Technische Dynamik. 


Rechnung für mehrere angenommene Außenzug- 
spannungen 0; durchführt und die Spannung am 
Innenrand in Abhängigkeit hiervon aufträgt Ebenso 
müßte dann der andere Ausdruck für diese Spannung 
(s. Hütte Bd. 1, S. 734) 
,w®r3 sr 
0 SE +0 (u 7 nn 

in Abhängigkeit von o; aufgetragen werden. Der 
Schnittpunkt ergäbe dann die gesuchte Spannung o,;. 


Ba % 
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Um diese nicht unerhebliche Rechenarbeit abzu- 
kürzen wird ein übersichtliches rasch konvergierendes 


'Näherungsverfahren vorgeschlagen, welches von dem 


für die ebene Scheibe mit Kranz und Nabe (s. Hütte) 
gültigen Gleichungssystem ausgeht. 

Für die dort auftretende Außenzugspannung o; 
läßt sich der Ausdruck errechnen 


.2 R2 2 y2 
Han =i—=0) + 0% Sn Aa 


a 98 

re 5 Bon 4 
N 
Prpatmon 


Hierin bedeuten die ersten beiden Glieder im Zähler 
die Spannung am Innenrand bei Vorhandensein einer 
Zugbelastung o, am Außenrand und der dritte Aus- 


druck im Nenner den Einfluß, mit dem die auf- 


tretende Außenzugspannung o; wegen Vorhanden- 
Fe Nabe die Spannung am Innenrand herab- 
rückt. 


Man berechnet daher zunächst für die gegebene 
Zugbelastung o, am Außenrand die Spannung am 


Innenrand der konischen Scheibe 5,;(0,) und damit 
einen besseren Wert von o; aus: 


> w?2r? 

Enlo) 
.Gi = Fi 
R sn, Ktr 
er, 


Mit diesem Wert o, rechnet man einen ersten Nähe- 
rungswert ol; (0a, 05) für die konische Scheibe aus, 
der sich analog der ebenen Scheibe auch schreiben 
läßt, 

0,; (0a, 01) = Cri(0a) — € Gi, > 


wobei e der unbekannte konische Einfluß ist, der dem 
a R?2 2 
Quotient —_ 2 

ot, 
Aus dieser Beziehung folgt e und mit diesem Wert 
wiederum ein besserer Wert von o; 


yw?r, 
oo = 


g 
2 Ar) 


ey he 


bei der ebenen Scheibe entspricht. 


9,,(04) — 


Hiermit findet sich ein besserer Wert für die Span- 
nung 6//(o,,0;) am Innenrand Dies ließe sich mit 
der Beziehung 

91 (09, 03) = Sri(00) — E Ci, 


solange fortsetzen, bis o; und 9;;(0,, 0;) sich nicht 
weiter ändern, jedoch ist dies meistens nicht not- 
wendig, weil das Verfahren sehr schnell konvergiert 
und schon der zweite Wert [74 (0a, 0) bzw. Gi, sich 


von dem ersten nicht sehr unterscheidet. 


Verfasser: Dr.-Ing. J. Meyer, Kassel, Kölnische 
Str. 95 1/2 


On the Extension of W. Nusselt’s Heat 
Transfer Theory for the Film Condensation 
of Steam 


W. Nusselt investigated the heat transfer for the 
film condensation of steam by using the approximate 
theory [1], but this theory is not applicable for the 
case of steam condensation on the flat horizontal 
plate having definite extension. 

Here the author wish to describe the extension of 
W. Nusselt’s theory under the consideration of the 
existence of the statie pressure distributions in the 
water film of condensed steam. 

The Navier-Stokes equation of motion can be 


ö 
written as follows under the assumption of = =, 


which shows that ”’there are no accelerations in the 
water flow in the film“, 
\ dp du 
ynO— +taga=0.... 
where 
y spezific weight of water in th= film 
p static pressure in the water film 
u velocity of water in the film 
© inclination angle of the surface of the cooling 
plate to the horizontal plane (small angle) 
4 coefficient of viscocity of water in the film. 


Using the water depth Y and Y, in the water film 
as shown in Fig.1, the static pressure p can be written 
as follows 


dp __dY, 3 
P=29%+Yy(% 3% TE 0-0 
where p; shows the pressure of saturated steam. 
„2 
Sarurated' 


coo/ing 
Surface 


Bild 1 


Using above relation, the equation (l) can be 
written as follows 


: dY, deu 
a — 0 a2): 
By integrating above equation (2) for y, the velocity u 


can be written as follows under the boundary con- 


d 
dition of (1) u = Oaby= Oand (2) 7, = Oaty= 9. 
Era -| = 

27, sin © Y Yo 9 a) (3). 


The mean velocity can be calculated as follows 
Yo dY 5 
=. Er Ei erh I 
in, eu, de sin © 3° 
The quantity of water flow for unit width of the 
water film is 
= a Fz i je 
G =UmYyıYy = Air sin © 3 
and further we can write as follows 
alle ] 
dG, = rnler sin 0) Yaıde. 
On the other hand, the condensed quantity is 


for unit width of the water film, where 


A  coefficient of thermal conductivity 
r _latent heat of vaporization of water 


0; temperature of saturated steam \ 
@,u temperature of the surface of cooling plate 


The relation of d@, = d@, can be written as follows 
Fl mo)n]= 4-98 


3u del\de Yo r 
In the Fig. 1, we can write as follows 
ee ar, __1 m 
0 008 9’ de cos® de 


Using above relation, the equation of (4) can be 
written as follows 


av Mm. ‚ ı|=-22:0,- 0 A 
sad sin © | = 5, (9% u) 
This differential equation can be written as follows 

by using the constant values of A, Band C 


a (avi) a 
ee 


de "Y 

where i 

4.0080’ ae se 3 

f6) 3 u (Ou— 95) 
Pr 
This differential equation can be written as follows by 
tti a ah dh dag 

ie et air ne 

An I. BihueGg OH —=0, ee (5). 


The differential equation of this type will be solved 
by the graphical and the numerical method as shown 
in the book edited by N. W. MeLachlan [2]. 

We can write many different q-h curves for the 
given different values of dg/dh, and we can choose the 
solution of the equation of (5) by the given boundary 
conditions.. 

The film condensation of steam on the definite 
horizontal plane was investigated theoretically by 
S. Takatama [3]. The author wish to express many 
thanks to Mr. S. Takatama, who is the author of 
the most useful paper for writing this paper. 
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Elementarer Beweis für die Produktent- 
wicklung des Sinus und die Partialbruchzer- 
legung des Cotangens 

1. Wir stützen uns im folgenden auf die elementaren 
Ungleichungen und Formeln: 


x 


x 7 de T 
Mer <-——-s 7 für 0 <u <mSs = 
sin%; SINX, 2 ER 
b (1); 
1l=lim-— für 20 (#0) 
sin x 
ir >, fürr0O <r< & 5) 
zZcar 2 9 ur 0osıs Bu > (2); 
1 a3 1 SUR } 
- — < DES Sn 
2» y! po—1 pz=» =4q ( )s 
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‚g natürlich; 2=p», q 2e2a—l (denn es ist 
5 B. für den angegebenen Zeigerbereich 


1 
Di>S,.cn- 26H 


1 1 HA eRIRL 2 
ar DR a EB 
/ 1 
e2: sl—xz se für Izr=25 (4); 
RE na 
sin« = 2sin > cos. — 2sinz- sin z-- (5); 
sinz-sinp= sinn + P in TR 23" FOR 


2. (5) mit dem +-Zeichen rechts ergibt für = x 


a 2 N zs+l1l 
sinszx = 2sinz Z sin z 2 


mE 


(7) auf jeden Faktor rechts angewendet, liefert 


} ge Ef N A a 

sinzxz = 2 sin 2," Sinn ng 

s+2 . z+3_ 
23 - sin # PIE 


X sinz 


Wiederholung dieser Umformung ergibt nach (n—2) 
weiteren Schritten 


0,20 —1 


: n_ II 4 z+rv 
sinn x = 2? 1 sın x 
v 


9° 


Bis auf den Faktor x fungieren rechts die Argumente 


z|ic+v za M12 42 —y 
aa m a 2" 
»—=172,..:202=1] = 1,2,..,2H RT 
Es ist 
. (2 + 2n) —y z v—xt 
— =sın3 
sin z e sin z = 
ferner 
; x x + 20-1 
sın m In sın 7 — uns 
. x x 2 x 
=sinz,,' c087 a} sin x Ina 
Mithin: 
sinztx 
1,291] 
on, 1. 8 1 ‚ z(w+2) . 2-2) 
=2 5 Sin an Bin —— an af: 
= 2 


Die {} darin ist nach (6) 


sin? — 

EN‘ ‚9 IV 6% ru 2n 
or sın“ — sin? = sin? ee 1,0 \ 

n 9n yn nv 

. ä s sin? — 

2’n 


Nimmt man (8) für 0O<x< 1 und dividiert beide 
Seiten durch x, so liefert der Limes x — 0 nach (1) 
ee 
20 — Yan a 


9709 
sin? zu > 
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wodurch sich (8) so schreibt: 


IT 
en on | sin? >= 
Bin. a Gr 202. sin 1l— 9). 
vi l Sag) 
sın an 
Führt man rechts gliedweise den Limes x — 0 aus, 


so entsteht das unendliche Produkt 


Wir setzen 


Sei jetzt x >ganze Zahl und k, so gewählt, daß 


IT 
Be] 
rss 3’ k>2 >I Re) a (10); 


ferner seien im folgenden k, n beliebig, jedoch so, daß 
kZk, und 2712 —-1>2(k+1)—1, alon=5 (11) 
ist, und schließlich sei 


97% 
Ber, Ban 
1) sin an J] De rn 
v sin? — 
NZi-, Terier 
en k+1,2 1 Sn. In 
Pr(x) = II 1— zu 
in? ——_ 
& sın 9n 


Dann ist 
Prix) _ Pr) 1 


= (n) 7) 2 
NR) 


) 
Wir schätzen Pelz) nach oben und unten ab. 
Wegen (10) und (1) ist 


sin T® «| 
n n 
- - z “ . (12), 
Ber Fr 
cr 1 gan 
Bea er ur 
an 2r [a] 2” je| (13) 
= p) > 
sin Sn ins, = sin SE A 
und also wegen (4) und (5) 
ren 
sin? zz 22 Een 
er ET (14), 
sin? T# e 
9n 
FE = x) ( 
1 Zudeh AL = 1 2 ’ > e 2 [2 (15), 
sin?” & = 
PAL 
mithin wegen (11) 
Ba ln . 
0) z w ; 
Pr(x) < e a a 
ae | 
m a re 
Pr(x) > e i > . 
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Man hat also 
x* ; aa, 
ar, [Pro] _ Pie) [2 u 
- | (m) sinzz in . x 
Pr): Pr) 


Führt man darin den Limes n — oo aus, so wird 
lim []=1, und man erhält 
a: 2 2° 
SER Prix 
et) < - *(x) <e2n, 
sin m x 


Der Limes k > oo ergibt jetzt 


= oder lim P;(2) =sinnz, 
ko 


. u x2 
sinze=a2: ]]]1-% ..„>2.(16), 
v 


gültig für jedes x ganze Zahl; ist x = ganze Zahl 
so ist die linke und die rechte Seite Null, so daß 
(16) für jedes x gilt. 

8. Sei jetzt 0 <x<1l. Dann darf man von (9) 
den Logarithmus nehmen, und die anschließende 
Differentiation liefert 


T% 
ER Dar = 
neigne—= . 
2 2nT 
sin zn: 
ae, TORE a6 
em De E72 
Am re 
1,221 sin? 
2n 
u — ku 
9X 
z Sn zu 
1 
97V 
sin In 


Führt man rechts den Limes’n — 00 gliedweise aus, 
so ergibt sich die unendliche Reihe 


1, © 


A 


rv 


Wir setzen 
ok 
1 Se 
Sr(lz) = 
ia) = — > Fr rap 
r 
TE, 00 TER 
sin = s 
9 ® 
RE “= 97V 
m ee = ’ sın an 
EZ zz m > LRT- ..n® ” 
0% gg 
sin mi 5 sin? 
ee 
97V 
sın DI 
: ne 5% 
sın In cos Im In 
m ee = 
2 u 
el sin? Zn j 
aa) = — as) 
: sin? 5 
I —PE Tu 
gl 
sin? 
Dann ist 


(n) e (n) 
Sy(x) — a ctg a © = [Irl2) — SHl@)] + (— s(x)) (19). 


80 


(k) n 
Wir schätzen — s;(x) nach oben und unten ab. Für 
die {} in (18) ergibt sich wegen (12), (13), (1), (2) 
und (10) und (11): 


nez mu 
u nx 1 
tr TE BIOR I 
sin 5 sin zn 
7 =. 
ur RR NE Le s;4 
u 2 um" 
na nu 
non a u | 
Bi nee erg 
sin zn sin zn 
I: 
ER S 1022 
u 2 u u 


Damit ergibt sich für den einzelnen Summanden [ ] 
in (18) die Abschätzung 


2 
Dr 2 2 
_ z® —® 
= ee Bee 
2 52 
Fe 


x 
ie 85 x 
U BT wm—ar "u 
7 
Man bekommt wegen (3): 
Kr, 2R 1 
m) 12 1 
— (2) < Z ie 
262) 5? : rpy 
® a up! = i 
=. BR An 
7 
en kHLaNI—ı | 1 
— srIT . 
ER - aa 
und weiter aus (19) 
RS " ia —ın en N 
I an) > € r ja € ” 
Par ee) > #2 


Der Limes n — 00 liefert wegen lim [ ] = 0 


ee r "x 1 
nl — nctgma}s 
EN = {S,() netgna,s 5 L—3’ 
und der Limes k — 


lim {S;2) — mn ctgmx}—=0 oder lim S.() = rn cetgnx 


ko IS 
d.h. 
1,0 d 
negne= — + Be ah 
1, 
1 | 1 1 
=. D) 
2 rt ut « (20), 
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gültig für 0<x<1l. Ersetzt man rechts x durch 
z+1, so konvergiert — wie sofort zu sehen = die 
neue Reihe wieder, und zwar zum selben Limes, 
m.a. W. auch die rechte Seite hat, ebenso wie die 
linke, die Periode 1. Mithin gilt (20) für jedes -# 
ganze Zahl. 

Verfasser: Prof. Dr. E.Mohr, Berlin-Wannsee, Tri- 

stanstr. 1 


Verfahren zur zeichnerischen Ermittl des 


äquatorialen Trägheitsmomentes —. rmi- 
ger Flächen 


Das Trapez abcd hat die Grundlinien b und B. Der 
Abstand e des Schwerpunktes wird in bekannter 
Weise bestimmt. Er liegt auf der Achse 0 —0. 
Parallel zu ihr zieht man im Abstand e von ab die 
Linie <&— x, die die Seiten adund be in den Punkten 
fund g schneidet. Die Strecke fg = b’ ist die Grund- 
linie eines Rechtecks mit der Höhe h. Es wird be- 
wiesen, daß das Trägheitsmoment Jr des Rechtecks, 
bezogen auf dessen Mittelachse gleich dem Trägheits- 
moment Jr der Trapezfläche, bezogen auf dessen 
Schwerachse 0—0 ist. Aus der Zeichnung ergibt sich 
mit B—b=b, 


B—b b_W—b, ya Y 
a u Be as EZ (1). 
Setzen wir für y=h—e und für e= ee 
(bekannte Formel) dann wird 3 2b+b 
eo 3b +2b, 
Umgeformt ergibt sich 
2 2 
he nn aa © (2). 


3(25+b,) 
Dieser Wert in die Formel für Jr eingesetzt: 
m 65 +65b, +}, 

7" Berne 
Der rechte Teil der Gl. (3) ist aber die Formel für Jr. 


Jr bezogen auf die Mittelachse parallel ab. ist also 
gleich Jr bezogen auf die Schwerachse 0 — 0. 

Das zeichnerische Verfahren bietet zweifellos Vor- 
teile. Das Maß b’ läßt sich auf dem Reißbrett schnell 
mit genügender Genauigkeit bestimmen und hierauf 
Jr = Jr ebenso schnell ausrechnen. 


Verfasser: Karl Borth, Dresden A 29, Unkers- 
dorfer Str. 10 


Ir= 


Zur numerische 
Funktionen 


In Z. angew. Math. Mech. 36 bemerkte W. Loh- 
mann [1] folgende Eigenschaft periodischer Funk- 
tionen: Integriert man eine genügend oft differenzier- 
bare periodische Funktion f(x) = f(x + p) über eine 
volle Periode p numerisch, so kann man bei An- 
wendung der Rechteckregel den dabei entstehenden 
Fehler mit einer Fehlerabschätzung der Simpsonregel 
abschätzen. 


Integration periodischer 


. ER ae 5 
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Hier soll gezeigt werden, daß die Rechteckregel, 


angewandt zur numerischen Integration periodischer 


Funktionen, im allgemeinen noch wesentlich schärfere 
Fehlerabschätzungen erlaubt. 


In diesem Zusammenhang sei zunächst auf ein 
Lemma zur numerischen Integration einer perio- 
dischen Funktion von G. Birkhoff, D.M. Young 
und E. H. Zarantonello [2] hingewiesen, das aus- 
sagt: Sei F(g) eine periodische Funktion der Periode 
2x, so werde das Intervall <0,2x)> durch die 
Punkte q,, 9»...,.qgv in N gleiche Teilintervalle 
unterteilt. Dann ist die Trapezregel 


= 2 


wa 
Fddr 2 Fa) dgr 
{ 


eine „bestmögliche“. 

‚„Bestmöglich‘ ist hier in folgendem Sinne zu ver- 
stehen: Führt man die numerische Integration mit 
Hilfe einer beliebigen Quadraturformel 


2r 

. 

J Fla)dg x Dar Fig), Pia) = Plan), 
2 Il) 

U 


N-mal durch, wobei als Anfangspunkte die Teilpunkte 
91» 92 » - - Qv genommen werden, dann ist der Mittel- 
wert der so erhaltenen bestimmten Integrale mit 
einem Fehler behaftet, der das arithmetische Mittel 
der bei den N Integrationen entstehenden Einzel- 
fehler ist. Der Fehler bei Verwendung der Rechteck- 
regel ist ebenfalls gleich dem arithmetischen Mittel 
der Fehler, die bei N-maliger Verwendung einer 
(N + 1)-Punkteregel entstehen. Dazu wird lediglich 
die Tatsache benötigt, daß Fa, = 2 ist, eine For- 
derung, die ja nur besagt, daß die Funktion F(g) =1 
exakt integriert wird. Der Beweis dieses Lemmas 
verläuft analog wie die im folgenden zu zeigende 
Erweiterung. 

Unsere weiteren Betrachtungen beziehen sich auf 
solche Abschätzungen des Fehlers R einer Quadratur- 
formel & x; F(g;), die von der Form 


|R| < © Max |&(F)| 
sind; die bekanntesten Fehlerabschätzungen gehören 


zu dieser Klasse. Dabei ist das lineare Funktional 
Ö(F) etwa 


d’F 
= - <y<s 
®(F) Ber Isv<N+l 
oder (Ph. Davis [3]) 
DIE) His 


Im letzteren Fall ist F in einem geeigneten Gebiet 
zu betrachten, in das die Funktion F analytisch fort- 
gesetzt wird. 

Sind nun R” die Fehler der bei g, beginnenden 
numerischen Integration, so ist 

| 


1 ee 
n23 RO = w2 |R®|< Max |R®|< C Max |@(F) 


die Abschätzung des arithmetischen Mittels der Fehler 
bei Verwendung einer (N + 1)-Punkteregel, und für 
die Einzelfehler gilt ebenfalls 


|R®| < C Max |&(F)|. 


Nach dem genannten Lemma gilt dieselbe Abschät- 
zung auch für den Fehler bei Verwendung der Recht- 
eckregel. 

Deshalb enthält dieses Lemma auch die folgende 
Feststellung: Die Rechteckregel — die ja bei der 
Integration einer periodischen Funktion mit der 
Trapezregel identisch ist — bringt einen Fehler mit 
sich, der mit derjenigen Fehlerabschätzung ab- 


geschätzt werden kann, die für eine beliebige Qua- 
draturformel der Gestalt 


2x 
[2 N 
IE ig Zap Flq). 


gilt. hy BER 

‚Die Aussage von W. Lohmann ist für N=2 
hierin enthalten. Für N >2 besagt die obige Aus- 
sage mehr, da die Fehlerabschätzung einer (N + l)- 
Punkteregel im allgemeinen engere Schranken liefern 
wird als die der Simpsonregel. en 

Es ist jedoch möglich, hier eine weitergehende 
Aussage zu gewinnen, die im allgemeinen noch engere 
Schranken liefern wird. Betrachten wir nämlich auch 
Quadraturformeln Hermitescher Art, die. neben den 


- Funktionswerten an den Stützstellen auch die Ab- 


leitungen heranziehen, also von der Form 
%+P 


Sy ıyP 
J Ka)da= I atom + kB) 


+ Zar (2 =. EP) + Bü | 


n 


= N (ar fe + ayf;) + En 
la 
sind; dabei ist R„ = O(n) D(f). 
Dann wird 


%+Pp i 
IT ))da=ahtagh+'' +anfn 
e +tahtrteahkt tom tn. 


n-+1i 
To oe, 
I= (ie ))ak=-mhtauf + aufn 
4 +aftafte tat Re 
2n—l 
: a 
I= 3) )da=ah tft: + aı fn 


7 PE 
k 


nI= Eh2y +22, +ERP ). 
k 


Ebenso wie man von einer linearen Quadratur- 
formel der Gestalt 


b 


n 
IK de > > 0% Pk 
Fi v 
mindestens verlangt, daß sie die Funktion o(x) = c 
exakt integriert, soll eine Formel 


b 
| plz) de =D) (ar Pr + 2} Pr) 


) 
a 


mindestens den Wert des Integrals über 


pa)=at+%r 
exakt liefern. 
Daraus folgt zunächst wieder, wenn wir der Ein- 
fachheit halber « = 0 annehmen, I ar = b. 
Es ist 
b 


p2 
Jatana@=ar+az 


0) 
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e ! f 
Q =. (X Pk + &% Pr) 
0 
= ZT 2 +2 + +n—1)a,.ıtNn&n} 
n 
+6 & an. 


Wir machen nun Gebrauch von der Eigenschaft, 
daß in einer Quadraturformel allgemeiner Art die 
Funktionswerte mit symmetrischen Gewichten ver- 
sehen sind, also &; = «&,_;ist. Dann ist 


N 
g=uadb+aNa 
0 


b n n n 
nt Fan t -+yatry 2} 
VD 
=u4ub+,d% az 
f) 


‚0 d? er 
=ab+, ta %- 
0 


Daraus folgt £«, = 0, und wir erhalten aus (*) 
nI=-pZh + ERD, 


ud EL. 
I=Sht BR. 
Nun ist R® = C(n) ©) | 
IR] = |C(n)| Max |O(/)| 
un 
2 7 
I= n 2 lh + En 
mit 4 
|En| < C(n) |Max |@(f)| ; 


hier ist etwa 
dn+2 7 


—  dan+2 


DU) 


Das ist aber die Rechteckregel mit der Fehlerab- 
schätzung einer Hermiteschen Quadraturformel für 
(n + 1) Stützstellen. 

Ergebnis: Integriert man eine periodische Funktion 
f(x) = f(x + p) über eine volle Periodenlänge p mit 
der Rechteckregel unter Verwendung der Funktions- 
werte an n äquidistanten Stützstellen, so kann man 
den dabei entstehenden Fehler mit derjenigen Fehler- 
abschätzung abschätzen, die bei Verwendung einer 
Hermiteschen Quadraturformel mit (n + 1) Stütz- 
stellen gilt. Die Aussage bezieht sich dabei auf solche 
allgemeinen Hermiteschen (r -+ 1)-Punkteregeln, in 
denen die Gewichte der Funktionswerte symmetrisch 
zur Intervallmitte sind. 

Es ist der besonderen Bemerkung wert, daß damit 

. freilich nicht ausgesagt wird, die Rechteckregel 
liefere durchweg ebenso genaue Ergebnisse wie eine 
Quadraturformel höheren Approximationsgrades. Die 
höhere Genauigkeit einer solchen Formel ist aber des- 
halb nur von beschränktem Nutzen, weil keine ge- 
naueren Fehlerschranken als bei der Rechteckregel 
zur Verfügung stehen. 
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Zur Identität von Gleitlinien und Charak- 
teristiken 2 

Es interessiere hier eine durch eine Mohrsche 
Hüllkurve gegebene Plastizitätsbedingung 


Tm=k+mtge, 


bei der Kohäsion k und Gleitwinkel g sich in Ab- 
hängigkeit von der Schubspannung ändern. Die 
Tangenten an die Hüllkurve sind dann gegeben durch 


w=k+oitge, re 
wobei 

k= Km), e= Mm): 
Funktionen der Schubspannung „im Berü punkt 


von Tangente und Mohrscher Hüllkurve sind. Be- 
deutet ein beigefügter Strich die Ableitung nach dem 


“ Argument, so lautet die Parameterdarstellung 


der Mohrschen Hüllkurve: 
k’ cos? go 


EEE ER ( 
n r ) 
mit der Identität 
we a BE: - 
14 
Bezeichnet man nun die unter den Winkeln I _ z 


gegen die größte Hauptspannung gerichteten Kurven 
als Gleitlinien, so lautet mit k, = k cotg o und 
K=/kd(lnk), 
sowie 
_ At 
2 
die Bedingung dafür, daß die Gleitlinien mit den 


Charakteristiken des auf sie bezogenen Systems der 
Gleichgewichtsbedingungen identisch sind!): 


v=e0o7, o 


+ kı 


8t=F+i=0. 


Daß diese Bedingung stets erfüllt ist, sieht man ein, 
wenn man sie in der Form schreibt 


ktge=—eo’o 
oder 
R 
K ( = 2 ge: 
er sin 0 cos @ e . nn 
Aus 
si an = sing 
2 
und 
Tn= n 9— cos 0 
folgt 
Tn 


 sin_coso ’ 


womit (2) in die Gleichung (la) übergeht, also stets 
erfüllt ist. 
Aus (2) entnimmt man, daß bei Auffassung von k 


als Funktion von o längs der Hüllkurve zu festem 0 
die halbe Spannungssumme 


gehört. 


!)H.Schlechtweg, Z. angew. Math. Mech. 38 (1958), S. 139/148. 


Verfasser: Prof. Dr. H. Schlechtweg, Kettwig/ 
Ruhr, Brederbachstr. 13 
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BUCHBESPRECHUNGEN 


Dr.-Ing. A.Weigand, Einführung indie Berech- 
nung mechanischer Schwingungen. Bandll. 
176 S.m.123 Abb. Berlin 1958. VEB Verlag Technik. 
Preis 21,— DM. 

Der nunmehr vorliegende zweite Band enthält die 
Berechnung der Schwingungen von mechanischen 
Systemen mit zwei und n Freiheitsgraden, wobei je- 
weils die freien ungedämpften und gedämpften, die 
erzwungenen Schwingungen und Einschaltvorgänge 
sehr sorgfältig, und durch klare Abbildungen illu- 
striert, behandelt werden. Bei der Berechnung des 
Einschwingvorganges wird von der Laplace-Trans- 
formation Gebrauch gemacht. 


Freiberg/Sa. D. Rüdiger 


Dr. W. Macke (Prof. a. d. TH Dresden), Wellen. 
Ein Lehrbuch der theoretischen Physik. XII + 4658. 
m. 160 Abb. Leipzig 1958. Akademische Verlags. 
gesellschaft Geest & Portig K.-G. Preis geb. 29,50 DM, 


Das Buch eröffnet, wie der Verfasser ankündigt, 
eine Lehrbuchreihe über theoretische Physik, der 
eine vom bisherigen Gebrauch völlig abweichende 
Einteilung des Lehrstoffes zugrundeliegt. Während 
die meisten Lehrbücher die auf der Sinneswahrneh- 
mung beruhende, klassische Einteilung der Physik 
in Mechanik, Akustik usw. bevorzugen, bilden hier 
gewisse, für den gegenwärtigen Erkenntnisstand 
grundlegende Begriffe, wie etwa der Wellenbegriff 
das ordnende Prinzip. s 

Es zeigt sich, daß die Verwendung des in der Wel- 
lengleichung formulierten Wellenbegriffes als Leit- 
motiv eine sehr vorteilhafte Methode darstellt. Sie 
erleichtert das Eindringen in die Fülle des Stoffes 
und fördert die Übersicht über den Zusammenhang 
zahlreicher Gebiete der Physik, die nach klassischer 
Auffassung weit auseinanderliegen. 

Während sich der Anfänger bisher die Wellen- 
gleichung praktisch mehrmals aneignen mußte, je 
nachdem in welchem phänomenologischen Gewande 
sie auftrat, wird er hier einmal gründlich mit ihr ver- 
traut. Er wächst gleichsam in den Formalismus 
hinein, der hier das Wesentliche ist und später in der 
Erinnerung das Bleibende sein wird. Das aber braucht 
gerade der künftige Experimentator. 

Mancher reizvolle Gedankengang könnte für den 
Aufbau der Vorlesung in Experimentalphysik An- 
regungen geben. 

In den ersten 3 Kapiteln geht zunächst, begin- 
nend beim anschaulichen Beispiel des harmoni- 
schen Oszillators, der Ausbau des mechanischen 
Modells und die Erweiterung der Theorie Hand in 
Hand. Es folgen Kapitel über skalare Wellen (Aku- 
stik), wobei u. a. die Beugung ausführlich behandelt 
wird, und über Strahlen als Grenzfall kurzer Wellen, 
u. a. mit der Theorie der Abbildung. Danach schafft 
eine Ergänzung des mathematischen Rüstzeuges die 
Grundlage für die Behandlung transversaler (elektro- 
magnetischer) Wellen sowie der Lichtausbreitung in 
Materie. 

Die Wellenausbreitung im bewegten Koordinaten- 
system erfordert weiter die Darstellung der speziellen 
Relativitätstheorie. Schließlich wird im letzten Kapi- 
tel nach Einführung der vierdimensionalen Vektor- 
rechnung in einer systematischen Wellentheorie dem 
Leser eine Vorstellung von der Tragweite der Theorie 
vermittelt. 

Zahlreiche Übungsaufgaben zu allen Abschnitten 
ergänzen und vertiefen den Stoff. 

Nach der Lektüre dieses Buches erwartet man mit 
Spannung die angekündigten weiteren Bände der 
Reihe über Teilchen, Quanten, Felder, Statistik, 
Relativität. 


Dresden H. Zimmer 


I. N. Bronstein und K. A. Semendjajew, Taschen - 
buch der Mathematik für Ingenieure und 
Studenten der Technischen Hochschulen. 
XII + 548 S. m. 427 Abb. Leipzig 1958. B. G. Teub- 
ner Verlagsgesellschaft. Preis geb. in Plastikfo- 
lie 22,50 DM. 


Es ist außerordentlich begrüßenswert, daß der Ver- 
lag sich der Mühe unterzogen hat, das bekannte 
Taschenbuch der Mathematik in deutscher Über- 
setzung herauszubringen. Es wird damit eine große 
Lücke im deutsch-sprachigen Schrifttum ausge- 
füllt. Das Taschenbuch überrascht den Leser durch 
die Fülle des auf relativ kleinem Raum gebotenen 
Stoffes. Die Darstellung ist für ein Taschenbuch sehr 
günstig gewählt, denn neben einer übersichtlichen 
leicht verständlichen Textgestaltung erleichtern zahl- 
reiche sorgfältig ausgewählt Beispiele sowie einwand- 
freie Abbildungen das Verständnis sehr. 

Um einen Überblick über die behandelten Stoff- 
gebiete zu vermitteln, sollen zunächst die Hauptteile 
des Buches angeführt werden: 

Teill, Tabellen und Kurven. Teil2, Elementar- 
mathematik (Näherungsrechnung, Algebra, Geometrie 
Trigonometrie). Teil3, Analytische Geometrie und 
Differentialgeometrie. Teil4, Grundzüge der Ana- 
lysis (Differentialrechnung, Integralrechnung,Diffe- 
rentialgleichungen). Teil5, Ergänzende Kapitel der 
Analysis (komplexe Zahlen- und Funktionentheorie, 
Vektorrechnung, Fourierreihen, Variationsrechnung. 
Teil6, Auswertung von Beobachtungsergebnissen 
(Grundzüge der Wahrscheinlichkeitsrechnung, empi- 
rische Formeln und Interpolation). Es wird also der 
Stoff vermittelt, der von Ingenieuren und Ingenieur- 
studenten der technischen Hochschulen sowohl beim 
Studium wie auch in der praktischen Arbeit benötigt 
wird, wobei natürlich weit mehr Stoff geboten wird 
als in den mathematischen Vorlesungen für Ingeni- 
eure an technischen Hochschulen im allgemeinen ge- 
‚bracht wird. Der Auffassung der Verlagsgesellschaft 
B.G. Teubner, daß das Buch auch für Physiker 
und Mathematiker, die häufig mit numerischen 
Rechnungen zu tun haben, ein brauchbares Nach- 
schlagewerk darstellt, kann sich der Referent nicht 
voll anschließen, da vor allem die mathematischen 
Grundlagen und Zusammenhänge erörtert werden, 
dagegen auf rein numerische Methoden nur in einigen 
Abschnitten eingegangen wird. 

Besonders hervorzuheben sind die umfangreichen 
und zweckmäßig ausgesuchten Tabellen das Werkes. 
So findet man Tabellen der Zylinder- und Kugelfunk- 
tionen sowie elliptischer Integrale. Im Abschnitt ‚‚Inte- 
gralrechnung‘“ befinden sich Integraltafeln von 515 
unbestimmten und 45 bestimmten Integralen. Es wäre 
jedoch begrüßenswert, wenn die Tabellen für Quadrate 
Kubikzahlen usw. (S. 3—21)in der Stellenzahl erwei- 
tert würden, was drucktechnisch keinerlei Schwierig- 
keiten bereiten dürfte. 

Warum der Verlag glaubte, das Werk in der 
deutschen Ausgabe durch ein Kapitel über Variations- 
rechnung vervollständigen zu müssen, ist nicht ein- 
zusehen. Mit der gleichen Berechtigung hätte ein Ab- 
schnitt über Integralgleichungen oder, um den Belan- 
gen der praktischen Mathematik mehr Rechnung zu 
tragen, ein Abschnitt über Nomographie oder 
dergleichen eingefügt werden können. Für eine weitere 
Auflage wäre es wünschenswert, wenn das Literatur- 
verzeichnis erweitert und das Sachregister noch 
stärker der deutschen Terminologie angepaßt würde. 
So findet man etwa im Sachregister das Stichwort 
„Additionstheorem‘“ nicht. Ferner wäre es günstig, 
wenn man z. B. die Bezeichnung der Begriffe ‚‚sinus- 
oidal“ und ‚‚solenoidales Vektorfeld‘‘ ändert. Die 
Anmerkung auf Seite 252 dürfte für die vorliegenden 
mathematischen Zwecke überflüssig sein. Die Defini- 
tion der unendlich kleinen Größen auf Seite 239 
ist unklar übersetzt. 
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Doch sollen diese Einwände den Gesamteindruck 
des Werkes in keiner Weise schmälern. Das Buch 
kann vor allem Ingenieuren und Ingenieurstudenten 
empfohlen werden. Dem Verlag gebührt Dank für 
die sonst gute Ausstattung des Werkes. 


Dresden H. Adler 


Dr. R. Sauer (o. Prof. f. Mathematik u. anal. Me- 
chanik a. d. TH München), Anfangswertprobleme 
beipartiellenDifferentialgleichungen. (Grund- 
lehren der mathematischen Wissenschaften, Band 
LXII). Zweite erweiterte Auflage. XVI + 284 S. m. 
68 Abb. Berlin/Göttingen/Heidelberg 1958. Springer- 
Verlag. Preis geb. 41,— DM. 

Der Inhalt der ersten Auflage des bekannten Buches 
über Anfangswertprobleme bei partiellen Differential- 
gleichungen und Differentialgleichungssystemen vom 
hyperpolischenTypus wurde im wesentlichen unver- 
ändert (hinzugekommen ist u. a. ein Abschnitt über 
Differentialgleichungen vom gemischten Typus) in 
die vorliegende zweite erweiterte Auflage übernom- 
men. Neu hinzugefügt wurde ein fast 50 Seiten langes 
Kapitel, in dem im Anschluß an eine bekannte Arbeit 
von F. Penzlin die Theorie der Distributionen (in der 
Darstellung von H. König) zur Behandlung von An- 
fangswertproblemen herangezogen wird. Der Ab- 
schnitt geht in seinen ersten beiden Paragraphen 
zunächst auf die Grundlagen des Kalküls ein (Distribu- 
tionen und ihre Ableitungen,Konvergenz von Folgen 
von Distributionen usw.) und auf die damit in sehr 
durchsichtiger Weise mögliche Behandlung von 
Sprungfunktionen (Riesz-Distributionen, Beziehun- 
gen zur Hadamardschen Verallgemeinerung der 
Riemannschen Integrationsmethode usw.). Nach 
einem weiteren Paragraphen über Faltungsgleichun- 
gen, derinsbesondere auch auf Anfangswertprobleme 
der Wellengleichung eingeht, wird der Kalkül auf 
wichtige Problemstellungen der dreidimensionalen 
Überschallströmung [Problem des flachen Körpers 
(Tragflügel),Problem des schlanken Körpers (drehsym- 


metrischer Rumpf)] angewandt. Abschließend wird, 


die Erweiterung der Laplace-Transformation auf Dis- 
tributionen und ihre Verwendung zur Lösung von 
Anfangswertproblemen dargelegt. Das schöne und 
vom Springer Verlag in gewohnter Weise vorzüglich 
ausgestattete Buch wird auch in seiner zweiten Auf- 
lage — nicht zuletzt wegen des sehr interessanten 
Kapitels über den Distributionskalkül — wieder 
viele Freunde finden. 


Dresden M. Landsberg 


Dr. K. Strubecker (o. Prof. d. Mathematik a. d. 
TH Karlsruhe), Vorlesungen über Darstellende 
Geometrie. (Mathematische Lehrbücher, Bd. 12) 
324 S. m. 202 Abb. Göttingen 1958. Vandenhoeck & 
Ruprecht. Preis geb. 16,80 DM. 

Vorliegendes Buch ist aus langjährigen Vorlesun- 
gen des Autors über Darstellende Geometrie 
für Studierende der Mathematik und Physik, des 
Maschinenbaus und der Elektrotechnik sowie des Bau- 
ingenieurwesens und der Architektur an der TH Karls- 
ruhe hervorgegangen. Das Werk ist eine hervorragende 
Einführung in die Grundlehren der Darstellenden 
Geometrie aus der Feder eines bedeutenden Fach- 
mannes und Pädagogen. Der Verfasser bemüht sich 
dem Verstänunis der Darstellenden Geometrie den 
Weg zu bereiten und den Leser durch diese hohe 
Schule des räumlichen Denkens zu führen. Dadurch 
sind eine gewisse Breite der Darstellung und gelegent- 
liche Wiederholungen bedingt, die dem Lernenden 
aber das Erfassen des gebrachten Stoffes wesentlich 
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Nachrichten 


erleichtern, insbesondere da alle Konstruktic a 
ausführlich begründet und genau beschrieben werden. 
Da der Verfasser auch analytische Verfahren und 
Betrachtungsweisen verwendet, wird dadurch der 
innere Zusammenhang der Darstellenden Geometrie 
und der Geometrie überhaupt mit den analytischen _ 
Disziplinen der Mathematik betont, ohne daß dadurch 
die anschaulich-geometrische Seite in den _Hinter- 
grund gedrängt würde. Die Sätze aus der en- 
theorie, die gelegentlich benötigt werden, sind - 
Schlußkapitel ohne Beweis, aber unter entsprechen- Si 
den Literaturhinweisen zusammengestellt. “. 
Nach einer kurzen Einleitung mit interessanten 
historischen Bemerkungen werden zunächst die 
wichtigsten Eigenschaften der Parallelprojektion und 
der perspektiven Affinität besprochen und die Grun 
aufgaben der Lage und des Maßes im Grund- und 
Aufrißverfahren behandelt. Sodann werden eben- 
flächig begrenzte Körper betrachtet und beim ebenen 
Schnitt einer Pyramide die perspektive Kollineation 
eingeführt. Der nächste Abschnitt behandelt das 
Schrägrißverfahren, den Satz von Pohlke und die 
schiefe Axonometrie (Aufbau- und Schnellrißverfah- 
ren). Daran schließen sich ausführliche Betrachtungen 
gewisser krummer Flächen, wie vor allem von Kegeln, * 
Zylindern und Kugeln. Nach einigen Sätzen über 4 
algebraische Kurven und Flächen werden Dreh- 
flächen (insbesondere 2. Ordnung) besprochen. Die 
Durchdringungskurven zweier Drehflächen werden 
stets punkt- und tangentenweise konstruiert, wobei 4 
auch auf das Aufsuchen der Scheitelkrümmungs- 
kreise Wert ‚gelegt wird. Dann folgt ein Abschnitt 
über Schraublinien und Schraubflächen, wobei der von 3 
Th. Schmid eingeführte Begriff der Drehflucht ver- 7 
wendet wird. Eine Zusammenstellung der benötigten 4 
Sätze aus der Flächentheorie beschließt, wie bereits 
erwähnt, dieses ganz ausgezeichnete Lehrbuch, das 
gewiß dem Unterrichtsfach ‚‚Darstellende Geometrie“ 
erneuten Auftrieb geben wird. Es ist nur zu hoffen, 
daß diesem Band bald ein zweiter folgen wird, in dem 
die hier noch nicht behandelten Gebiete der Dar- 
stellenden Geometrie wie z.B. normale Axonometrie, 
Perspektive usw. gebracht werden und in dem sich 
auch Hinweise über den Stoff der konstruktiven Übun- 
gen, getrennt nach den einzelnen Fachrichtungen, 
befinden, die ja unentbehrliche Ergänzungen der 
Vorlesungen darstellen. 


Dresden R.B ereis 


Die besprochenen und angezeigten Bücher sind durch den 
Buchhandel zu beziehen. 


Nachrichten 
1. Die wissenschaftliche Jahrestagung 1959 der 
Gesellschaft für angewandte Mathematik und Me- 
chanik findet in der Zeit vom 19.—23. Mai in Han- 


nover unter der örtlichen Tagungsleitung von Herrn 
Prof. Dr. W, Quade statt. 


2. In der Zeit vom 15.—20. Juni 1959 veranstal- 
tet die UNESCO in Paris eine internationale Kon- 
ferenz über Informationsverarbeitung. 


3. An der Technischen Hochschule Dresden fand 
in der Zeit vom 5.—15. Januar 1959 im Institut für 
Maschinelle Rechentechnik unter Leitung von Prof, 
Dr.-Ing. Lehmann für die mathematischen Hoch- 
schullehrkräfte der DDR ein Lehrgang über Maschi- 
nelle Rechentechnik statt. Es schloß sich eine zwei- 
tägige Exkursion mit Besichtigung von einschlägigen 
Betrieben und Anlagen in Karl-Marx-Stadt, Jena und 
Ilmenau an. 


Verlag: Akademie-Verlag 
: P und Verlagsnummer dieses Heftes: 
erscheint monatlich. Bezugspreis: vierteljährlich DM 15,—. Zu- 


Abbestellungen können nur bis 4 Wochen vor Quartalsende anerka i 
das folgende Quartal noch geliefert. Veröffentlicht unter Lizenznummer ZLN 5011 des ee en a 


ns für Kultur, Hauptverwaltung Verlags- 


wesen. Gesamtherstellung: VEB Druckerei ‚Thomas Müntzer‘‘ Bad Langensalza (V/12/6) (1). Printed in Germany, 


